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Ice  Formation  In  Frequently  Transited  Navigation  Channels 


ROBERT  ETTEMA  AND  HUNG-PIN  HUANG 


INTRODUCTION 

Among  ice  navigators  an  old  adage  runs  "an 
ice-breaker  is  an  ice-maker"  (Keating  1970).  What 
is  implied  is  that  a  vessel  moving  through  an  ice- 
covered  channel  creates  an  open  water  track  over 
which  ice  may  regrow .  If  an  ice-covered  channel  is 
transited  many  times  along  more  or  less  the  same 
track  so  that  numerous  cycles  of  icebreaking  and 
regrowth  occur,  then  considerably  more  ice  may 
grow  in  the  channel  than  if  it  were  not  transited. 

The  effects  of  frequent  vessel  transit  on  ice- 
cover  formation  pose  important  concerns  for  winter 
navigation  of  heavily  trafficked  river,  harbor  and 
coastal  channels.  Of  prime  concern  is  the  likeli¬ 
hood  that  frequent  transits  may  cause  the  channels 
to  become  clogged  and  frozen  over  with  thick 
accumulations  of  broken  ice,  usually  termed  brash 
ice,  possibly  until  they  become  unnavigable. 
Another  concern  is  that  transits  may  prevent  or 
disrupt  the  initial  establishment  of  an  ice  cover 
over  a  river,  with  potentially  troublesome  conse¬ 
quences  if  the  river  water  becomes  supercooled 
and  generates  large  amounts  of  frazil  ice.  Also, 
when  an  ice  cover  breaks  up,  the  increased  volume 
of  ice  grown  as  a  result  of  transiting  may  increase 
the  potential  for  ice  jamming.  Increased  ice  vol¬ 
ume  may  also  exacerbate  difficulties  associated 
with  passing  broken  ice  at  dams. 

This  study  examines  the  effects  of  frequent 
vessel  transit  on  ice  formation  in  navigation  chan¬ 
nels.  It  was  conducted  under  the  U.S.  Army  Corps 
of  Engineers'  (COE)  River  Ice  Management  (RIM) 
program.  The  RIM  program  is  aimed  at  develop¬ 
ing  engineering  methods  for  managing  navigable 
rivers  that  become  ice  covered  during  winter.  In 
accordance  with  this  aim,  this  study  places  special 
focus  on  transits  of  river  channels  by  tow  vessels, 
such  as  those  that  ply  navigable  channels  in  the 
upper  basin  of  the  Mississippi  River.  Tow  vessels 
typically  comprise  a  towboat  and  one  or  more  flat- 
bottomed  barges  (Fig.  1).  The  results  of  the  study 


are  sufficiently  generic,  however,  so  as  to  apply  to 
ice-covered  channels  in  general. 

The  study's  principal  contributions  include 
detailed  descriptions  and  a  predictive  model  of 
the  effects  of  frequent  vessel  transit  on  ice-cover 
formation  over  navigable  channels.  The  descrip¬ 
tions  are  founded  on  computer-based  numerical 
modeling  and  the  results  of  extensive  laboratory 
experiments  conducted  in  an  ice  tank  with  scale- 
model  hulls  used  to  simulate  tow-barge  and  ship- 
form  hulls.  The  numerical  model  can  be  used  for 
predicting  ice  formation  in  channels  under  ranges 
of  meteorological  conditions  and  transit  frequen¬ 
cies.  It  can  give  information  on  volumes  of  ice 
grown  as  well  as  on  distributions  of  broken  ice  at 
typical  channel  cross  sections.  Its  development 
and  calibration  were  guided  by  the  results  of  the 
ice-tank  experiments.  Though  prior  studies  have 
used  numerical  modeling  of  ice  formation,  to  the 
best  of  the  authors'  knowledge,  this  is  the  first  for 
which  frequent-transit  effects  were  studied  using 
an  ice  tank.  The  numerical  modeling  and  ice-tank 
experiments  were  also  conducted  to  determine  if 
there  are  optimum  schedules  of  vessel  transits  that 
would  reduce  or  ameliorate  some  of  the  problems 
attributable  to  frequent  transiting. 

A  troublesome  feature  of  transiting  to  ws  is  their 
propensity  to  accumulate  broken  ice  beneath  their 
flat-bottomed  hulls.  Such  accumulationsmaysig- 
nificantly  increase  transiting  resistance  and,  for 
shallow  channels,  may  cause  tows  to  ground.  Also, 
when  moving  through  thinly  refrozen  brash  ice  or 
brash  that  is  not  refrozen,  tows  can  develop  and 
shove  extensive  false  bows  of  broken  ice.  This 
study  also  investigated  the  mechanics  of  ice  accu¬ 
mulation  beneath  and  ahead  of  tows. 


LITERATURE  REVIEW 

The  literature  on  ice  formation  in  navigation 
channels  is  not  extensive,  though  recently  it  has 


Figure  1.  Tow  moving  through  a  brash-ice<overed  channel  in  the  Mississippi 
River. 


burgeoned  in  response  to  heightened  interests  in 
developing  or  improving  the  use  of  navigation 
channels  during  frigid  winter  conditions.  The 
potential  concerns  incurred  with  vessel  transit  of 
ice-covered  channels  are  fairly  well  recognized, 
but  few  in-depth  studies  have  been  conducted  to 
address  them.  Especially  lacking  are  adequate 
descriptions  and  knowledge  of  the  ways  that  fre¬ 
quency  of  transiting,  hull  forms  of  vessels,  vessel 
maneuvering  activities,  channel  geometry  and 
water  currents  affect  ice-cover  formation. 

This  review  gathers  the  sparse  information  on 
ice  formation  in  navigation  channels.  It  begins 
with  the  standard  formulation  of  ice-cover  growth 
leading  to  the  well-known  Stefan  equation,  which 
is  the  basis  of  prior  predictive  models  of  ice  forma¬ 
tion  in  frequently  transited  channels.  A  digest  of 
field  observations  on  ice  formation  in  navigation 
channels  follows  for  both  river  channels  and  har¬ 
bor  and  coastal  channels.  The  review  also  exam¬ 
ines  the  formulations  of  prior  predictive  models  of 
ice  formation  in  vessel  tracks  opened  through  ice- 
covered  channels. 

Useful  assessments  of  the  engineering  concerns 
faced  in  developing  and  maintaining  navigation 
of  ice-covered  channels  are  given  by  Ashton  (1974), 
Mellor  et  al.  (1978),  Vance  (1980),  Eranti  et  al. 
(1983a,b),  Beurket  and  Argiroff  (1984),  Tronin  et 
al.  (1984),  Argiroff  and  Weigum  (1986)  and  Solve 
(1986).  These  publications  are  not  reviewed  here, 
although  some  are  referred  to  when  they  provide 
insight  on  how  vessel  transits  affect  ice  formation. 
Also  not  reviewed  are  publications  dealing  with 


the  strength  properties  of  brash  ice  or  with  the 
resistance  encountered  by  vessels  transiting  chan¬ 
nels  filled  with  brash  or  broken  ice.  These  topics 
are  germane  to  vessel  transiting  of  navigation 
channels,  but  they  exceed  the  scope  of  the  present 
study.  Appendix  A,  however,  does  deal  briefly 
with  the  resistance  encountered  by  tows  transiting 
layers  of  brash  ice.  Brash-ice  strength  properties 
are  discussed  by  Tatinclaux  and  Cheng  (1978), 
Mellor  (1980)  and  Urroz  and  Ettema  (1987).  Fur¬ 
ther  information  on  brash-ice  resistance  is  avail¬ 
able  from  Voelker  and  Levine  (1972),  Ashton  et  al. 
(1973),  Kotras  et  al.  (19 77),  Mellor  (1980),  Eskola 
(1983),  Ettema  et  al.  (1985)  and  Foreman  and 
Sandkvist  (1986). 

Ice-cover  growth 

The  rate  at  which  an  ice  cover  thickens  ther¬ 
mally,  or  "statically,"  as  an  ice  sheet  is  determined 
mainly  by  the  net  flux  of  heat  through  it,  as  is 
indicated  in  Figure  2.  Formulation  of  thickening 
(or  growth)  rate  usually  involves  the  following 
assumptions: 

•  The  temperature  profile  through  the  ice 
cover  is  linear; 

•  The  heat  fluxes  of  water  and  ice  are  negli¬ 
gibly  small  compared  with  heat  fluxes 
through  solid  ice  covers;  and 

•  The  heat-transfer  coefficient  between  ice 
and  air  h.  is  large  relative  to  the  ice  conduc¬ 
tivity  per  unit  ice  thickness  fc./rj. 

In  accordance  with  these  assumptions,  the  heat 
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Figure  2.  Energy  fluxes  associated  with  the  thermal  growth  of  ice  covers. 


balance  at  the  interface  between  an  ice  cover  and 
water  can  be  expressed  as 


4>i  -4>wi  =  pi  Li  ^ 


(1) 


where  <t>wi  =  thermal-energy  flux  from  water  to  ice 
(Jr  =  flux  by  conduction  through  the  ice 
cover 

p.  =  density  of  ice 


Also, 


L.  =  latent  heat  of  ice  fusion 
q  =  ice  sheet  thickness  at  time  t. 


<t>i=-fcj 


<TS-Tf) 


(2) 


where  k.  =  thermal  conductivity  of  the  ice  cover 
r  =  surface  temperature  of  the  ice  cover 
Tf  =  freezing  temperature  of  water 
(0°C  or  32°F,  for  pure  water). 

If  the  insulating  effects  of  snow  precipitation  on 
ice  growth  can  be  neglected,  is  usually  equated 
to  the  flux  of  thermal  energy  from  ice  to  air  <j>ia: 


<t>i  =<t>ia  =*ia  (Ts-Ta) 


(3) 


where  h  As  the  coefficient  of  heat  transfer  between 
ice  and  air,  and  Tais  the  air  temperature.  Combin¬ 
ing  eq  1,  2  and  3  and  eliminating  T  yields 

T-T 


-  f 


0  .  -<t>  •  = 

1  W1  ^ 


(4) 


+  -x 
k.  h. 


If  is  negligibly  small,  the  following  relation¬ 

ship  is  arrived  at  for  ice  cover  thickening: 


<t>i  = 


Tf  Ta  _  ,  dx\ 

i  p'  1  it 


(5) 


■  + 
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Integration  of  eq  5  gives 
2k.  2k 
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0  (6) 


where  C  is  a  constant  of  integration.  If  q  =  0  at  time 
f  =  0,  eq  6  is  reduced  to 


k j 
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(7) 


where  a  is  a  growth  coefficient  [sometimes  called 
air  freezing  index  (Lunardini  1981)],  which  is 
usually  less  than  its  face  value  of  (21c/ pi..) 0  5  =  34 
mm/(°C-days)0,5;  and  S  Sd  is  the  accumulated 
degree-days  of  air  temperature  below  the  tem¬ 
perature  at  which  water  freezes. 
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eq  7  is  reduced  to  the  well-known  Stefan  equation: 

I  Ni  \0'5 

n=«(lSd)  .  (9) 

Equation  9  is  the  basis  for  all  prior  predictors  of  ice 
growth  in  navigation  channels. 

Ashton  (1986)  and  Hausser  and  Parkinson  (1986) 
pointed  out  that  the  presence  of  a  "thermal  bound¬ 
ary  layer"  of  air  above  a  thin  ice  cover  should  be 
considered  in  the  foregoing  formulation.  This  re¬ 
finement  is  not  widely  used,  and  it  diminishes  in 
importance  once  an  ice  cover  is  more  than  a  few 
centimeters  thick.  The  usual  practice  is  to  assign 
values  that  have  been  proven  representative  of 
particular  conditions  of  ice  growth.  Michel  (1971 ), 
for  example,  recommended  using  an  a  value  of  28 
mm  /  (°C-days)°  5  for  windy  lakes  without  a  snow 
cover  (the  thermal  boundary  layer  is  thin).  For 
average  rivers  with  snow  and  sheltered  small  riv¬ 
ers  (the  thermal  boundary  layer  is  thick),  he 
suggested  a  values  of  17-13  mm/(°C-days)°  5  and 
14-7  mm/(°C-days)05,  respectively. 

Ice  formation  in  navigation  channels 

Ice  covers  over  frequently  transited  navigation 
channels  do  not  form  or  remain  as  monolithic 
covers.  Instead,  they  become  riven  with  the  tracks 
of  vessels,  sometimes  to  the  extent  that  they  con¬ 
sist  only  of  broken  ice,  which  is  usually  in  the 
brash-ice  size  range.  Vessel  tracks,  and  occasion¬ 
ally  entire  channels,  may  fill  with  brash  ice,  which 


consists  of  ice  fragments  that  typically  are  about 
0.02-2  m  in  major  dimension  (Mellor  1980,  Welsh 
and  Kingsbury  1975).  Brash  ice  normally  results 
when  broken  ice,  often  originating  as  broken  sheet 
ice,  is  further  abraded  and  fragmented.  Few  pro¬ 
files  have  been  measured  of  brash-ice  accumula¬ 
tions  in  navigation  channels.  Fewer  measurements 
have  been  made  of  the  size  distribution  and  micro¬ 
morphology  of  brash  ice  in  vessel  tracks  and  navi¬ 
gation  channels.  What  is  clear,  however,  is  that,  in 
addition  to  meteorological  conditions,  ice  forma¬ 
tion  in  navigation  channels  is  influenced  by  the 
hull  geometry  and  speed  of  transiting  vessels,  the 
frequency  of  transits,  the  channel  geometry  and 
water  currents.  None  of  these  influences  appears 
to  have  been  examined  extensively. 

For  channels  of  quiescent  water,  initial  ice  cov¬ 
ers  may  form  as  level  ice  sheets  that  thicken  ther¬ 
mally.  For  channels  with  significant  water  cur¬ 
rents,  ice  covers  may  establish  initially  as  accumu¬ 
lations  of  drifting  ice  floes,  broken  ice  or  frazil  ice 
(Ashton  1986)  that  consolidate,  then  thicken  ther¬ 
mally.  The  initial  establishment  of  an  ice  cover 
over  a  heavily  trafficked  channel  may  be  accompa¬ 
nied  by  the  simultaneous  formation  of  a  track  or 
tracks  opened  by  transiting  vessels.  Vessel  tracks 
through  moderately  or  lightly  trafficked  channels 
may  form  through  already  established  ice  covers. 
In  either  case  the  track  width  and  the  amount  of 
broken  ice  left  in  it  depend  on  icebreaking  patterns 
created  by  the  transiting  vessel,  as  well  as  the 
speed  and  hull  form  of  the  vessel. 

Several  detailed  measurements  of  icebreaking 
patterns  and  ice-fragment  sizes  have  been  made 
for  vessels  moving  through  level  sheets  of  ice. 
Kashteljan  et  al.  (1968),  Lewis  and  Edwards  (1970), 
Enkvist  (1972)  and  others  have  published  descrip¬ 
tions  of  ship  tracks  through  ice  sheets.  Ashton  et  al. 
(1973)  described  tracks  typically  formed  by  tows 
navigating  ice-covered  rivers.  The  track  left  by  a 
vessel  moving  through  an  ice  sheet  is  slightly 
wider  than  the  vessel  and  for  most  ship-form  hulls 
may  have  cuspate  edges.  The  tracks  left  by  tows 
are  characteristically  relatively  straight-edged  and, 
for  thin  ice,  can  have  widths  equivalent  to  tow 
beams  (Ashton  et  al.  1973),  although  this  depends 
on  the  tow  configuration  and  the  occurrence  of 
wave-induced  icebreaking.  The  size  of  the  ice  frag- 
ments  left  in  a  vessel  track  depends  on  ice  thick¬ 
ness,  ice  strength  and  hull  speed,  as  well  as  hull 
bow  and  form. 

Generally  ice  fragments  increase  in  size  (plan 
dimensions)  with  increasing  thickness  and  strength 
of  the  ice  sheet  transited  by  vessels.  Somewhat 
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smaller  fragments  can  occur  for  greater  hull  speeds 
(e.g.  Enkvist  1972),  and  larger  fragments  occur  for 
larger  bow  flare  angles.  McKindra  and  Lutton 
(1981)  statistically  analyzed  measurements  of 
broken  ice  in  tracks  formed  by  an  icebreaking  tug 
during  a  single  transit  through  level  ice  sheets  and 
ice  covers  composed  of  refrozen  brash  ice  on  Lake 
Superior.  They  found  that  refrozen  brash  ice  broke 
into  pieces  about  30%  smaller  than  did  ice  sheets. 
The  mean  size  of  ice-sheet  pieces  decreased  with 
hull  speed,  but  the  mean  size  of  refrozen  brash-ice 
fragments  increased  slightly  with  hull  speed.  For 
both  types  of  ice  cover,  they  found  the  size  distri¬ 
butions  of  broken  ice  to  be  log  normal.  Tuovinen 
(1978)  reached  the  same  conclusion  during  a  study 
of  channels  cut  through  ice  sheets  on  the  Baltic  Sea . 

There  are  few  published  studies  dealing  with 
ice  formation  in  frequently  transited  navigation 
channels  or  vessel  tracks.  The  most  informative 
studies  are  those  by  Ashton  et  al.  (1973),  Kannari 
(1983)  and  Sandkvist  (1986).  However,  these  are 
limited  to  field  investigations  of  particular  chan¬ 
nels  traversed  at  irregular  frequencies  by  vessels 
of  quite  different  hull  forms.  The  field  study  re¬ 
ported  by  Sandkvist  constitutes  perhaps  the  most 
detailed  set  of  field  observations  on  ice  formation 
in  a  frequently  transited  navigation  channel.  His 
study  was  conducted  with  an  icebreaker  that,  at 
irregular  frequencies,  plied  two  channels  at  Lulea 
Harbor  in  coastal  northern  Sweden.  Kannari's 
study  involved  several  icebreakers  and  ice-wor¬ 
thy  cargo  ships  that  transited  coastal  channels  in 
Finnish  waters  off  the  Gulf  of  Bothnia.  Ashton  et 
al.  gave  an  extensive  account  of  tow  transits  through 
diverse  conditions  of  ice-covered  channels  in  the 
upper  Mississippi  River. 

Significant  differences  in  the  observations  re¬ 
ported  by  the  three  studies,  mainly  brought  about 
by  the  radically  different  hull  forms  and  channel 
conditions  that  each  involved,  suggest  that  river 
channels  should  be  viewed  separately  from  har¬ 
bor  and  coastal  channels.  River  channels  are  char¬ 
acteristically  shallow  and  often  narrow.  They 
contain  flowing  water,  are  prone  to  becoming 
supercooled  and  generating  frazil  ice,  and  may  be 
punctuated  by  barriers  such  as  locks  and  dams. 
Furthermore,  they  nearly  always  are  transited  by 
shallow-draft  vessels  that  are  relatively  bluff- 
bowed  and  flat-bottomed,  such  as  tows,  barges, 
scows  and,  for  the  St.  Lawrence  and  St.  Marys 
rivers.  Great  Lakes  bulk  carriers.  Harbor  and  coastal 
channels,  on  the  other  hand,  are  usually  (though 
not  exclusively)  comparatively  deep,  adjoin  large 
water  bodies  that  are  less  readily  supercooled  than 


are  rivers,  and  are  transited  by  relatively  deep- 
draft,  ship-form  hulls.  There  is,  of  course,  some 
overlap  of  the  two  categories  of  channel,  as  they 
have  many  common  features  of  ice  formation.  On 
occasion,  barges  are  used  for  coastal  transporta¬ 
tion,  and  sea-going  ships  may  transit  portions  of 
the  St.  Lawrence  River  (Canada),  the  Saimaa  Canal 
(Finland)  or  the  Trollhatte  Canal  (Sweden). 

Transit  of  river  channels 

The  only  substantial  account  of  ice  formation  in 
river  channels  transited  by  tows  is  given  by  Ash¬ 
ton  et  al.  (1973).  Their  description  of  the  effects  of 
repeated  tow  transiting  on  ice  formation  is  compli¬ 
cated  by  the  fact  that  their  observations  were  made 
during  comparatively  few  transits  by  eight  differ¬ 
ent  configurations  of  tow.  Nonetheless,  their  study 
indicates  that,  as  the  result  of  water  current  and 
tow  entrapment  and  shoving  of  brash  ice,  brash  ice 
accumulates  very  nonuniformly  in  tracks  made  by 
tows.  They  observed  that  brash  ice  in  tracks  tended 
to  consolidate  at  certain  locations  separated  by 
extensive  areas  of  open  water,  as  illustrated  in 
Figure  3.  It  is  at  these  locations,  which  often  seem 
to  coincide  with  a  channel  constriction,  that  tows 
had  great  difficulty  in  transiting  thick  accumula¬ 
tions  of  brash  ice.  Ashton  et  al.  (1973)  found  that 
thick  accumulations  of  brash  ice  almost  halted  tow 
transit  at  two  especially  difficult  locations  (Fig.  4). 

Ashton  et  al.  (1973)  described  the  tracks  left  by 
tows  transiting  ice  sheets,  which  were  about 
0.20-0.25  m  thick,  as  being  "cleanly  edged"  and 
about  the  same  width  as  the  tow.  Much  of  the 
broken  or  brash  ice  created  during  the  transit  was 
either  cast  to  the  sides  of  the  tow,  with  very  little 
remaining  in  the  track,  or  accumulated  beneath 
barges.  They  detected  this  feature  from  an  exten¬ 
sive  trail  of  brash  ice  that,  on  one  occasion,  was  left 
behind  a  tow  once  it  had  entered  open  water.  They 
noted  that,  once  an  accumulation  was  in  place,  it 
required  increased  tow  speed  to  entrain  it. 

The  form  of  the  tracks  opened  by  tows  and  the 
amounts  of  brash  ice  accumulated  in  them  varied 
with  tow  arrangement,  tow  speed,  extent  of  ma¬ 
neuvering  and  ice-sheet  thickness.  Ashton  et  al. 
(1973)  reported  that  "while  an  initial  passage 
through  an  unbroken  ice  cover  may  result  in  a 
track  with  considerable  open  water  areas,  after 
several  passages  the  track  becomes  entirely  cov¬ 
ered  with  brash."  They  qualified  this  description 
by  noting  that  the  increased  areal  concentration  of 
brash  ice  was  due  in  great  part  to  the  additional 
breaking  of  ice  bordering  the  track  by  the  bow 
wave  produced  by  a  tow  during  a  subsequent 
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Figure  4.  Local  areas  of  thick  brash-ice  accumulations  in  the  Mississippi  River.  (From  Ashton  et  al. 
1973.) 


passage  of  the  track.  They  did  mention,  though, 
that  the  arrangement  involving  the  towboat  alone 
resulted  in  more  brash  ice  remaining  in  the  track; 
about  95%  of  the  track  was  covered  by  brash  ice. 
These  observations  su  ggest  that  for  tow  hulls,  hull 
length  influences  the  displacement  of  brash  ice 
from  a  track  and  therefore  affects  ice  formation  in 
a  frequently  transited  channel. 


Entrapment  and  shoving  of  broken  ice  seems  to 
be  a  characteristic  feature  of  tow  transit  through 
ice-covered  channels.  In  addition  to  collecting 
broken  ice  beneath  their  hulls,  tows  moving 
through  loose  (not  refrozen)  covers  of  brash  ice 
develop  false  bows  of  brash  ice.  False-bow  and 
under-hull  accumulations  of  ice  rubble  may  cause 
especially  severe  problems  for  tows  attempting  to 
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Figure  5.  Tow  shoving  brash  ice  while  approaching  a  lock. 


pass  through  locks.  Brash  ice  accumulated  ahead 
of  a  lead  barge  can  be  pushed  into  locks,  creating 
problems  for  operation  of  their  miter  gates.  Ice 
accumulated  beneath  barges  may  cause  them  to 
become  stuck  on  the  exit  sill  of  locks.  Figure  5 
illustrates  the  amount  of  ice  that  a  tow  can  shove 
ahead  of  itself,  especially  when  moving  through  a 
narrow  or  confined  channel.  Under-hull  accumu¬ 
lations  of  ice  not  only  increase  the  resistance  expe¬ 
rienced  by  tows,  but  they  may  also  cause  them  to 
ground  in  shallow  channels.  Reports  prepared  by 
the  U.S.  Army  Corps  of  Engineers  (1976, 1978)  and 


Borland  (1987)  cite  frequent  instances  of  tows 
grounding  on  ice  under  the  hull.  One  report  (U.S. 
Army  Corps  of  Engineers  1978)  mentions  an  inci¬ 
dent  in  which  a  tow  grounded  on  an  accumulation 
that  was  about  6.1  m  thick.  Another  (U.S.  Army 
Corps  of  Engineers  1976)  documents  an  incident 
in  which  a  tow  grounded  on  about  3.9  m  of  brash 
ice  when  attempting  to  pass  through  Lock  and 
Dam  No.  26  in  the  Mississippi  River.  Figure  6 
shows  ice  removed  from  beneath  a  tow  and  left  in 
a  lock  after  a  difficult  lockage.  Operators  of  tows  in 
thewaterwaysoftheupperMississippiRiver  report 


Figure  6.  Ice  removed  from  beneath  a  tow  as  it  left  a  lock. 


that  accumulation  problems  can  be  particularly 
severe  for  tows  with  empty,  open  barges  in  frigid 
conditions,  because  brash  ice  may  adfreeze  to 
them.  Incidents  of  grounding,  especially  in  locks, 
have  become  of  such  concern  to  the  Corps  of 
Engineers,  as  well  as  to  tow  operators,  that  at  Lock 
and  Dam  No.  26  the  Corps  of  Engineers  has  in¬ 
stalled  sonar  transducers  for  monitoring  ice  accu¬ 
mulation  beneath  tows  (Borland  1987). 

One  river  that  regularly  poses  problems  for 
winter  navigation  is  the  St.  Marys  River  linking 
Lakes  Superior  and  Huron.  Kotras  et  al.  (1977), 
Mellor  et  al.  (1978)  and  Vance  (1980)  described  the 
difficulties  experienced  by  Great  Lakes  bulk  carri¬ 
ers  in  transiting  relatively  thick  accumulations  of 
brash  ice  that  form  along  the  St.  Marys  River. 
Although  these  studies  did  not  explicitly  identify 
the  origin  of  the  brash-ice  accumulations,  they  did 
infer  that  the  accumulations  are  attributable  to  the 
water  current  and  the  shoving  of  ice  by  bluff- 
bowed  bulk  carriers.  Mellor  et  al.  ventured  some 
concepts  on  how  best  to  manage  brash-ice  accu¬ 
mulations,  as  did  Vance.  Scheduling  of  transits 
was  not  considered.  Kotras  et  al.  also  provided 
some  anecdotal  information  on  transits  of  brash 
ice  in  vessel  tracks  through  ice  covers  on  Lake 
Superior.  The  brash-ice  layers  were  thin  and  posed 
no  significant  problem  except  when  wind-driven 
ice  sheets  pressurized  the  layers. 

Transit  of  harbor  and  coastal  channels 

Field  studies  indicate  that  ice  formation  in  ship 
tracks  through  ice-covered  harbor  and  coastal 


channels  is  more  orderly  than  in  tracks  through 
ice-covered  river  channels  transited  by  tows.  The 
greater  orderliness  can  be  attributed  to  the  inabil¬ 
ity  of  typical  ship-hull  forms  to  entrap  and  shove 
broken  ice,  and  to  the  absence  of  significant  water 
currents,  which  would  cause  brash  ice  to  drift  and 
collect  as  locally  thick  accumulations. 

Detailed  field  studies  of  ice  formation  in  spe¬ 
cific  harbor  and  coastal  channels  have  been  re¬ 
ported  by  Sandkvist  (1978, 1980, 1981, 1982, 1986), 
Tuovinen  (1978)  and  Kannari  (1983).  These  studies 
show  that  frequent  transits  of  harbor  and  coastal 
channels  result  in  ship  tracks  that  fill  with  more- 
or-less  uniform  layers  of  brash  ice  flanked  by 
inverted  linear  mounds  or  ridges  of  brash  ice 
accumulated  beneath  the  ice-cover  edges  adjoin¬ 
ing  the  tracks.  Figure  7  depicts  a  typical  cross 
section  of  a  brash-filled  track.  Their  studies  indi¬ 
cate  that  layers  and  ridges  of  brash  ice  resulted 
largely  from  repeated  cycles  of  ice  regrowth  and 
icebreaking,  rather  than  from  brash  ice  conveyed 
and  collected  by  the  water  current  or  by  vessel 
shoving. 

Sandkvist,  who  studied  ice  formation  in  two 
200-m-long  channels  at  Lulea  Harbor,  Sweden, 
obtained  information  on  the  rates  at  which  layers 
of  brash  ice  formed  and  thickened  in  ship  tracks, 
their  cross-sectional  profiles,  size  distributions  of 
their  constituent  brash  ice,  and  layer  porosities.  He 
used  the  data  on  layer  thickening  (Fig.  8)  to  de¬ 
velop  a  fairly  simple  predictor  of  ice  formation. 
The  two  channels  were  ostensibly  transited  at 
different  frequencies,  but  there  were  no  major 
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Figure  7.  Typical  cross  section  of  a  vessel  track  filled  with  brash  ice. 
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Figure  8.  Field  data  on  brash  ice  and  ice-caver  thicknesses  in  ship  tracks  formed  in  Luled  Harbor,  Sweden  (Sandkvist  1986). 
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Figure  9.  Cross  section  of  a  navigation  channel  covered  with  brash  ice.  (From  Kannari  1983.) 


differences  in  ice  formation. 

Tuovinen  and  Kannari  studied  ice  formation  in 
coastal  channels  in  Finnish  waters  of  the  Gulf  of 
Bothnia.  They  characterized  the  dimensions  and 
mechanical  properties  (primarily  shear  strength 
and  friction)  of  brash-ice  layers  in  what  they  termed 
ice-clogged  ship  tracks.  The  brash-ice  accumula¬ 
tions  in  the  navigation  channels  monitored  by 
Sandkvist  and  Kannari  are  typified  by  the  cross- 
sectioned  profile  shown  in  Figure  7.  However, 
depending  on  the  relative  magnitudes  of  naviga¬ 
tion-channel  width  and  vessel  beam,  the  profiles 
may  be  more  irregular,  as  shown  in  Figure  9. 
Unsteady  speeds  or  maneuvering  may  cause  lon¬ 
gitudinal  profiles  of  brash-ice  accumulations  to  be 
somewhat  irregular  or  undular,  as  indicated  in 
Figure  10. 

With  each  transit  of  a  vessel  through  a  naviga¬ 
tion  channel,  brash  ice  is  displaced  sideways  be¬ 
neath  the  bordering  ice  cover.  There,  brash  ice 
accumulates  to  form  inverted,  linear  mounds  or 
ridges  that  run  more-or-less  parallel  to  the  track  so 
as  to  confine  brash  ice  in  the  track.  The  ridges 
solidify  due  to  freezing  from  above  and  local  freez¬ 


ing  at  contact  points  between  ice  pieces.  If  the 
channel  is  relatively  wide  compared  to  the  vessel 
beam  (maximum  width),  an  irregular  series  of 
furrows  may  form  between  the  ridges.  Brash-ice 
accumulations  may  be  particularly  irregular  if 
vessels  do  not  follow  exactly  the  same  track  or  if 
vessels  of  different  sizes  ply  the  channel. 

Detailed  measurements  of  the  sizes  of  individ¬ 
ual  brash-ice  pieces  led  Sandkvist  to  conclude  that 
the  distribution  of  brash-ice  thickness  was  not  log¬ 
normal,  as  had  been  found  for  the  distribution  of 
horizontal  lengths  of  freshly  broken  ice  pieces.  He 
also  deduced  that  additional  breaking  caused  the 
mean  plan  size  of  brash  ice  to  decrease  gradually 
with  increasing  number  of  transits  until  it  attained 
some  constant  value.  He  reported  further  that 
brash-ice  pieces  become  more  spherical  and  that 
the  average  thickness  or  vertical  dimension  in¬ 
creased  to  an  asymptotic  value,  called  the  equilib¬ 
rium  vertical  block  thickness  (Sandkvist  1982). 
Sandkvist's  measurements  of  brash-ice  fragment 
thickness  are  reproduced  in  Figure  11.  Tuovinen 
(1978),  Greisman  (1981)  and  Kannari  (1983)  also 
reported  that  brash  ice  in  ship  tracks  becomes 
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Figure  11.  Thickness  distribution  of  brash-ice  fragments.  (From  Sandkoist  1983.) 


rounded.  Greisman  coined  the  term  "cannonball 
ice"  for  such  ice.  Kannari  presented  data  on  meas¬ 
ured  sizes  of  brash  ice  in  channels  in  coastal  re¬ 
gions  of  theGulf  of  Bothnia  (Fig.  12).  He  concluded 
that  the  size  distribution  is  affected  by  XSd  and  by 
the  number  of  transits. 

Prior  formulations 

Five  semi-empirical  formulations  of  ice  forma¬ 
tion  in  ship  tracks  have  been  proposed  (Ashton 
1974,  Berenger  and  Michel  1975,  Vance  1980, 
Sandkvist  1981,  Hamza  1985).  Each  formulation  is 
based  on  eq  7  and  invokes  the  assumption  that  ice 
covers  regrow  over  a  porous  layer  of  uniformly 
distributed  brash  ice.  In  accordance  with  this  as¬ 


sumption,  each  formulation  yields  only  an  aver¬ 
age,  or  equivalent,  thickness  of  ice  accumulation 
in  ship  tracks. 

Ashton  (1974)  modified  eq  9  and  proposed  the 
following  expression  for  the  average  thickness  of 
ice  regrown  and  accumulated  in  an  ice-covered 
track: 


q=pa 


Ni 

X 


0.5 


+(1-P) 
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Figure  12.  Size  distribution  of  brash-ice  fragments.  (From  Kannari  1983.) 
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The  first  term  is  the  thickness  of  ice  grown  on  open 
water,  which  occupies  a  fraction  |3  of  the  track's 
surface  area;  the  second  term  is  the  equivalent 
thickness  of  the  ice  cover  (occupying  1-6  of  the 
track's  surface)  over  the  track  prior  to  the  i®  transit 
of  a  vessel  through  the  track;  and  the  third  term 
relates  to  the  additional  thickening  of  the  ice  frag¬ 
ments  floating  in  the  track.  Thickness  T).  is  actually 
a  length  index  that  is  equivalent  to  the  accumu¬ 
lated  volume  of  ice  grown  per  unit  surface  area  of 
track  prior  to  the  Ith  vessel  transit.  Ashton  suggested 
a  =  24  mm/(°C-day)0'5,  based  on  observations  of 
ice  growth  at  sections  of  the  Mississippi  River.  He 
also  suggests  using  (3  =  0.25,  at  least  for  initial 
transits.  Berenger  and  Michel  (1975)  developed  a 
considerably  more  complicated  expression  for  ice 
accumulation  in  a  track.  They  proposed  that  the 
increased  average  thickness  At|i  of  ice  formation 
between  transits  be  calculated  as 


+  C2  (1  -p)"1"1 

/  Ni  ]°-5  \ 

^  (""Is  +  ilaf+a2  !  Sd  -  (mTls  +Tla)  J 

(ID 

where  p  =  ratio  of  water  volume  to  total  volume 
(porosity) 

qs  =  thickness  of  individual  ice  fragments 
composing  a  brash-ice  layer 
Zw  =  distance  between  the  upper  ice  and  the 
water  surface 

T|a  =  additional  or  residual  thickness  of  ice 
formed  under  the  layer 
m  =  layer  number 
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C1  =p,  for  a  |  X  sdJ  >  frls  ~ZW> 


or 


C2  =  1,  for  m  S  K 


C2  =  0,  for  m  >  K 


K  =  maximum  number  of  layers  to  form  the 
residual  ice  {l/[p  (1-p)]  m}. 

In  eq  11,  the  first  term  is  ice  growth  over  open 
water  at  the  surface  of  the  track;  the  second  term  is 
the  ice  growth  within  the  pores  of  layers  of  ice 
fragments;  and  the  third  term  is  the  additional,  or 
residual,  ice  grown  beneath  a  layer  of  brash  ice. 
For  periods  of  regrowth  Pn,  or  M  (=1  n)  transits,  the 
total  volume  of  ice  grown  per  unit  plan  area  of  the 
track  is  estimated  from 


M 


(12) 
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A  difficulty  in  applying  Michel  and  Berenger's 
formulation  is  that  m  and  K  are  difficult  to  assess. 

Vance  (1980)  suggested  the  following  expres¬ 
sion  for  estimating  q.: 
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(13) 


The  first  two  terms  relate  to  continuous  thickening 
of  brash-ice  pieces;  the  third  term  is  ice  growth 
over  open  water.  A  comparison  of  eq  10  and  13 
shows  that  eq  13  is  misformulated  and  would  lead 
to  a  significant  overprediction  of  ice  regrowth.  In 
eq  13  the  summation  of  the  first  two  terms,  com¬ 
bined,  should  be  to  the  power  0.5,  as  is  the  case  for 
the  second  and  third  terms  in  eq  10.  Also,  the  third 
term  in  eq  13  should  be  multiplied  by  P,  as  is  the 
first  term  in  eq  10. 

Sandkvist  adapted  eq  9,  the  Stefan  equation, 
and  proposed  the  following  simple  formulation 


for  the  average  or  "equivalent"  thickness  of  ice 
formation  q.,  which  he  defined  as  the  resulting 
thickness  of  brash  ice  accumulated  in  a  vessel 
track,  of  width  equal  to  the  vessel  beam,  if  all  brash 
ice  is  contained  in  the  track,  as  depicted  in  Figure 
13: 


M  /  Ni  \°'5 


Tii=Ti0+I“  Xsd 

1=1  /= 1 


The  first  term  is  the  initial  ice-cover  thickness 
before  the  first  transit  (i = 0);  the  second  term  is  the 
equivalent  thickness  of  ice  formation  for  M  tran¬ 
sits  of  a  channel.  By  setting  a = 12  mm/(°C-days)0'5, 
Sandkvist  was  able  to  fit  eq  14  with  field  measure¬ 
ments  of  brash  ice  accumulation  in  two  test  chan¬ 


nels  at  Lulea  Harbor  in  northern  Sweden.  The 
result  is  shown  in  Figure  8.  Eranti  et  al.  (1983b) 
reported  using  eq  14,  with  a = 6.5  mm/(°C-days)05, 
in  estimating  ice  formation  in  a  30-m-wide  chan¬ 
nel  transited  by  12-m-beam  ships  (Saimaa  Canal, 
Finland). 

Hamza  (1985)  proposed  the  following  expres¬ 
sion  for  accumulated  brash-ice  thickness  Aq.  be¬ 
tween  two  successive  transits: 


%  =  «  I  sd, 
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Figure  13.  Sandkvist's  model  of  brash-ice  accumulation. 
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The  first  term  is  the  ice  growth  over  open  water, 
and  the  second  term  is  the  continuous  thickening 
of  ice  fragments  in  the  channel.  The  thickness  of 
accumulated  brash  ice  after  the  z'th  passage  is 


T,i=ifpK+ATii+- 


(16) 


Equation  15  is  faulty  and  leads  to  an  overestima¬ 
tion  of  ice  regrowth,  because  the  first  term  should 
be  multiplied  with  p,  and  the  second  term  with 
(1— P),  as  in  eq  10.  Equation  16  is  also  incorrect 
because  it  involves  (1— P)-1  as  a  coefficient  for  T|., 
which  has  involved  the  fraction  of  open  water  P  in 
eq  15.  Unlike  the  other  models  cited,  Hamza  con¬ 
sidered  that  p  (=  p)  is  not  constant  but  depends  on 
the  frequency  of  transit.  He  suggested  that  p  de¬ 
creases  with  the  increasing  frequencies  of  transit. 
The  accelerated  thickening  of  an  ice  cover  with  an 
underlayer  of  frazil  ice  is  similar  to  the  regrowth  of 
ice  in  a  navigation  channel  filled  with  brash  ice. 
Calkins  (1979)  modified  eq5  to  include  the  effect  of 
an  underlayer  of  frazil-ice  slush  on  a  statically 
thickening  ice  cover.  He  assumed  that  with  a  layer 
of  slush  present,  (J>wi  =  0,  and  proposed  allowing 
for  slush-layer  porosity  by  setting 

•  <17> 

dt 

Substituting  for  L.  from  eq  2  and  integrating  the 
resulting  expression  yields,  if  q  =  0  at  time  t  =  0, 
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A  problem  with  using  eq  18  for  estimating  ice 
growth  through  a  porous  layer  of  brash  ice  is  that 
values  of  a  would  have  to  be  different  than  those 
recommended  for  the  growth  of  smooth  mono¬ 
lithic  ice  covers.  The  geometry  of  jumbled  brash 
ice  would  call  into  question  the  use  of  eq  2  and  the 
assumptions  behind  it. 

A  shortcoming  in  the  foregoing  models  is  that 
they  do  not  take  into  account  the  influences  on  ice 
growth  of  brash-ice  displacement  during  transit¬ 
ing.  For  instance,  except  for  Hamza's  model,  the 
models  are  based  on  the  assumption  that  a  full 
layer  of  brash  ice  develops  immediately  across  a 
track  transited  by  vessels.  Furthermore,  it  is  as¬ 
sumed  for  all  models  that  all  ice  grown  remains  in 
vessel  tracks.  In  accordance  with  the  first 
assumption,  P  is  held  constant,  either  as  zero  or 


0.25  (Ashton).  Hamza,  however,  took  P  to  depend 
only  on  number  of  transits.  As  a  result  of  the 
second  assumptions,  none  of  the  models  enables 
prediction  of  ridge  thickness.  Although  Sandkvist's 
empirical  expression  gives  an  equivalent  thick¬ 
ness  of  the  brash-ice  layer,  the  other  models  give 
only  the  volume  of  solid  ice  produced.  They  do  not 
allow  the  prediction  of  the  bulked  thickness  of  a 
layer  that  is  brought  about  by  icebreaking  that 
converts  a  cover  of  solid  ice  to  a  porous  brash-ice 
layer. 

The  models  are  presented  as  being  capable  of 
predicting  ice  growth  in  ice-covered  channels 
generally.  However,  different  hull  forms  and  sizes 
may  displace  ice  differently,  result  in  different 
patterns  of  accumulation  and  thereby  influence 
ice  growth.  Also,  channel  geometry,  notably  width 
and  depth  relative  to  the  dimensions  of  transiting 
vessels,  may  influence  ice  formation,  as  may  water 
currents. 


FORMULATION  OF  A  PREDICTIVE  MODEL 

It  is  likely  that,  at  best,  deterministic  numerical 
models  can  simulated  ice  formation  in  frequently 
transited  navigation  channels  only  approximately. 
The  nature  of  transiting,  icebreaking  and  redis¬ 
tributing  broken  and  brash  ice  entails  stochastic 
processes  that  are  difficult  to  predict  and  may  vary 
from  vessel  to  vessel  and  channel  to  channel. 
Nonetheless,  a  comparatively  simple  numerical 
model  is  formulated  here  as  a  means  of  evaluating 
the  sensitivity  of  ice  formation  to  transit  frequency 
and  number,  air  temperature,  hull  form  and  chan¬ 
nel  geometry.  The  model  is  an  attempt  to  simulate 
ice  formation  more  realistically  than  do  the  prior 
formulations,  by  taking  into  account  the  effects  on 
ice  growth  of  brash-ice  redistribution  during  each 
transit  of  a  vessel.  Its  formulation  was  guided  in 
part  by  the  observations  of  brash-ice  accumulation 
and  from  insights  and  data  gained  from  the  ice- 
tank  experiments. 

The  predictive  model  is  designed  to  simulate 
ice  formation,  or  cycles  of  ice  regrowth  and  brash- 
ice  accumulation,  in  a  frequently  transited  naviga¬ 
tion  channel  such  as  is  illustrated  in  Figure  14.  An 
assumption  that  is  fundamental  to  the  model  is 
that  transiting  vessels  displace  brash  ice  only  side¬ 
ways  and  do  not  shove  significant  amounts  of  it 
along  their  tracks.  Also,  it  is  assumed  that  no 
fragments  break  from  the  ice  sheet  bordering  the 
track.  In  accordance  with  these  assumptions  and 
Figure  14,  ice  formation  can  be  described  as  fol- 
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lows.  The  first  few  transits  of  an  ice-covered  navi¬ 
gation  channel  leave  an  open  layer  of  broken  or 
brash  ice  uniformly  spread  over  a  track  opened 
through  it,  with  some  ice  being  cast  beneath  the 
adjoining  ice  cover  along  each  side  of  the  track 
(Fig.  14b).  With  subsequent  transits  and  periods  of 
ice  regrowth,  more  ice  is  broken,  causing  the  track 
to  become  increasingly  covered  by  a  thickening 
layer  of  brash  ice.  Brash  ice  cast  laterally  collects  as 
inverted  ridges  flanking  the  channel  (Fig.  14c). 
Once  the  track  is  completely  covered  with  a  layer 
of  brash  ice,  further  transits  cause  the  brash-ice 
layer  to  thicken  uniformly  but  lagging  behind  the 
thickening  of  the  ice  ridges  because  brash  ice  is 
always  displaced  sideways  to  ridges  beneath  the 
ice  cover  bordering  the  track  (Fig.  14d).  Also,  each 
transit  causes  the  entire  brash-ice  layer  to  be  dis¬ 
rupted  and  churned  over  such  that  refreezing 
commences  again.  However,  some  brief  period  of 
time  is  required  for  the  water  or  the  water  and  ice 
in  the  track  to  calm  sufficiently  for  ice  to  regrow.  In 
relatively  shallow  channels  the  ridges  may  even¬ 
tually  extend  down  to  the  channel  bed  (Fig.  14e), 
confining  brash  ice  within  the  track. 

The  predictive  model 

Consider  a  section  of  vessel  track  filled  with 
brash  ice,  as  shown  in  Figure  15.  Conservation  of 
thermal  energy  flux  between  ice,  water  and  frigid 
air  enables  the  volume  of  ice  in  the  track  per  unit 
area  of  track  surface  at  any  instant  q.  to  be  ex¬ 
pressed  as 
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The  first  term  is  the  value  of  q;  immediately  prior 
to  the  ;fh  vessel  transit;  the  second  term  is  the 
average  thickness  of  ice  regrowth  over  open  wa¬ 
ter;  and  the  third  term  is  the  equivalent  thickness 
of  ice  grown  through  and  under  accumulated 
brash  ice  occupying  the  track  at  areal  concentra¬ 
tion  (1-p).  In  the  second  and  third  terms,  the 
conventional  cumulative  degree-days  of  freezing 
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is  modified  to  allow  for  the  time  required  for  water 
in  the  track  to  become  sufficiently  calmed  after  a 
vessel  transit  such  that  ice  regrowth  can  com¬ 
mence.  In  accordance  with  this  condition, 
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Figure  15.  Modeled  geometry  of  brash-ice  accumulation  across  a  vessel  track. 
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is  modified  to 
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with  P,  =  1  //t  being  the  period  between  transits, 
and  t  is  the  time  for  the  water  or  brash-ice  layer  to 
calmDefore  regrowth  commences  after  a  transit. 
At  very  large  frequencies  of  transiting,  when  P  -t 
approaches  zero,  little  opportunity  occurs  for 
regrowth,  and  the  effect  of  transiting  is  to  maintain 
an  open  track.  During  each  transit  of  the  track, 
when  it  is  fully  covered  by  brash  ice  as  depicted  in 
Figure  14,  the  regrown  ice  cover  is  broken,  mixed 
with  existing  brash  ice,  and  bulked  with  water  to 
form  a  layer  of  uniform  thickness: 

= 

(1-£i)hbi--(1-P)(1-Pi)  +  K-T^i-ll] 

(20) 


in  which  e.  is  the  rate  of  brash  or  broken  ice 
displacement  from  the  track  to  the  ridges  flanking 
the  track. 

Evaluation  of  ice  volumes  grown  within  a  po¬ 
rous  accumulation  of  brash  ice  entails  a  calcula¬ 
tion  that  agrees  with  many  assumptions  concern¬ 
ing  the  thermodynamics  of  ice  growth  within  an 
assemblage  of  irregular  cavities  (namely,  eq  11, 
which  was  proposed  by  Berenger  and  Michel). 
Rather  than  becoming  enmeshed  in  a  net  of  tenu¬ 
ous  assumptions  and  calculational  complexities, 
such  as  ascertaining  when  growth  through  a  layer 
reaches  the  layer  underside  and  thermal  thicken¬ 
ing  of  layer  commences,  the  present  formulation 
invokes  the  following  simple  assumption:  for  the 
purpose  of  estimating  ice  growth  associated  with 
a  layer  of  brash  ice,  it  is  assumed  that,  between 
transits,  ice  growth  through  and  beneath  a  layer  of 
thickness  qbi  and  porosity  p  can  be  reasonably 
represented  by  the  growth  of  an  equivalent  mono¬ 
lithic  ice  cover  of  thickness  (l-p)qb.. 

So  far,  no  account  has  been  macle  of  brash  ice 
that  may  have  been  displaced  to  ridges  develop¬ 
ing  beneath  the  ice  cover  adjoining  the  track.  A 
simple  mass-conservation  relationship  can  be  for¬ 
mulated  to  account  for  brash  ice  that  accumulates 
in  the  track  and  collects  in  the  ridges  flanking  it. 
The  volume  of  ice  displaced  per  unit  length  of 
track  surface,  and  supplied  to  each  ridge,  is 


Each  ridge  is  assumed  to  thicken  and  broaden, 
as  shown  in  Figure  14,  in  a  manner  similar  to 
heaping  a  linear  pile  of  granular  material.  If  at  any 
instant  the  ridge  thickness  is  r|ri,  and  its  shape  is 
defined  in  terms  of  an  angle  of  static  repose  of 
brash  ice  0^  the  ridge  volume  is 


tan  <j>r  + 


2(! 


The  first  term  is  the  volume  of  ice  (per  unit  length 
of  track)  already  included  in  a  ridge;  the  second 
term  is  the  volume  of  brash  ice  (per  unit  length  of 
track!'  displaced  from  the  track;  and  the  last  term  is 
the  volume  of  ice  grown  (per  unit  length  of  track) 
at  each  ridge.  As  the  last  quantity  is  negligible 
comp  ;K  a  io  the  other  two,  it  is  discarded  hereaf- 
te, .  The  ridge  thickness  qri  and  width  VV.  after  the 
Ith  transit  can  be  expressed  as 


and 

W  .  =2tl  .  tan  0  . 

n  'n  r 


(22) 

(23) 


The  maximum  values  of  ririand  r^are,  of  course, 
limited  by  channel  depth  and  water  currents.  Brash 
ice  resting  along  a  ridge  may  be  entrained  by  water 
currents. 

In  the  formulation  of  ice  formation  in  transited 
channels,  several  quantities  are  difficult  to  know. 
These  include  the  areal  concentration  of  the  open 
water  along  a  track  (3,  the  porosity  of  the  brash-ice 
layer  p,  the  proportion  of  ice  displaced  from  the 
layer  to  the  ridges  e,  the  slope  angle  of  the  ridge 
accumulations  0^  and  the  time  required  for  ice 
regrowth  to  commence  after  a  transit  t  .  Among 
other  things,  (3,  e,  p,  0r  and  t  are  influenced  by  the 
hull  geometry,  draft  and  speed,  the  number  and 
frequency  of  transits,  the  channel  geometry  and 
water  currents,  and,  in  some  cases,  the  relative 
sizes  of  ice  fragments  and  vessel  hull.  The  value  of 
P  diminishes  with  increasing  number  of  vessel 
transits,  becoming  more-or-less  constant  once  a 
track  is  completely  covered  by  brash  ice.  The 
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Figure  16.  General  transit  schedule. 


proportion  of  ice  displacement  e,  in  addition  to 
varying  with  hull  geometry,  draft  and  speed  is 
influenced  by  the  ratio  of  ridge  thickness  and 
channel  thickness  Tir/t|b.  The  larger  T(r/r|fcf  is,  the 
less  ice  is  displaced  and  accordingly  the  smaller  is 
e.In  a  sense,  during  a  series  of  hull  transits,  T]r/rib 
and  e  are  cyclic.  A  goal  for  the  ice-tank  experi¬ 
ments  was  to  determine  how  |3,  e,  p,  0rand  t  vary 
with  some  of  these  factors. 

Transit  frequencies 

Vessels  may  transit  channels  at  a  single  fre¬ 
quency.  For  example,  a  channel  may  be  transited 
round-the-clock  at  hourly  or  multi  -hourly  frequen¬ 
cies.  Typically,  however,  transits  conform  to  some 
form  of  navigation  schedule.  The  simplest  are 
daily  transits  limited  to  daylight  hours.  Transits 
may  also  occur  randomly. 

Most  schedules  can  be  described  in  terms  of  a 
general  schedule  (Fig.  16).  After  an  initial  period  of 
ice-cover  formation  P.,  a  sequence  ensues  of  n 
transits  with  a  time  interval  P  for  a  sequence 
period  Pg.  If  n  =  1,  then  P{=  PJ  n;  if  n  >  1,  then  P( 


-PJ  (n-1 ).  The  sequences  of  transits  may  be  inter¬ 
mittent,  being  separated  by  a  period  of  no  navi¬ 
gation  P  .  A  fairly  simple  schedule  that  can  be 
defined  in  terms  of  three  periods,  Pv  Pn  and  Pq, 
gives  rise  to  three  transit  frequencies: 

•  Frequency  /  =  1/P(  defines  the  frequency  of 
transits  within  a  sequence  of  transits  (e.g. 
one  transit  every  two  hours); 

•  Frequency /  =  n/(P  +Pq)  =  «/Pn  defines  the 
frequency  of  transits  within  a  navigation 
period  (e.g.  six  transits  every  day);  and 

•  Frequency  /  =  1/Pn  -  fjn  defines  the  se¬ 
quences  witftin  a  navigation  period  (e.g.  one 
sequence  of  transits  daily). 

For  example,  consider  a  transit  schedule  requir¬ 
ing  20  days  of  daily  transits  with  five  transits 
during  eight  hours  of  daylight.  For  this  schedule, 
n  =  5,  ft  =  1  transit/(2  hours),  Pg  =  8  hours  of 
daylight,  Pq=  16  hours  of  darkness,  fs  =  5/ day,  Pn 
=  1  day,  and  /  =  1  sequence/day.  Another  sched¬ 
ule  may  entail  continuous  sequences  of  four  ves¬ 
sels,  0.5-hour  apart,  with  six  hours  between  se¬ 
quences;  i.e.  n  =  4,  /  =  2  transits/ hour,  Pg  =  1.5 
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hours,  Pq  =  6  hours,  /  =  4/(75  hours),  and  /  =  1 
sequence/ (7.5  hours).  Theoretically  the  general 
schedule  shown  in  Figure  14  is  part  of  a  larger,  self- 
similar  array  ‘hat  can  be  described  more  formally 
in  terms  of  fractal  geometry.  Typical  frequencies 
of  transiting  on  ice-covered  rivers  within  the  U.S. 
range  from  1  / (several  days)  to  about  6/day.* 

The  predictive  model  is,  after  being  calibrated 
against  the  laboratory  results,  exercised  later  with 
the  aims  of  determining  the  effects  of  frequent 
vessel  transit  on  river  ice  formation,  and  ascertain¬ 
ing  if  scheduling  of  vessel  transits  might  be  a 
feasible  means  for  minimizing  excessive  ice  growth 
in  navigation  channels.  The  exercises  involve  vessel 
transits  during  a  range  of  meteorological  condi¬ 
tions  and  channel  geometries. 

ICE-TANK  EXPERIMENTS 

The  ice  tank  at  the  Iowa  Institute  of  Hydraulic 
Research  (IIHR)  was  used  to  simulate  an  ice-cov¬ 
ered  navigation  channel  during  experiments  aimed 
at  determining  and  documenting  the  effects  of 
freque  t  vessel  transit  on  ice-cover  formation.  The 
ice  tank,  which  is  located  in  a  refrigerated  labora¬ 
tory,  enabled  these  effects  to  be  examined  under 
controlled  conditions  of  air  temperature,  channel 
depth,  hull  geometry  and  transit  frequency.  The 
experiments  entailed  towing  model-scale  hulls 
through  ice  covers  formed  on  still  water  (i.e.  there 
was  no  water  flow  through  the  ice  tank). 

The  program  began  with  a  brief  set  of  prelimi¬ 
nary  experiments  aimed  at  gaining  an  overview  of 
the  effects  of  tow  transit  on  ice  formation.  These 
experiments  used  a  l:30-scale  tow  configuration 
comprising  conventional  barge  and  towboat  hulls. 
The  insights  gained  from  these  experiments  were 
used  in  designing  a  subsequent  set  of  experiments 
that  was  conducted  using  two  larger  (l:15-scale) 
but  radically  simplified  hull  forms.  These  experi¬ 
ments,  which  are  here  termed  the  main  experi¬ 
ments,  were  conducted  to  obtain  detailed  infor¬ 
mation  on  the  principal  effects  of  frequent  vessel 
transit  of  ice-covered  channels. 

In  addition  to  the  transit-frequency  experiments 
described  in  this  section,  two  other  laboratory 
studies  were  conducted  to  examine  significant 
features  of  tow  transit  that  were  brought  to  light 
during  the  preliminary  experiments.  One  feature. 


‘Personal  communications  with  lockmasters  on  the  Missis¬ 
sippi  and  Ohio  rivers  and  the  Illinois  Waterway. 


which  was  related  to  the  facility  of  tows  to  entrap 
and  move  broken  ice  beneath  flat-bottomed  hulls, 
was  the  subject  of  an  extensive  sideline  study.  Its 
objective  was  to  ascertain  the  equilibrium  charac¬ 
teristics  and  dimensions  of  ice  accumulations 
formed  beneath  barges  moving  through  ice-cov¬ 
ered  channels.  This  study  involved  experiments 
that  used  a  l:30-scale  barge  towed  along  a  nonre¬ 
frigerated  flume,  as  well  as  along  the  ice  tank. 

Another  study  was  conducted  to  elucidate  the 
general  characteristics  of  tow  transit  through  lay¬ 
ers  of  brash  ice.  The  interest  here  was  on  the 
manner  by  which  tows  move  and  displace  brash 
ice,  on  the  levels  of  resistance  that  they  encounter 
when  transiting  layers  of  brash  ice,  and  on  the 
merits  of  convoying  as  a  means  of  reducing  prob¬ 
lems  produced  by  frequent  transiting.  This  study 
entailed  experiments  using  the  ice  tank.  The  lay¬ 
out,  conduct  and  results  of  this  study  are  pre¬ 
sented  in  Appendix  A. 

The  ice  tank 

HHR's  21  -m-long,  5-m-wide,  1 .3-m-deep  ice  tank 
was  used  to  simulate  a  length  of  ice-covered  chan¬ 
nel.  The  layout  of  the  ice  tank  and  the  refrigerated 
laboratory  housing  it  are  shown  schematically  in 
Figure  17,  which  also  indicates  the  alignment  of 
the  track  along  which  the  model-scale  hulls  were 
towed.  The  track  width,  for  most  of  the  experi¬ 
ments,  was  about  0.7-1 .0  m.  The  model-scale  hull" 
were  towed  by  means  of  the  tank's  motorized 
towing  carriage,  which  rims  along  rails  on  the 
tank's  walls.  The  ice  tank,  its  towing  carriage  and 
the  refrigerated  laboratory  can  be  seen  in  several 
of  the  illustrations  presented  later. 

The  air  in  the  refrigerated  laboratory,  which  is 
543  m3  in  volume,  was  chilled  and  circulated 
through  a  refrigeration  system  comprising  four 
air-chillers  and  eight  dispersion  ducts  running  the 
length  of  the  ice  tank  but  located  about  2.36  m 
above  it.  Chilled  air  was  discharged  along  the 
ducts,  which  dispersed  it  approximately  uniformly 
over  the  ice  tank  such  that  ice  sheets  grew  uni¬ 
formly  to  within  10%  of  the  mean  thickness.  Air 
was  recirculated  into  the  chiller  units  after  passing 
over  the  ice  tank. 

The  ice  tank  was  filled  to  depths  of 1 .07, 0.53 and 
0.22  m  with  freshwater  doped  with  urea  at  1.0% 
concentration  by  weight.  At  a  depth  of  1 .07  m,  the 
ice  tank  can  be  considered  as  a  deep  channel,  in 
relationship  to  the  draft  of  the  model  hulls.  At 
depths  of  0.22  and  0.53  m,  the  ice  tank  represented 
a  shallow  channel  and  a  channel  of  intermediate 
depth,  respectively. 
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Model  hulls 

The  experiments  were  conducted  using  model- 
scale  hulls  simulating  several  forms  of  vessel.  The 
preliminary  experiments,  which  provided  an 
overview  of  transit-frequency  effects,  involved 
l:30-scale  hulls  of  a  river  towboat  and  two  jumbo 
barges.  These  hulls  were  made  of  wood  and  had 
the  dimensions  and  geometry  shown  in  Figures  1 8 
and  19;  both  model  and  full-scale  dimensions  are 
given .  The  hull  forms  of  the  models  are  identical  to 
those  of  the  full-scale  towboat  and  barges  used  in 
the  field  study  performed  by  Ashton  et  al.  (1973). 
The  towboat  and  barges  were  configured  such  that 
the  towboat  pushed  the  two  barges  in  parallel 
configuration  (Fig.  20).  This  configuration  is  termed 
the  1 :30-scale  tow.  It  was  moved  along  the  channel 
by  means  of  a  towing  dynamometer  mounted 
beneath  the  towing  carriage.  A  guide  pin  fixed  to 
the  tow's  bow  and  a  guide  bracket  that  extended 
from  the  carriage  were  used  to  ensure  that  the 
model  held  course  along  the  centerline  of  the 
channel. 

The  main  set  of  experiments  used  two  simpli¬ 
fied  hull  forms  that  are  nominally  1 :15-scale  analo¬ 
gies  of  tows  and  ship-form  hulls.  One  hull,  which 
was  flat-bottomed,  is  termed  the  l:15-scale  tow 
hull  as  it  is  intended  to  simulate  a  tow,  albeit  a 
shortened  and  greatly  simplified  one.  Its  dimen¬ 
sions  and  geometry  are  shown  in  Figure  21,  and  it 
is  illustrated  in  Figure  22.  The  other  hull  was 
wedge-shaped  and  is  termed  the  1: 15-scale  wedge 
hull.  It  is  intended  to  simulate  a  ship-form  hull  that 
does  not  entrap  and  shove  broken  ice.  Alterna¬ 
tively  it  can  be  viewed  as  an  idealized  tow  whose 
underside  has  developed  a  rounded  accumulation 
of  brash  ice  such  that,  when  transiting  an  ice  cover, 
it  no  longer  accumulates  and  transports  additional 
brash  ice.  The  dimensions  and  geometry  of  the 
wedge  hull  are  shown  in  Figure  23,  and  it  is  illus¬ 
trated  in  Figure  24.  The  two  hulls  are  considered  to 
be  1:15  in  scale  because  their  drafts  and  beams  are 
at  1:15  scale  of  the  full-size  towboat  and  barges 
depicted  in  Figures  18  and  19.  The  wedge  hull  has 
the  same  waterline  beam  as  the  tow  hull.  Both  the 
tow  and  the  wedge  hulls  were  fabricated  from  3- 
mm-thick,  mild-steel  plate  so  that  they  could  with¬ 
stand  unsealed  ice  loads  and  ice  abrasion  when 
transiting  the  ice  tank. 

The  tow  hull  was  fitted  with  two  0. 13-m-diame- 
ter  propellers  for  secondary  breaking  of  ice  and  for 
creating  water  currents  around  the  rear  of  the  hull. 
The  propellers  were  driven  by  a  186-W  electric 
motor,  which  is  visible  in  Figure  22a.  The  propul¬ 
sion  system  was  not  adequate  to  drive  the  hull 


through  the  ice,  though  it  could  propel  the  hull 
through  open  water.  The  tow  hull  was  connected 
to  the  towing  carriage  by  way  of  two  rigid  towing 
posts  bolted  to  the  towing  carriage  (Fig.  22a). 
Hinged  joints  between  the  hull  and  these  posts 
enabled  the  hull  to  pitch  but  not  to  move  in  any 
other  mode.  The  wedge  hull  was  rigidly  fitted  to 
the  towing  carriage  by  way  of  a  single  steel  post 
(Fig.  24).  It  could  not  pitch.  The  posts  supporting 
both  hulls  were  telescopic,  so  the  hulls  could  be 
raised  and  lowered  out  of  the  channel  and  so  that 
they  could  be  towed  at  different  depths.  Between 
transits  the  hulls  were  hoisted  from  the  channel  to 
ensure  that  they  did  not  become  frozen  into  it. 
During  each  experiment,  especially  during  peri¬ 
ods  when  the  hulls  were  suspended  in  frigid  air, 
the  hulls  were  kept  warm  with  a  1 .5-kW  warm-air 
fan  operating  inside  them.  This  was  adequate  for 
keeping  ice  from  freezing  to  the  two  hulls. 

Instrumentation 

The  air  temperature  and  the  temperature  of  the 
urea  solution  were  measured  using  platinum  re¬ 
sistance  thermometers  (RTDs)  arrayed  over  and  in 
the  channel.  Each  probe  was  of  0.01  °C  precision.  A 
post  bearing  four  RTD  probes  was  located  half¬ 
way  along  the  channel  and  close  to  the  tracks  cut 
through  it,  as  indicated  in  Figure  17.  Two  probes 
for  measuring  the  temperature  of  the  urea  solution 
were  positioned  on  the  post  at  0.30  and  0.10  m 
(deep  channel  only)  below  the  water  level.  The 
remaining  two  probes  were  positioned  0.05  and 
0.20  m  above  the  water  level.  Three  additional 
RTD  probes  were  placed  1.86  m  above  the  water 
level  in  the  ice  tank.  Temperature  data  from  the 
temperature  probes  were  automatically  logged  at 
quarter-hour  intervals  using  IIHR's  HP1000  com¬ 
puter. 

The  thicknesses  of  brash-ice  accumulations,  as 
a  layer  covering  the  track  or  as  ridges  flanking  it, 
were  measured  over  a  2-m-long  area  halfway  along 
the  channel.  This  area,  designated  in  Figure  17  as 
the  monitoring  area,  is  adjacent  to  the  temperature 
post  in  the  tank.  The  thicknesses  of  the  ice  sheet  q 
and  the  ice  regrowth  t|n  were  measured  with  a 
vernier  caliper  of  0.1 -mm  precision.  Each  meas¬ 
urement  entailed  cutting  and  removing  a  sample 
from  the  ice  sheet  or  track. 

At  each  time  of  measurement,  the  thicknesses  of 
the  layer  qb  and  ridge  T|r  were  averaged  from  at 
least  four  measurements  made  at  different  sec¬ 
tions  in  the  monitoring  area.  The  measurements 
were  made  using  a  specially  developed,  hook¬ 
shaped  pointer  fitted  to  a  vernier  scale  that  was 
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2.76 

(0.09) 


Figure  19.  Dimensions  and  geometry  of  the  l:3()~scale  barges.  All  dimensions  are  in  meters. 
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a.  Transiting  the  ice  tank. 


b.  Stowed  beneath  the  towing  carriage. 


Figure  20.  l:30-scale  tow  model. 
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Figure  21.  Dimensions  and  geometry  of  the  l:15-scale  tow  hull.  All  dimensions  are  in  centimeters. 
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BOW  VIEW  STERN  VIEW 

Figure  23.  Dimensions  and  geometry  of  the  l:15-scale  wedge  hull.  All  dimensions  are  in  centimeters. 
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b.  Stowed  beneath  the  toiving  carriage. 


Figure  24.  l:15-scale  wedge  hull. 


mounted  on  a  rigid  cross  beam.  The  precision  of 
the  vernier  scale  was  0.1  mm.  For  the  deep-channel 
experiments,  the  beam  was  placed  across  two 
working  platforms  hung  in  front  of  the  towing 
carriage.  For  the  shallow-channel  experiments, 
the  beam  was  placed  across  concrete-block  sup¬ 
ports  seated  on  the  tank  bottom. 

During  the  preliminary  experiments  with  the 
l:30-scale  tow,  and  subsequently  with  the  1:15- 
scale  tow  hull,  brash  ice  collected  in  relatively 
thick  accumulations  at  the  ends  of  the  channel. 
Such  accumulations,  here  termed  plugs,  were 
found  to  be  an  important  consequence  of  frequent 
transits  as  they  eventually  blocked,  or  plugged, 
tracks.  The  maximum  thickness  of  plug  was 
measured  using  the  hooked  point  gauge. 

The  porosity  p  of  brash-ice  accumulations  was 
measured  by  removing  a  portion  of  each  accumu¬ 
lation  with  a  perforated  0.20-m-diameter  canister 
of  known  volume  and  weight.  The  accuracy  with 
which  p  could  be  determined  varied  with  the  size 
of  individual  pieces  of  brash  ice  relative  to  the0.20- 
m-diameter  opening  of  the  canister;  finer  brash 
resulted  in  more  accurate  estimates. 

The  areal  concentration  of  the  open-water 
surface  (3  over  tracks  formed  by  the  hulls  was 
determined  after  each  transit  by  planimetry  of 
photographs  taken  above  the  monitoring  area. 

The  size  distribution  of  brash  ice  in  the  track 
was  determined  at  the  end  of  each  experiment.  A 
0.30-m-diameter  bucket  was  used  to  obtain  samples 
of  brash  ice,  which  were  then  visually  sorted  into 
several  size  classes  and  photographed.  Appendix 
B  contains  photographs  of  brash-ice  samples.  The 
orthogonal  dimensions  (dx,  d  and  d  )  of  ice  pieces 
were  measured  so  as  to  obtain  an  equivalent 
dimension: 

de  =  3Jdxdydz  (24) 

for  each  class.  Finally,  each  class  was  weighed 
such  that  a  histogram  and  a  cumulative-frequency 
curve  of  Rvalues  could  be  prepared.  These  were 
assessed  for  each  experiment,  in  order  to  charac¬ 
terize  the  resulting  brash  ice  in  terms  of  a  median 
dimension: 

Je  =  iflde  (25) 

N  n  =  1 


in  which  d£&4  and  dfil6  are  values  of  dg  for  which 
84%  and  16%,  respectively,  of  the  brash-ice  pieces 
are  finer  in  size. 

During  and  after  several  experiments  (includ¬ 
ing  all  the  preliminary  ones),  profiles  of  brash  ice 
accumulated  along  the  towing  track  were  meas¬ 
ured  and  recorded  using  an  ultrasonic  depth- 
sounder.  The  sonar  head  of  the  sounder  was  at¬ 
tached  to  an  underwater  towing  carriage  that  ran 
along  the  video-camera  trench  indicated  in  Figure 
17.  The  underwater  carriage  was  connected  to  the 
towing  carriage  via  a  cable  and  pulley  tow  mecha¬ 
nism.  Its  use  was  eventually  abandoned.  Although 
the  sounder  enabled  profiles  to  be  determined 
with  acceptable  precision  (to  within  2  mm),  it  had 
to  be  towed  slowly  (14  mm/s).  Furthermore,  as 
transits  by  the  tow  hulls  produced  brash-ice  accu¬ 
mulations  that  were  of  irregular  thickness,  the  in¬ 
formation  yielded  by  the  sounder  proved  to  be  of 
little  use. 

Model  ice 

As  the  primary  focus  of  this  study  is  on  ice 
growth,  or  ice  formation,  the  experiments  were 
conducted  using  an  aqueous  model  ice  that  can 
grow  thermally  in  a  similar  manner  as  ice  on 
rivers,  lakes  and  seas.  However,  a  major  impedi¬ 
ment  to  any  laboratory  or  scale-model  experiment 
involving  ice  growth  as  well  as  ship  (or  structure) 
interaction  with  ice  is  that  both  ice  growth  proc¬ 
esses  and  structural  interactions  cannot  be  mod¬ 
eled  simultaneously.  In  fact,  there  is  no  model  ice 
that,  with  complete  fidelity,  can  be  used  for  model- 
scale  simulations  of  either  the  growth  or  strength 
properties  of  ice. 

To  overcome  this  impediment,  these  experi¬ 
ments  were  conducted  without  strict  adherence  to 
similitude  criteria  for  either  ice  growth  or  ice- 
vessel  interaction.  Instead,  urea  ice,  which  is  a 
common  model  ice  that  can  be  grown  in  ice  tanks, 
was  used  to  replicate  the  major  features  of  ice 
grown  on  rivers,  lakes  or  seas.  Urea  acts  as  a  brine 
included  in  ice  sheets.  Warming,  or  tempering,  is 
a  means  of  weakening  ice  sheets  so  that  they  can  be 
used  in  scale  models  of  vessel  and  structure  inter¬ 
actions  with  ice  sheets.  A  limitation  for  the  present 
study,  however,  is  that  the  ice  could  not  be  tem¬ 
pered  prior  to  each  transit  of  a  model  hull. 


and  a  geometric  standard  deviation: 
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Scaling  relationships 

In  accordance  with  the  ice-growth  and  ice-ves¬ 
sel-interaction  facets  of  the  experiments,  two  time 
(26)  scales  were  required.  One  time  scale  is  based  on 
conventional  Froude-  and  Cauchy-number  scal- 
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ing  of  ice  interactions  with  vessels  or  structures. 
The  other  scale  is  based  on  ice  growth  rates  in  the 
ice  tank  relative  to  typical  rates  of  ice  growth  on 
rivers. 

In  accordance  with  the  Froude-  and  Cauchy- 
number  scaling  criteria,  the  time  scale  associated 
with  hull  and  ice-piece  motion  during  vessel  tran¬ 
sit  is 


(27) 


in  which  X^is  the  geometric  scale  of  the  hull. 

The  time  scale  associated  with  ice  growth  in  the 
ice  tank  can  be  evaluated  using  eq  4.  Ice  growth 
rate  can  be  expressed  as 


^H  =  _l - 41 —  .  (28) 

dt  PiLi  R  ,  1 

^  ia 


Integrating  eq  28,  and  using  the  initial  condition  t| 
=  0  at  t  =  0,  as  well  as  assuming  that  T  is  constant, 
leads  to 
-.2 
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Normalizing  q  with  barge  draft  D  or  another 
appropriate  vessel-related  length,  and  rearrang¬ 
ing  eq  29,  yields 


(q\ 2  (  q2fcj  _2fcj  AT  t 
W  D2hia  PjhjD2 


(30) 


From  the  first  term,  the  scale  of  ice-sheet  thickness, 
=  Xp  =  X(.  From  the  first  two  terms,  X^  =  Xhia  Xp. 
From  all  three  terms,  with  Xp.  =  X^  =  1,  the  time 
scale  associated  with  ice  growth  is 


Xt  = 


(31) 


The  selection  of  an  exact  value  for  the  time  scale 
associated  with  ice  growth  clearly  entails  subjec¬ 
tive  assessment  of  the  representative  heat-flux 
conditions  for  channels  in  nature.  If,  to  consider 
the  simplest  case,  X^  and  XATare  each  taken  to  be 
unity,  eq  31  implies  tkat  the  time  scale  associated 
with  ice  growth  is  equivalent  to  the  geometric 
scale.  In  other  words,  for  geometric  scales  of  X(=  30 
and  15,  one  hour  of  ice  growth  in  the  laboratory 
channel  is  equivalent  to  30  and  15  hours,  respec¬ 
tively,  at  full  scale. 

In  nature,  h.  and  AT  are  unsteady,  at  least  for 


periods  in  excess  of  one  day,  and  they  are  difficult 
to  assess.  Shen  (1985)  estimated  an  h.  value  of 
about  20  W/m2oC  for  the  St.  Lawrence  River. 
Ashton  (1986)  recommended  that  heat  fluxes  from 
rivers  under  frigid  conditions  be  assessed  as  200 
W/m2.  Values  of  heat  flux  from  the  laboratory 
channel  were  estimated  by  eq  5  to  be  about  70  W/ 
m2  when  air  in  the  refrigerated  laboratory  was 
maintained  at  -10°C.  Therefore,  if  the  product 
X^.  XAJis  taken  as  200/70,  X(  would  be  about  10 
and  5  for  X(  =  30  or  15,  respectively. 

Program  of  experiments 

A  two-part  program  of  experiments  was  con¬ 
ducted  .  The  first  part  comprised  three  preliminary 
experiments  conducted  with  the  l:30-scale  tow 
hulls.  The  second  part,  or  main  experiments, 
comprised  thirteen  experiments  conducted  with 
the  l:15-scale  tow  and  wedge  hulls.  Details  of  the 
experiments  are  summarized  in  Table  1. 

All  but  two  experiments  involved  transits  at  a 
single  frequency,  such  as  one  transit  every  two 
hours.  The  two  experiments  simulate  a  general 
schedule  of  transiting  as  depicted  in  Figure  16.  For 
all  experiments  the  air  temperature  above  the  ice 
tank  was  maintained  at  -10  ±  1.5°C,  except  for 
short  designated  periods  (as  explained  later).  The 
±  1 .5°C  variation  in  air  temperature  was  caused  by 
automatic  defrosting  of  the  chiller  units,  which 
took  place  at  two-hour  intervals. 

The  preliminary  experiments  were  intended  to 
gain  an  overview  of  the  effects  of  tow  transit  on  ice 
formation  in  a  navigable  channel.  The  results  of 
these  experiments  guided  the  design  and  program 
of  subsequent  experiments.  The  first  experiment 
involved  11  transits  of  the  l:30-scale  tow  at  a 
transit  frequency  /  of  one  every  four  hours.  The 
second  experiment  was  conducted  with  the  tow 
moving  alternately  along  two  overlapping  tracks, 
with  /  =  1/hour.  The  third  experiment  involved 
only  the  towboat. 

The  objective  of  the  main  set  of  experiments 
was  to  obtain  detailed  information  on  the  effects  of 
vessel  transit  on  ice-cover  formation  for  a  range  of 
transit  frequencies  and  channel  conditions.  The 
principal  variables  were  transit  frequency  hull 
shape  and  channel  depth  y  .  Values  of  /  ranged 
from  1/(0 .5  hour)  to  1/(8  kours)  at  1:15  model 
scale.  Each  experiment,  except  one,  lasted  about 
three  to  four  days  and  ended  when  the  ice  sheet 
had  attained  a  thickness  of  about  60-70  mm 
(0.90-1.05  m  at  full  scale).  A  brief  experiment 
lasting  4.2  hours  was  performed  at  the  high  tran¬ 
siting  frequency  of  1/(0.167  hours).  Its  purpose 
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Table  1.  List  of  experiments  on  ice  regrowth. 


a 


Experiment 

number 

p; 

(hour) 

n 

P 

S 

(hour) 

P» 

(hour) 

Pr 

(hour) 

In 

T 

(°C) 

y0 

(m) 

I mm 1 

( °C-day)0Sl 

Note 

Part  1 .  Preliminary  Experiments  with  1:30-Scale  Hulls 

1  4.0  1  4.0  0 

2.04 

11 

-9.27 

1.07 

10.0 

Tow 

2 

1.0 

i 

1.0 

0 

1.07 

12 

-10.88 

1.07 

9.8 

Alternate 

3 

2.0 

i 

2.0 

0 

1.44 

13 

-9.17 

1.07 

9.7 

tracking 
of  tow 

Towboat  only 

Part  2.  Main  Experiments  with  1 :15-Scale  Hulls:  Deep  Channel 

1  2.0  1  2.0  0  4.34 

29 

-10.05 

1.07 

11.8 

Tow 

2 

4.0 

i 

4.0 

0 

2.98 

15 

-10.38 

1.07 

13.0 

" 

3 

0.5 

4 

15 

6 

2.33 

32 

-10.56 

1.07 

13.0 

" 

4 

8.0 

1 

8.0 

0 

3.00 

9 

-10.82 

1.07 

125 

" 

8 

4.0 

1 

4.0 

0 

3.67 

18 

-8.79 

1.07 

13.0 

Wedge 

9 

2.0 

1 

2.0 

0 

4.00 

30 

-11.74 

1.07 

14.8 

10 

0.5 

4 

1.5 

6 

4.00 

36 

-12.78 

1.07 

13.0 

" 

11 

8.0 

1 

8.0 

0 

3.10 

9 

-10.81 

1.07 

13.2 

" 

12 

0.5 

1 

0.5 

0 

3.10 

98 

-10.88 

1.07 

13.0 

// 

13 

0.17 

1 

0.17 

0 

0.18 

25 

-10.00 

1.07 

13.0 

" 

Part  3.  Main  Experiments  with  l:15-Scale  Hulls:  Shallow  Channel 

5  4.0  1  4.0  0  3.00 

15 

-9.41 

0.22 

13.6 

Tow 

6 

4.0 

1 

4.0 

0 

2.94 

14 

-9.53 

053 

12.9 

u 

7 

2.0 

1 

2.0 

0 

2.75 

24 

-11.02 

0.22 

12.7 

u 

*P(  =  period  between  transits; 
n  =  number  of  transits  during  navigation  period; 

Ps  =  period  associated  with  a  sequence  of  transits; 

Po  =  period  between  sequences  of  transits; 

PT  =  P.+Z(Ps+Po)  =  total  period  of  navigation,  or  ice  growth; 
Pj  =  initial  period  of  ice  growth  here  set  as  P(; 

In  =  number  of  transits  in  the  navigation  period; 

Ta  =  air  temperature 
yo  =  water  depth 
a  =  Stefan  coefficient. 


was  to  determine  how  p  varied  with  n  at  high 
values  of /(. 

The  effects  of  tow  transits  on  ice  formation  were 
investigated  for  three  channel  depths;  1.07,  0.53 
and  0.22  m.  The  1.07-m-deep  channel  is  consid¬ 
ered  to  be  deep  because  the  channel  bottom  was 
distant  from  the  underside  of  the  tow  hull.  The 
0.22-m-deep  channel  was  used  to  simulate  a  shal¬ 
low  channel,  with  a  minimum  navigation  depth  of 
3.05  m  (10  ft).  The  water  depth  was  not  varied  for 
experiments  with  the  wedge  hull  because  no  ice 
passed  beneath  it,  so  the  water  depth  would  not 
significantly  affect  the  movement  of  broken  ice 
around  its  inclined  sides. 

Experimental  procedure 

At  the  beginning  of  each  experiment,  the  air 
above  the  ice  tank  was  chilled  to  an  average  tem¬ 


perature  of  -10°C.  The  urea  solution  in  the  tank 
simultaneously  cooled  to  -0.1  ±  0.02°C.  The  sur¬ 
face  of  the  tank  was  then  wet-seeded  by  means  of 
a  fog  of  fine  ice  crystals.  The  nascent  ice  sheet  was 
allowed  to  grow  for  an  initial  period  P.,  the  same 
as  the  interval  between  vessel  transits  /^.  Its  thick¬ 
ness  was  measured  immediately  prior  to  each 
transit.  The  model  hull  for  each  experiment  was 
launched  from,  and  coupled  to,  the  towing  car¬ 
riage  at  one  end  of  the  ice  tank.  At  this  stage  the 
propellers  of  the  1 : 1 5-scale  tow  were  set  in  motion. 

Each  hull  was  then  accelerated  and  towed  at  a 
constant  speed  of  0.15  m/s  (about  0.58  m/s,  or 
about  2  ft/s,  at  full  scale)  along  the  simulated 
navigation  channel.  While  each  hull  was  being 
towed  through  the  channel,  the  track  behind  it  was 
wet-seeded  so  as  to  initiate  regrowth  of  ice  over 
the  track.  For  several  of  the  experiments,  tow 
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transits  and  ice  formation  were  videotaped. 

On  completion  of  a  transit,  each  hull  was 
winched  from  the  channel  and  held  beneath  the 
towing  carriage  until  the  subsequent  transit.  The 
experiments  were  conducted  continuously  on  a 
"round-the-clock"  basis  for  about  three  to  four 
days  until  the  ice  cover  over  the  tank  had  grown  to 
an  equivalent  prototype  thickness  of  about  1  m 
(thicknesses  greater  than  1  m  are  extreme  for  most 
navigable  channels,  particularly  those  of  the  up¬ 
per  Mississippi  River  and  its  tributaries).  Transit¬ 
ing  was  interrupted  only  for  the  period  from 
approximately  0000  hours  to  0500  hours,  during 
which  the  air  temperature  was  adjusted  so  as  to 
give  the  same  incremental  value  of  as  for  the 
regular  interval  between  transits. 

Throughout  each  experiment  the  following 
quantities  were  measured: 

•  Thickness  of  the  ice  sheet  t|; 

•  Thickness  of  the  ice  regrowth  over  open  water 

•  Average  thickness  of  the  ice  accumulation  in 
the  navigation  channel  T|b; 

•  Average  thickness  of  the  ice  accumulation  as 
a  ridge 

•  Thickness  of  plugs  of  brash  ice  n  (for  transits 
of  the  tow-hull); 

•  Porosity  p  of  brash  ice  accumulations; 

•  Areal  concentration  of  open  water  over  the 
track  P; 

•  Size  distribution  of  brash  ice  (on  completion 
of  an  experiment). 

The  first  two  measurements  were  made  imme¬ 
diately  prior  to  each  transit,  whereas  the  remain¬ 
der  were  made  after  each  transit.  For  a  few  experi¬ 
ments  the  variation  of  the  average  plan  dimension 
of  brash  ice  on  the  surface  of  the  track  was  meas¬ 
ured  after  each  transit.  Also  measured  occasion¬ 
ally  were  the  ridge  geometry  and  the  time  re¬ 
quired  for  water  and  brash  ice  in  the  track  to  calm 
sufficiently  so  that  ice  regrowth  could  commence. 

Experimental  errors 

One  significant  source  of  error  arises  because 
the  experiments  were  performed  without  strict 
adherence  to  similarity  criteria  for  modeling  ves¬ 
sel  transit  through  ice.  These  criteria  had  to  be 
relaxed,  as  explained  earlier,  but  the  penalties  of 
doing  so  have  to  be  kept  in  mind.  The  following 
penalties,  or  scale  effects,  are  perhaps  the  most 
important  ones: 

•  Thermodynamic  terms  that  are  significant 
when  describing  thin  ice  (e.g.  the  thermal 
boundary  layer)  are  not  taken  into  account. 


•  Brash-ice  fragments  formed  in  the  ice  tank 
may,  for  several  reasons,  be  larger  than  the 
equivalent  brash  ice  formed  at  full  scale. 
Incomplete  scaling  of  ice  strength  and  defor¬ 
mation  characteristics  may  have  resulted  in 
the  formation  of  initially  larger  sizes  of  bro¬ 
ken  ice  and  brash  ice  that  may  have  been  less 
readily  broken  during  impact  with  the  hull. 
Lack  of  direct  control  over  the  thermody¬ 
namics  of  ice  growth  at  the  crystal  scale  may 
have  caused  contacting  ice  fragments  to  fuse 
and  form  larger  fragments. 

•  At  the  highest  transit  frequency  [/t=  1  /(0.167 
hour)],  the  resulting  brash  ice  behaved  cohe¬ 
sively,  similarly  to  frazil  ice.  These  ice  frag¬ 
ments  tended  to  bunch  as  floes  (akin  to  frazil 
ice)  rather  than  to  remain  discrete  ice  frag¬ 
ments.  Therefore,  this  ice  may  not  replicate 
the  behavior  of  full-scale  ice,  which  would  be 
less  inclined  to  form  floes  (an  analogy  here  is 
the  dissimilar  behavior  of  clay  and  sand-size 
sediments).  However,  intensely  transited 
channels  may  produce  mush  ice  that  can 
behave  cohesively. 

•  Incomplete  scaling  of  ice  forces  and  hull 
motions  may  have  affected  brash  ice  dis¬ 
placement  and  distribution.  (The  ice  resis¬ 
tance  exerted  against  the  rigid  model  hulls 
was  strong  enough  for  thick  ice  to  scrape 
paint  off  them). 

Despite  concerns  for  the  significant  scale  effects 
inherent  in  scale  modeling  of  ice  growth  as  well  as 
ice  interaction  with  vessels,  the  experiments  do 
illuminate  most  of  the  more  important  features  of 
ice  formation  in  frequently  transited  ice-covered 
channels.  In  a  simpler  sense,  one  that  sidesteps  the 
nettlesome  issue  of  scaling,  the  main  experiments 
can  also  be  viewed  as  a  physical  analogy  of  tow 
transit.  In  accordance  with  this  view,  the  steel 
hulls  are  essentially  greatly  simplified  tows,  or 
tow-like  hulls,  that  transit  an  element  of  ice-cov¬ 
ered  navigation  channel. 

Errors  attributable  to  instrument  precision  are 
much  smaller  than  errors  resulting  from  the  gen¬ 
eral  design  and  conduct  of  the  experiments.  For 
example,  the  mechanical  operation  of  the  refrig¬ 
eration  plant  caused  the  air  temperature  above  the 
ice  tank  to  fluctuate  within  a  tolerance  of  1.5°C, 
which  exceeded  the  0.01  °C  precision  of  the  RTD 
temperature  probes.  However,  a  major  source  of 
error  resulted  from  the  nonuniform  distribution  of 
brash  ice  accumulated  in,  or  next  to,  the  tracks 
formed  by  the  hulls.  This  error  was  because,  for 
some  of  the  experiments,  the  size  of  individual 
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pieces  of  brash  ice  in  the  accumulations  comprised 
up  to  50%  of  accumulation  thickness.  Therefore, 
the  accuracy  of  accumulation  thicknesses  reported 
herein  is  typically  of  the  order  of  the  mean  dimen¬ 
sion  of  brash-ice  fragments. 

ICE-TANK  OBSERVATIONS  AND  DATA 

Frequent  transits  of  the  ice  tank  by  both  the  1:30- 
and  l:15-scale  tow  hulls  revealed  a  pattern  of  ice 
formation  that  is  markedly  different  from  that 
portrayed  in  Figure  14  and  that  assumed  in  the 
development  of  prior  predictive  models.  Ice  for¬ 
mation  in  channels  transited  by  the  tow  hulls  was 
dominated  by  the  propensity  of  these  hulls  to 
entrap,  shove  and  drag  ice  along  the  tracks  that 
they  opened.  The  result  was  the  formation  of 
localized  thick  accumulations  of  brash  ice,  here 
termed  ice  plugs.  Frequent  transits  of  the  ice  tank 
by  the  wedge  hull,  which  did  not  move  brash  ice 
along  its  track,  resulted  in  essentially  the  same 
sequence  of  ice  formation  portrayed  in  Figure  14. 
The  experiments  demonstrated  that,  among  other 
factors,  hull  shape  significantly  affects  ice  forma¬ 
tion  in  frequently  transited  ice-covered  channels. 
The  experiments  also  demonstrated  that  accumu¬ 
lations  of  brash  ice  become  especially  troublesome 
when  shallow  ice-covered  channels  are  transited. 

Overview  of  ice  formation  in  channels  transited 
by  tows 

When  the  tow  hulls  first  transited  the  channel 
covered  by  a  thin,  newly  formed  ice  sheet  (Fig. 
25a),  most  of  the  resulting  brash  ice  was  displaced 
from  the  track  opened  by  the  tow.  Some  of  it  was 
cast  beneath  the  edges  of  the  adjoining  ice  sheet, 
forming  inverted  ridges  that  flanked  both  sides  of 
the  track.  Much  of  it  became  entrapped  and 
dragged  along  by  the  tow  until  it  was  deposited  at 
some  location  along  the  tow's  track  (Fig.  25b). 
Only  a  relatively  small  proportion  of  the  initial  ice 
cover  remained  as  brash  ice  strewn  over  the  track, 
some  of  it  having  been  drawn  back  into  the  track 
by  propeller  wash.  At  the  end  of  the  track  (in  actual 
channels,  the  equivalent  locations  could  be  bars, 
submerged  rock  outcrops,  bends  or  channel  con¬ 
strictions  such  as  are  evident  in  Figure  3),  broken 
ice  was  sloughed  from  beneath  the  tow  and  depos¬ 
ited  as  a  plug  lodged  in  the  track. 

Further  transits,  following  periods  of  ice  re¬ 
growth  over  the  track,  displaced  additional  ice 
sideways  to  the  ridges  or  caused  it  to  be  dragged 
and  shoved  to  the  plugs  formed  in  the  track.  The 


areal  concentration  of  brash  ice  spread  over  the 
track  usually  remained  small,  and  the  track  itself 
remained  straight-edged,  while  the  channel  was 
covered  by  a  thin  ice  sheet.  The  experiments  sug¬ 
gest  that,  when  the  ice  sheet  is  less  than  about  0.30 
m  thick  at  full  scale,  its  edges  become  reamed  by 
the  sides  of  the  tow. 

The  ice  sheet  continued  to  thicken  and 
strengthen  and  encroached  laterally  inwards  (as  a 
form  of  border  ice  growth)  along  the  track.  Near 
the  track,  thickening  of  the  ice  sheet  was  aug¬ 
mented  by  freezing  of  water  washed  over  it  during 
transits.  The  edges  of  the  track  thickened  and 
strengthened  such  that,  when  impacted  by  the 
tow's  bow,  comparatively  large  cuspate  pieces  of 
ice  broke  flexurally  from  the  ice  sheet.  These  pieces 
drifted  into  the  track  (Fig.  25c).  Their  compara¬ 
tively  large  size  caused  them  to  pass  less  readily 
beneath  the  tows.  Instead,  they  were  either  dis¬ 
placed  beneath  the  edges  of  the  ice  sheet,  becom¬ 
ing  part  of  the  ridge,  or  were  shoved  to  plugs 
formed  in  the  track.  Although  the  amount  of  ice  in 
the  track  initially  increased  with  the  entry  of  bro¬ 
ken  sheet  ice,  it  was  reduced  by  subsequent  tran¬ 
sits,  which  cleared  it  from  the  track  (compare  Fig. 
25c  and  d).  Regrowth  of  ice,  together  with  the 
freezing  of  wave-generated  wash  from  subsequent 
transits,  smoothed  serrated  edges  and  sealed  cracks 
in  the  ice  sheet  (Fig.  25e)  at  some  locations.  At 
others,  fresh  breaks  occurred. 

Fed  with  brash  and  broken  ice  during  succes¬ 
sive  tow  transits,  the  plugs  thickened  and  ex¬ 
tended  along  the  track  (Fig.  25f).  In  the  ice  tank, 
plugs  grew  as  ice  was  dislodged  and  then  depos¬ 
ited  when  the  tow  hulls  overrode  them.  In  river 
channels  with  substantial  water  current,  brash  ice 
may  be  conveyed  along  a  track  and  collected  at 
plug  sites.  Eventually,  after  numerous  transits, 
and  depending  on  plug  spacing,  the  track  becomes 
filled  with  broken  ice.  In  a  sense  the  ice  formation 
processes  illustrated  in  Figure  25f  represent  a 
pattern  of  lateral  accumulation  of  brash  and  bro¬ 
ken  ice  along  a  track,  as  opposed  to  the  pattern  of 
vertical  accumulation,  or  thickening,  occurring 
simultaneously  along  a  track  (Fig.  14). 

The  experiments  revealed  thattransport  of  brash 
ice  and  plug  formation  are  features  that  may 
dominate  tow  transits  of  ice-covered  channels. 
However,  a  concern  arises  as  to  whether  this  find¬ 
ing  is  merely  a  quirk  of  the  experiments  or  whether 
it  also  holds  at  full  scale.  Brash-ice  accumulations 
beneath  tows  that  have  traveled  long  distances 
may  have  attained  an  equilibrium  such  that  they 
no  longer  entrap  and  transport  significant  quanti- 
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Figure  25.  Ice  formation  in  channels  transited  by  tows  that  move  ice  along  tracks. 
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ties  of  brash  ice.  The  ice-tank  experiments,  how¬ 
ever,  were  conducted  with  an  initially  ice-free  tow 
transiting  a  relatively  short  channel;  it  was  only  1 0 
tow-lengths  long  (nominally  600  m  of  channel)  at 
1 :30  scale,  or  20  tow-lengths  long  (nominally  1 200 
m  of  channel),  if  the  1 : 15-scale  tow  hull  is  assumed 
to  simulate  a  60-m-long  tow. 

To  address  this  concern,  experiments  were 
conducted  using  the  wedge  hull.  This  hull  was 
intended  to  simulate  a  hull  form  that  does  not 
entrap  or  shove  ice.  It  could  be  viewed  as  being 
either  a  tow  that  is  enveloped  by  an  equilibrium 
accumulation  of  broken  ice  or  more  simply  as  a 
ship-form  hull.  Transits  of  the  laboratory  channel 
by  the  wedge  hull  resulted  in  practically  the  same 
sequence  of  ice  formation  illustrated  in  Figure  14, 
and  they  could  be  simulated  by  the  numerical 
model  presented  earlier.  In  accordance  with  this 
sequence,  most  of  the  ice  broken  during  transits  by 
the  wedge  hull  was  passed  beneath  the  hull  and 
remained  in  the  track,  or  it  was  displaced  to  ridges 
formed  beneath  the  adjoining  edges  of  the  ice 
sheet.  The  wedge  hull  did  not  shove  significant 
quantities  of  brash  ice  along  the  channel. 

Though  frequent  transits  of  the  ice  tank  by  the 
tow  and  wedge  hulls  generally  resulted  in  in¬ 
creased  volume  of  ice  growth,  the  increase  was 
never  strikingly  large.  Atmost,  transitingresulted 
in  a  doubling  of  the  volume  of  solid  ice.  For  both 
hulls,  ridges  formed  along  the  tracks  were  about 
three  times  the  ice-cover  thickness,  whereas  the 
l:30-scale  tow  model  formed  ridges  that  were 
about  five  times  the  ice-cover  thickness.  The  tracks 
formed  by  the  tow  hulls  usually  contained  negli¬ 
gible  amounts  of  brash  ice,  most  of  it  having  been 
pushed  to  the  plugs  that  developed  at  the  ends  of 
the  tracks.  The  plug  thickness  can  be  on  the  order 
of  ten  (or  more)  times  the  ice  sheet  thickness. 

Plug  and  ridge  formation  pose  problems  in  ice- 
covered  channels  because  they  are  resistant  to 
transit.  They  may  also  be  troublesome  when  ice 
covers  break  up  if  they  precipitate  ice  jamming. 
The  problems  associated  with  ridge  and  plug  for¬ 
mation  are  exacerbated  in  shallow  channels  that 
provide  little  clearance  beneath  vessels.  In  such 
channels,  ridges  and  plugs  develop  more  rapidly 
than  in  deep  channels,  and  they  may  extend  over 
the  full  channel  depth.  The  sketches  in  Figure  26 
illustrate  how  the  lack  of  clearance  beneath  a  tow 
causes  brash  ice  to  be  displaced  sideways  to  ridges 
that  flank  the  track,  or  to  be  accumulated  as  a 
jumbled  ice  mass  shoved  ahead  of  the  tow.  Under 
these  conditions,  a  tow  may  become  grounded  on 
brash  ice  lodged  against  the  channel  bottom. 


Obviously  the  most  difficult  channels  for  tows  to 
transit  are  those  that  are  maintained  to  the  COE 
minimum  navigation  depth  of  3.05  m  (10  ft). 

Preliminary  experiments 

The  photographs  in  Figure  27  show  a  track 
opened  by  the  1 :30-scale  tow  transiting  the  ice  tank 
at  a  frequency  of /t=  1/(4  hours).  The  results  show 
how  little  ice  remained  in  the  track  opened  by  the 
model  tow  and  how  rapidly  a  brash-ice  plug  formed 
at  the  end  of  the  channel. 

After  the  first  two  transits,  a  negligible  amount 
of  brash  ice  remained  in  the  track  (Fig.  27a).  Most 
of  it  had  been  entrapped  beneath  the  tow,  dragged 
along  by  it,  and  deposited  as  a  plug  at  the  end  of 
the  channel  (the  plug  is  evident  in  the  foreground 
of  Figure  27a)  or  cast  sideways  beneath  the  edges 
of  the  ice  sheet  bordering  the  track.  There,  inverted 
linear  piles  of  brash  ice  accumulated  as  ridges.  The 
tows  transported  brash  ice  by  entrapping  it  along 
the  tow's  flat-bottomed  underside,  beneath  its 
relatively  bluff  bow,  and  in  small  wakes  immedi¬ 
ately  behind  the  barges  comprising  the  tow.  The 
amount  of  brash  ice  remaining  in  the  channel  after 
four  transits  was  also  negligible  (Fig.  27b).  Even 
after  ten  transits,  little  brash  ice  had  accumulated 
in  the  channel,  although  by  then  the  plug  extended 
back  along  the  channel  for  about  one  third  of  its 
length  (Fig.  27c).  The  experiment  was  halted  after 
1 1  transits,  as  negligible  amounts  of  brash  ice  had 
accumulated  in  the  track. 

The  thicknesses  of  the  ice  sheet  q  and  the  accu¬ 
mulation  as  ridges  qr  and  the  thickest  plug  q  are 
presented  in  Figure  28  vs  cumulative  cooling  rime 
ZSj.  At  the  monitoring  area,  qb/q  =  0.  Although 
the  ridges  were  by  no  means  uniformly  thick 
(being  instead  jagged),  the  average  ridge  thick¬ 
ness  normalized  with  ice  sheet  thickness  qr/q  was 
about  1.3.  As  plug  formation  was  unforeseen,  its 
thickness  was  not  measured  until  the  end  of  the 
experiment,  when  it  was  found  that  q  /q  was 
about  8. 

The  second  experiment  conducted  with  the  1 : 30- 
scale  tow  resulted  in  essentially  the  same  channel 
condition  as  illustrated  in  Figure  27.  For  this 
experiment,  /  =  1/hour,  but  the  tow  model  was 
transited  alternately  along  two  overlapping  tracks 
as  indicated  in  Figure  29  (the  overlap  being  one 
barge  width).  The  tracks  overlapped  to  replicate 
the  more  realistic  transit  situation  whereby  tows 
move  along  several  overlapping  or  adjoining  tracks 
rather  than  along  a  single  track.  However,  this 
experiment  had  essentially  the  same  outcome  as 
the  first  one.  As  indicated  in  Figure  30,  negligible 
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Figure  28.  Ice  thickness  vs  cumulative  cooling  time;  1 :30-scale  tow;  ff  =  1/(4 
hours). 


Figure  29.  Ice  formation  in  the  ice  tank  transited  by  the  1:30- 
scale  tow  along  overlapping  tracks;  f(  =  1/ hour. 
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Figure  30.  Ice  thickness  vs  cumulative  cooling  time;  l:30-scale  tow;  f,  =  1/  hour. 


brash  ice  accumulated  in  the  track,  qb/q  =  0,  T^/q 
=  3,  and  qp/ri  a  10. 

The  third  experiment  was  conducted  using  only 
the  1 :30-scale  towboat  (its  dimensions  are  given  in 
Figure  18),  which  was  transited  through  the  chan¬ 
nel  at/{=  1/(2  hours).  The  intent  of  this  experiment 
was  to  indicate  the  effects  of  tow  length  on  ice 
formation.  The  relative  shortness  and  the  rounded 
form  of  the  towboat  hull  shoved  and  dragged 
considerably  less  brash  ice  than  did  the  tow  hull. 
However,  the  towboat  left  more  brash  ice  in  its 
track  and  produced  a  smaller  plug,  as  can  be  seen 
by  comparing  Figures  31  and  27,  or  by  comparing 
the  data  in  Figures  32  and  28.  This  finding  is  likely 
identical  to  that  observed  at  full  scale  by  Ashton  et 
al.  (1973). 

Although  the  preliminary  experiments  with  the 
l:30-scale  tow  and  towboat  hulls  were  useful  in 
delineating  the  general  features  of  tow  transit 
through  ice-covered  channels,  further  experiments 
with  larger  hulls  fitted  with  propellers  were  re¬ 
quired  to  produce  more  detailed  and  quantitative 
information.  For  example,  a  drawback  of  the  1  .-30- 
scale  towboat  was  that,  without  propellers,  its 
rounded,  flat-bottomed  hull  produced  a  track  filled 
with  dubiously  long  strips  of  broken  ice,  as  is 
evident  in  Figure  31 . 

Deep-channel  experiments 

Presented  here  are  observations  and  data  from 
the  experiments  conducted  with  the  1 : 1 5-scale  tow 


a.  n  =  2. 

Figure  31 .  Ice  formation  in  the  ice  tank  transited  by  the 
l:30-scale  towboat;  ff  =  1/(2  hours). 
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b.  n  =  6. 


c.n  =  12. 


Figure  31  (Continued).  Ice  formation  in  the  ice  tank  transited  by  the  l:30-scale  towboat;  f(  =  2/(2  hours). 


Figure  32.  Ice  thickness  vs  cumulative  cooling  time;  l:30-scale  towboat;  ff  =  1/(2  hours). 
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and  wedge  hulls  transiting  the  ice  tank  at  the  deep- 
channel  condition;  the  ratios  of  water  depth  to 
vessel  draft  yJD  were  5.62  and  3.75  for  the  tow  and 
wedge  hulls,  respectively.  As  listed  in  Table  1  the 
transit  frequencies  investigated  with  both  the  tow 
and  wedge  hulls  wer ef  =  1/(2  hours),  1/(4  hours) 
and  1/(8  hours).  An  additional  set  of  experiments 
was  conducted  with  transits  scheduled  at  f  =  1  / 
(0.5  hour)  and  fn  =  1  /(7.5  hours).  The  wedge  was 
used  for  two  further  experiments  at  /  =  l/(0.5 
hour)  and  1/(0.17  hour).  The  latter  experiment 
was  performed  for  a  short  duration  in  order  to 
assess  the  variation  of  (3  with  the  number  of  tran¬ 
sits.  Full-scale  equivalent  frequencies  can  be  as¬ 
sessed  using  eq  31.  Included  in  the  results  pre¬ 
sented  for  each  experiment  are  the  following: 

•  A  brief  illustrated  description  of  the  salient 
features  of  ice  formation; 

•  The  effect  of  £Sd  on  the  thicknesses  of  the 
brash-ice  layer  over  the  track  qb,  the  ridge  q^ 
the  plug  (for  the  tow  hull)  q  and  the  ice  sheet 

b;  P 

•  The  effect  of  £Sd  and  number  of  transits  n  on 
the  areal  concentration  of  open  water  P  and 
the  porosity  p  of  brash-ice  accumulations; 
and 

•  The  size  distribution  of  brash  ice  in  the  chan¬ 
nel. 

Each  data  point  in  plots  of  qb,  qrand  q  vs  ESd 
represents  a  transit.  The  brash-ice  sizes,  given  in 
terms  of  dQ  and  og,  are  summarized  in  Table  2. 


Photographs  of  brash-ice  samples  taken  at  the  end 
of  each  experiment  are  presented  in  Appendix  B. 
The  volumes  of  ice  grown  along  each  track  tran¬ 
sited  by  the  wedge  hull,  and  the  thickness  of  brash- 
ice  accumulation  in  each  track,  both  normalized 
with  q,  are  presented  in  Appendix  C. 

Transits  at  f(  =  1/(0. 5  hours) 

Figure  33  illustrates  the  formation  of  an  ice- 
filled  track  along  the  laboratory  channel  transited 
at  0.5-hour  intervals  by  the  wedge  hull.  Ice  formed 
along  the  track  in  approximately  the  same  se¬ 
quence  as  illustrated  in  Figure  14.  The  first  few 
transits  left  brash  ice  strewn  along  the  track  at  high 
areal  concentrations.  Most  of  the  ice  was  cast  side¬ 
ways  beneath  the  ice  sheet,  and  no  ice  plug  formed 
at  either  end  of  the  channel.  Further  transits  by  the 
wedge  hull  led  to  increased  coverage  of  the  track 
with  brash  ice,  and  eventually  the  track  became 
fully  covered  by  a  brash-ice  layer  that  thickened 
more  or  less  continuously  with  subsequent  tran¬ 
sits.  Values  of  qb,  qr  and  q  increased  monotoni- 
cally,  though  somewhat  erratically,  with  increasing 
ISdand  n  (Fig.  34),  reflecting  the  nonuniformity  of 
the  accumulations.  Also,  as  shown  in  Figure  33d, 
the  layer  of  brash  ice  accumulated  in  the  track  was 
occasionally  marked  by  a  furrow  [similarly  the  ice 
thickness  is  much  less  at  the  center  of  the  profile 
measured  by  Kannari  (Fig.  9)].  The  rate  at  which 
the  track  became  covered  with  brash  ice  can  be 
determined  from  Figure  35,  which  shows  that  the 


Table  2.  Summary  of  brash-ice  dimensions. 


Hull 

form 

(hours) 

y„ 

(m) 

d 

c 

(nun) 

Wedge 

0.5 

1.07 

15 

2.3 

2.0 

1.07 

93 

2.5 

4.0 

1.07 

135 

1.5 

8.0 

1.07 

124 

2.4 

Tow 

2.0 

1.07 

76 

3.7 

4.0 

1.07 

179 

1.5 

8.0 

1.07 

149 

1.7 

4.0 

0.53 

69 

1.9 

2.0 

0.22 

24 

3.8 

4.0 

0.22 

36 

4.5 

Wedge 

* 

1.07 

66 

3.9 

Tow 

* 

1.07 

48 

5.9 

'Multifrequency  transits \n  =  4,  /(=  1  /  (0.5  hour),  /  =  1/(75  hours)] 
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Figure  34.  Ice  thickness  vs  cumu¬ 
lative  cooling  time ;  l:15-scale 
wedge  hull,  f(  =  l/(0.5  hour);  the 
total  number  of  transits  is  98. 
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Figure  36.  Size  distribution  of 
brash  ice;  1 :15-scale  wedge  hull;  f 
=  1/(05  hour). 
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areal  concentration  of  open  water  p  steeply  de¬ 
clined  with  increasing  number  of  transits  n. 
Throughout  this  and  the  other  experiments,  ap¬ 
proximately  10%  of  the  track  remained  open  wa¬ 
ter.  As  shown  in  Figure  35,  the  layer  was  very 
porous,  with  values  of  porosity  p  generally  rang¬ 
ing  from  about  65  to  45%. 

On  completion  of  the  experiment,  when  ri  =  60 
mm  (0.9  m  at  full  scale)  and  n  =  98,  the  brash-ice 
layer  had  an  average  thickness  such  that  T)b/r) =2.5 
to  3  (App.  C).  The  ridges  adjoining  the  track  at¬ 
tained  an  average  thickness  such  that  r|r/r|  s  3.5. 
The  estimates  (in  Appendix  C)  of  the  volume  of  ice 
grown  per  unit  area  of  the  track,  or  the  equivalent 
to  ice  thickness  expressed  as  r\.Jr\,  was  about  2. 

The  transits  produced  comparatively  small  ice 
pieces  that  resulted  from  broken  regrowth  ice  and 
from  ice  that  was  crushed  or  sheared  from  the  ice- 
sheet  edges.  Figure  36  shows  that  at  least  90%  of 
the  brash  ice  in  the  layer  consisted  of  fine  ice  pieces 
with  d  less  than  10  mm  (0.150  m  at  full  scale).  The 

e 

remaining  10%  of  the  brash  ice  comprised  pieces 
with  deup  to  65  mm  (about  1  m  at  full  scale,  for  X( 
=  15).  These  larger  sizes  were  formed  from  smaller 
pieces  that  had  fused  together  and  grown  during 
periods  of  regrowth.  During  this  experiment  no 
large  pieces  of  the  ice  sheet  were  broken  and 
drawn  into  the  channel. 

Transits  at  ff  =  1/(2  hours) 

Figures  37  and  38  illustrate  ice  formation  in  the 
channel  when  transited  every  two  hours  by  the 
wedge  and  tow  hulls,  respectively.  Figure  39  gives 
values  of  qb  and  qr  for  thickness  vs  ZSd- 

Frequent  transits  of  the  channel  by  the  wedge 
hull  resulted  in  practically  the  same  pattern  of  ice 
formation  as  portrayed  in  Figure  14  and  as  de¬ 
scribed  for  /t=  l/(0.5  hour).  A  layer  of  brash  ice 
covered  the  track,  then  thickened,  as  did  the  ad¬ 
joining  ridges,  with  increasing  number  of  transits. 
At  the  end  of  the  experiment,  T\b/ r\  and  T|r/r)  were 
both  approximately  1 .7  to  2.3,  with  q. /rj  s  1 .5  to  1 .7 
(Fig.  39a,  App.  C). 

Figure  40  shows  the  decline  of  p  with  increasing 
number  of  transits,  and  itshows  that,  once  formed, 
the  layer  of  ice  pieces  was  comparatively  porous, 
with  values  of  p  of  about  30-65%.  A  notable  differ¬ 
ence,  made  evident  by  comparing  Figures  37  and 
33,  is  that  transits  at  f  =  1/(2  hours)  produced 
larger  ice  pieces  than  did  transits  at  0.5-hour  inter¬ 
vals.  Because  individual  ice  pieces  were  larger, 
correspondingly  larger  variations  occurred  in  qb 
and  r|r  (Fig.  39a).  For  /( =  1/(2  hours),  de=  93  mm 
(1.40  m  at  full  scale),  with  og=  2.5.  The  ice  pieces 


b.n  =  4. 


Figure  37.  Ice  formation  in  the  ice  tank  transited  by  the 
l:15-scale  wedge  hull;  ff  =  1/(2  hours). 
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e.  n  =  22.  /.  n  =  30. 

Figure  37  (Continued).  Ice  formation  in  the  ice  tank  transited  fry  the  l:15-scale  wedge  hull;  f(  =  Zf(2  hours). 
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e.  n  =  24. 


/.  n  =  29. 


Figure  38.  Continued. 
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a.  l:15-scale  wedge  hull. 


CUMULATIVE  COOLING  TIME  (°C-DAYS) 


b.  l:15-scale  tow  hull. 

Figure  39.  Ice  thickness  vs  cumulative  cooling  time;  f(  =  1/(2  hours). 
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Figure  40.  Variation  of  areal  concentration  of  open  water  (3  and  porosity  p  with 
number  of  transits  n;  ff  =  1/(2  hours). 


were  typically  formed  of  broken,  platy  pieces  of 
regrowth  ice  fused  together. 

Frequent  transits  of  the  channel  by  the  tow  hull 
resulted  in  the  same  manner  of  ice  formation  indi¬ 
cated  in  Figure  25.  During  the  first  six  to  ten 
transits,  the  tow  broke  and  reamed  a  straight- 
edged  track  (Fig.  38a,b)  that  retained  small  pieces 
of  regrowth  ice  at  high  areal  concentrations.  Most 
of  the  brash  ice  was  either  conveyed  to  a  plug  that 
formed  at  one  of  the  channel,  or  it  was  displaced 
sideways  beneath  the  adjoining  ice  sheet.  During 
b  ar  transits  the  tow  hull  rode  against  and  broke 
sections  of  the  ice  sheet  bordering  the  track,  caus¬ 
ing  relatively  large  pieces  of  ice  sheet,  and  portir  = 
of  the  ridge,  to  enter  the  track  (Fig.  38c).  Ai 
several  more  transits  the  major  portion  of  the  ice  in 
the  track  was  cleared  and  shoved  to  the  plug  (Fig. 
38d)  or  displaced  beneath  the  ice  sheet.  The  broken 
section  of  ice  edge  became  smoothed  and  re-reamed 
during  further  periods  of  regrowth  and  transit 
(Fig.  38d,e).  However,  other  portions  of  ice  sheet 
were  broken  and  resulted  in  more  broken  ice 
bursting  into  the  track  (Fig.  38e).  Still  further  peri¬ 
ods  of  regrowth  and  transit  enabled  the  edges  of 
the  track  to  become  smoothed  (Fig.  38f).  Conceiva¬ 
bly  a  track  could  experience  several  cycles  of  this 
process  of  sheet  breaking  and  self-healing,  until 
the  edge,  reinforced  by  its  ridge,  could  not  be 
flexurally  broken  by  a  tow;  then  tows  would  be 
unable  to  transit  the  channel. 

The  fluctuations  in  rj.  and  qr  shown  in  Figure 
39b  coincide  with  periods  when  pieces  of  ice  sheet 


and  segments  of  ridge  entered  the  channel  and 
were  eventually  cleared  from  it.  The  wedge  hull 
produced  a  single  population  of  brash-ice  sizes  in 
its  track,  but  the  tow  hull,  by  mixing  regrowth  ice 
and  pieces  of  broken  ice  sheet,  produced  at  least 
two  populations  of  brash-ice  size.  The  mixture  of 
the  two  populations  led  to  a  nonuniform  distribu¬ 
tion  of  brash-ice  sizes,  as  shown  in  Figure  41;  d  = 
76  mm  (1.14  m  at  full  scale)  and  ag  =  3.7. 

Transits  at  ff  =  1/(4  hours) 

Ice  formation  and  tracks  formed  by  the  wedge 
hull  and  the  tow  hull  transiting  the  laboratory 
channel  at  a  frequency  of  ft=  1/(4  hours)  are 
shown  in  Figure  42.  Figure  43  presents  values  of  tlb 
and  T|  for  wedge  and  tow  transits  at  4-hour  inter¬ 
vals.  figure  44  presents  values  of  p  and  p. 

For  the  wedge  hull  the  reduction  in  transit 
frequency  from  1  /(0.5  hour)  to  1  /(4  hours)  led  to 
a  mild  decrease  in  the  thickness  of  the  brash-ice 
layer  formed  over  the  track  and  accumulated  as 
ridges  next  to  it.  A  reduction  in  the  volume  of  ice 
formed  is  expected  in  accordance  with  the  formu¬ 
lation  presented  earlier.  However,  the  compara¬ 
tively  slight  decrease  in  regrowth  with  decreasing 
frequency  (by  a  factor  of  16)  is  somewhat  unex¬ 
pected,  especially  in  light  of  what  current  litera¬ 
ture  might  lead  one  to  expect  (e.g.  Greisman  1981 ). 
Transits  of  the  wedge  hull  resulted  in  a  brash  ice 
layer  for  which  T|.  /q = 2.0  to  2.5  and  T| .  /t|  =1 .5  to  1 .7 
(App.  C).  The  ridge  thickness  was  such  that  T|b/T| 
=  2.3  to  2.6. 
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Figure  41.  Size  distribution  of  brash  ice;  f(  =  1/(2  hours). 


a.  1:1 5 -scale  wedge  hull  (n  =  18). 


b.  l:15-scale  tow  hull  (n  =  14). 


Figure  42.  Ice  formation  in  the  ice  tank  transited  at  f(  =  lj(4  hours). 
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a.  l:15-scale  wedge  hull. 


b.  l:15-scale  tow  hull  (n  =  14). 


Figure  43.  Ice  thickness  vs  cumulative  cooling  time;  ff  =  1/(4  hours). 
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Figure  44.  Variation  of  areal  concentration  of  open  ivater  p  and  porosity  p  with 
number  of  transits  n;  =  1/(4  hours). 
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Figure  45.  Size  distribution  of  brash  ice;  f(  =  1/(4  hours). 


Reductions  in  the  frequency  of  tow  transit  from 
1/(2  hours)  to  1  /  (4  hours)  resulted  in  only  slight 
changes  to  the  manner  by  which  ice  formed  in  the 
tow's  track,  and  they  did  not  result  in  dramatic 
variations  in  the  thickness  of  ice  accumulated  as  a 
layer  in  its  track  or  in  ridges  flanking  the  tracks. 
The  accumulations,  however,  were  less  uniform  in 
thickness  for  the  lesser  transit  frequency.  At  the 
end  of  tht  tow-hull  experiment,  q.  /q  =  2.3  (at  the 
monitoring  area,  but  elsewhere  qb/T|  =  1 .5)  and  r\r/ 
q  =  3.0.  For  these  two  patterns,  ice  formation  was 
dominated  by  the  facility  of  the  tow  to  move  ice 
along  the  channel.  However,  less  frequent  transits 
led  to  less  movement  of  ice  along  the  channel, 
proportionately  less  plug  development,  and  larger 
pieces  of  brash  ice  in  the  track.  Fewer  occasions  to 
clear  ice  and  therefore  longer  periods  of  regrowth 
produced  thicker,  more  rounded  (by  regrowth)  ice 
pieces  (compare  Fig.  38  and  42)  that  were  less 
easily  shoved  by  the  tow;  they  more  readily  rolled 
around  and  beneath  it.  Comparison  of  Figures  45 
and  41  shows  that,  for  tow  transits  at  a  frequency 
of  1  / (4  hours),  de=179  mm,  which  is  about  two  to 
three  times  the  value  of  de resulting  from  a  transit 
frequency  of  1/(2  hours).  Also,  the  lesser  transit 
frequency  resulted  in  more  uniform  brash  ice. 

Transits  at  f  =  1/(8  hours) 

The  condition  of  the  channels  transited  by  the 
wedge  and  tow  hulls  at  a  frequency  of  1/(8  hours) 
is  shown  in  Figure  46.  Figure  47  presents  values  of 
q,  q.  and  qrvs  ISd-  Transits  by  the  wedge  and  tow 
hulk  at  a  frequency  of  1  /(8  hours)  had  about  the 
same  effect  on  ice  formation  as  at  a  frequency  of  1  / 


(4  hours).  Figures  48  and  49  present  values  of  P,  p 
and  de  for  these  experiments.  At  the  end  of  the 
experiment,  qb/q  and  q./q  =  1.8  to  2.0  and  q./q  = 
1.2  to  1.3. 

Transits  at  ff  =  1/(0. 5  hour)  and  fn  =  1/(7. 5  hours) 

Figures  50  and  51  show  the  main  features  of  ice 
formation  in  the  channel  when  transited  by  the 
tow  and  wedge  hulls  in  a  multifrequency  transit 
schedule  involving  sequences  of  four  transits  at  f 
=  l/(0.5  hour),  with  the  sequences  having  a 
frequency  fn  =  1  /(7.5  hours).  (At  full  scale,  one 
could  envision  four  transits  every  8  hours  with  3-4 
day  lapses  between  transits.)  Values  of  qb,  qrand 
q  vs  ISd  are  presented  in  Figure  52  for  the  wedge 
and  tow  hulls. 

These  experiments  show  that  when  a  schedule 
of  transits  involves  sequences  of  transits  at  rela¬ 
tively  high  frequency,  then  displacement,  move¬ 
ment  and  clearance  of  brash  ice  are  important.  For 
transits  by  the  wedge  hull,  ice  was  displaced  side¬ 
ways  beneath  the  bordering  ice  sheet,  giving  the 
track  a  furrowed  appearance  as  shown  in  Figure 
50.  For  transits  by  the  tow  hull,  ice  was  shoved  and 
dragged  toward  plugs  formed  at  either  end  of  the 
track.  Figure  51  illustrates  the  clearance  of  ice 
during  three  consecutive  transits  by  the  tow  hull. 

Displacement  of  ice  from  the  tracks  formed  by 
the  wedge  hull  caused  cyclic  fluctuations  in  the 
recorded  values  of  qb  and  qK  as  shown  in  Figure 
52.  Typically  the  first  transit  of  a  series  reopened 
the  track  and  displaced  broken  ice  sideways  to  the 
ridge  beneath  the  ice  sheet.  Some  of  the  displaced 
ice  was  drawn  back  into  the  track.  Subsequent 
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b.  l:15-scale  tow  hull. 

Figure  46.  Ice  formation  in  the  ice  tank  transited  at  ff  =  1/(8  hours). 
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THICKNESS  (MM)  THICKNESS  (MM) 


CUMULATIVE  COOLING  TIME  (°C-DAYS) 


a.  l:15-scale  wedge  hull. 


CUMULATIVE  COOLING  TIME  (°C-DAYS) 


b.  l:15-scale  tow  hull. 

Figure  47.  Ice  thickness  vs  cumulative  cooling  time ;  f,  =  1/(8  hours). 
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'ure  48.  Variation  of  areal  concentration  of  open  water  p  and  porosity  p  with 
mber  of  transits  n;  f(=  1(8  hours). 


Figure  49.  Size  distribution  of  brash  ice;  f(  =  1/(8  hours). 
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b.  n  =  35. 

Figure  50.  Ice  formation  in  the  ice  tank  transited  by  the  l:15-scale  tow  hull;  ff  =  1/(0. 5  hour)  and  fn  =  1/(7. 5  hours), 
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a.  1:1 5-scale  wedge  hull. 


b.  l:15-scale  tow  hull. 

Figure  52.  Ice  thickness  vs  cumulative  cooling  time;  ff  =  1/(0. 5  hour)  and  f  =  1/(7. 5  hours). 


transits  in  quick  succession,  and  without  much 
opportunity  for  regrowth,  displaced  more  broken 
ice  to  the  ridges,  causing  them  to  enlarge  and 
thinning  the  layer  in  the  track  as  well  as  creating  a 
furrow  along  its  center  line.  This  process  is  evident 
in  Figure  52a  for  transit  sequences  until  ISd  = 
12°C-days.  Thereafter,  it  is  obscured  in  the  data 
because  of  poor  accuracy  of  measurement;  the 
thickness  of  individual  ice  pieces  was  a  significant 
proportion  of  the  overall  thickness.  Throughout 
the  experiment  the  track  filled  with  brash  ice  such 
that,  on  average,  q.  /q  s  2.  The  series  of  three 
transits  at  /j  =  1  / (0.5  hour)  depicted  in  Figure  51  is 
representative.  The  first  transit  following  a  6-hour 


interval  between  transits  (Fig.  51a)  left  the  track 
strewn  with  broken  ice.  The  two  _  ubsequent  tran¬ 
sits  at  0.5-hour  intervals  (Fig.  5lb,c)  cleared  much 
of  the  broken  ice  from  the  channel,  shoving  it  into 
the  plug  or  displacing  it  to  the  ridges.  Data  on  |3,  p 
and  de  are  presented  in  Figures  53  and  54. 

Shallow-channel  experiments 

The  three  experiments  conducted  with  the  tow 
hull  transiting  shallow,  ice-covered  channels  (y  / 
D  =  1.16  and  2.78)  indicate  that  ridges  and  plugs 
develop  more  rapidly  than  they  do  during  fre¬ 
quent  transits  of  deep  channels.  The  experiments 
also  demonstrated  that  ridge  and  plug  formation 
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Figure  53.  Variation  of  areal  concentration  of  open  water  p  and  porosity  p  with 
number  of  transits  n;  ff =  l/(0.5  hour)  and  fn  =  1/(7. 5  hours). 


Figure  54.  Size  distribution  of  brash  ice;  f(  =  1/(05  hour)  and  (n  =  1/(75  hours). 


pose  potentially  severe  problems  for  navigation  of 
ice-covered  shallow  channels  because  they  may 
extend  over  the  full  depth  of  the  channel,  as  illus¬ 
trated  in  Figure  26. 


Transits  at  f  =  1/(2  hours)  and  1/(4  hours) 

Figures  55  and  56  show  the  shallow  channel 
(0.22  m  deep)  after  being  transited  by  the  tow  hull 
at/,=  1  /  (2  hours)  and  1/(4  hours),  respectively.  It 
is  apparent  from  these  figures  that  little  ice  accu¬ 
mulated  in  the  tracks  opened  during  transiting. 
Figure  57,  which  shows  t|b,  r)rand  q  vs  ZSd,  shows 
that  little  brash  ice  remained  in  the  track.  The 
constricted  clearance  (only  0.03  m)  between  the 
tow  and  the  channel  bottom  forced  ice  to  be  dis¬ 


placed  sideways  to  the  ridges  or  shoved  along  the 
track.  Figure  57  indicates  the  resulting  rapid  growth 
of  the  ridges,  especially  when  compared  with  the 
rates  of  ridge  growth  indicated  in  Figures  39b  and 
43b.  As  the  ridges  approached  the  channel  bottom, 
less  ice  could  be  displaced  sideways  and  more  was 
confined  in  a  jumbled  mass  shoved  ahead  of  the 
tow  hull,  as  shown  in  Figure  58.  The  mass  grew  as 
the  tow  moved  along  the  track,  collecting  addi¬ 
tional  ice.  The  ice  mass  was  shoved  to  plugs  at 
either  end  of  the  track.  While  shoving  ice,  the  tow 
overrode  and  frequently  grounded  on  ice.  The 
frequency  of  grounded  increasing  as  the  ridges 
neared  the  channel  bottom. 

An  obvious  consequence  of  increasing  confine- 
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a.n-1. 


b.n  =  3. 


Figure  56.  Ice  formation  in  the  ice  tank  transited  by  the  l:15-scale  tow  hull;  ff  =  1/(4  hours);  yo/D  =  1.16. 
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a.  f t  =  1/(2  hours). 


b.  ff  =  1/(4  hours). 

Figure  57.  Ice  thickness  vs  cumulative  cooling  time;  y0/D  =  1.16. 


ment  of  ice  to  the  track  was  increasing  the  resis¬ 
tance  encountered  by  the  tow  hull.  Although  the 
resistance  was  not  measured,  the  levels  of  resis¬ 
tance  were  obvious  from  the  vibrations  of  the 
towingarmsand  the  noise  generated  during  break¬ 
ing  and  shoving  of  ice.  The  resistance  increased 
noticeably  when  the  ridges  extended  about  three- 
quarters  of  channel  depth  and  when  the  ice  sheet 
thickness  was  of  similar  magnitude  to  the  clear¬ 
ance  beneath  the  tow.  From  then  on,  ice  was  in¬ 
creasingly  confined  and  shoved  ahead  of  the  tow 
or  forced  sideways  to  the  ridges,  which  thickened 


rapidly  until  they  reached  the  channel  bottom  or 
were  overridden  by  the  tow.  During  such  transits 
the  tow  exerted  considerable  force  against  the 
bordering  ice  sheet,  occasionally  breaking  large 
cusps  (underpinned  by  ridge),  as  shown  in  Figure 
55a.  Further,  when  the  tow  and  the  ice  beneath  it 
occasionally  ground  against  the  channel  bottom,  it 
would  bend  the  aluminum  false  floor  that  was 
used  to  simulate  part  of  the  channel  bed.  At  full 
scale  this  grounding  would  gouge  the  channel  bed 
and  possibly  damage  the  tow.  The  level  of  resis¬ 
tance  was  such  that,  after  about  a  dozen  transits  (at 
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Figure  58.  l:15-scale  tow  hull  shoving  confined  broken  ice  while  transiting  the  shallow 
channel;  yo/D  =  1.16. 


Figure  59.  Size  distribution  of  brash  ice;  yo/D  =  1.16. 
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either  frequency),  a  full-scale  tow  would  find  the 
channel  unnavigable. 

Brash  ice  remaining  in  the  shallow  tracks  was 
considerably  finer  in  size  than  the  brash  ice  in 
tracks  through  the  deep  channel,  as  can  be  seen  by 
comparing  values  of  dein  Figures  59  and  41  and 
Table  2.  The  average  sizes  of  brash  ice  in  the 
shallow  channel  were  about 20-30%  of  those  in  the 
deep  channel.  Confined  clearance  between  the 
tow  underside  and  the  channel  bottom  screened 
large  ice  from  passing  beneath  the  hull.  Addition¬ 
ally,  ice  in  the  shallow  channel  underwent  much 
more  abrasion  and  breaking  than  did  ice  in  the 
deep  channel. 

Analysis  of  results 

The  major  value  of  the  ice-tank  experiments  is 
that  they  enabled  detailed  observation  and  de¬ 
scription  of  ice-formation  processes  in  a  frequently 
transited,  ice-covered  channel.  Before  now,  many 
of  these  processes  were  either  unknown  or  inade¬ 
quately  understood  and  documented.  Quantita¬ 
tive  evaluation  of  the  data,  notably  in  the  accurate 
evaluation  of  volumes  and  rates  of  ice  growth,  was 
complicated  by  substantial  variations  in  the  way 
that  brash  and  broken  ice  accumulated  in  tracks 
cut  through  the  channel,  especially  those  created 
during  transits  by  the  1:30- and  the  1:1 3-scale  tow 
hulls.  Consequently  the  numerical  model  could 
only  be  calibrated  with,  and  tested  against,  the 
results  obtained  from  experiments  conducted  with 
the  wedge  hull.  Those  experiments,  however,  were 


also  fraught  with  variations  in  the  pattern  of  ice 
formation,  though  not  to  the  same  extent  as  for  the 
tow  hull.  The  results  of  the  ice-tank  experiments 
suggest  that  precise  simulation  of  ice  formation  in 
frequently  transited  channels  is  complicated  by 
the  occasionally  irregular,  or  stochastic,  manners 
by  which  brash  ice  forms  and  is  redistributed  in 
vessel  tracks.  For  example,  large  fragments  of 
sheet  ice  bordering  a  track  may  drift  into  it  and 
temporarily  alter  the  local  pattern  of  ice  regrowth 
and  brash-ice  accumulation. 

The  experiments  confirm  the  conclusion  tenta¬ 
tively  drawn  earlier  that  hull  characteristics,  espe¬ 
cially  bow  form  and  hull  length,  significantly  in¬ 
fluence  ice  formation  in  navigation  channels,  and 
that  they  do  so  by  affecting  the  patterns  of  brash- 
ice  accumulation.  The  disparate  patterns  of  ice 
formation  produced  by  ships  and  ice-entrapping 
tows,  two  very  different  hull  forms,  are  summa¬ 
rized  in  Figures  14  and  25,  respectively.  The  essen¬ 
tial  feature  of  transiting  that  causes  the  dissimilar 
patterns  is  the  propensity  of  tows  to  entrap  and 
move  brash  ice  along  the  tracks  that  they  open.  As 
a  consequence,  ice  plugs  form  as  thick  local  accu¬ 
mulations  of  brash  ice.  The  two  patterns  would  be 
more  similar  if  the  tow  hull  form  were  modified  by 
ice  accumulation  so  as  to  become  rounded  and 
entrap  less,  or  no  additional,  ice. 

In  the  ensuing  discussion  the  influences  of  tran¬ 
sit  frequency /t  on  ice  formation  are  evaluated  first. 
The  influences  of  water  depth  on  ice  formation  are 
then  appraised.  Finally,  an  attempt  is  made  to 
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relate  the  areal  concentration  of  open  water  (1  to 
the  number  of  transits  and  the  transit  frequency. 
Also  assessed  is  the  variation  of  ice-displacement 
proportion  as  a  brash-ice  layer  forms. 

Influence  of  transit  frequency 

For  the  range  of  frequencies  examined,  increas¬ 
ing  transit  frequency  /.  led  to  increased  volume  of 
brash-ice  growth  in  the  ice  tank.  The  equivalent 
thickness  of  solid  ice  grown  in  tracks  transited  by 
the  wedge  hull,  normalized  with  ice  sheet  thickness, 
q./q,  is  shown  in  Figure  60  vs  fa  q.is  the  volume  of 


brash  ice  per  unit  surface  area  of  track,  with  the 
track  width  taken  as  being  equivalent  to  the  hull 
beam.  Each  value  of  q./q  is  estimated  using  an 
average  cross-sectional  profile  and  an  average 
value  of  porosity  of  accumulated  brash  ice.  The 
maximum  value  of  qr/q  attained  in  the  ice  tank  is 
about  2.  For  the  wedge  hull  the  thickness  of  brash- 
ice  layers  and  ridges,  expressed  as  qb/q  and  qr/q, 
increased  mildly  with  increasing  f^  as  shown  in 
Figure  61a.  The  propensity  of  the  tow  hull  to  clear 
ice  from  its  track  resulted  in  somewhat  different 
trends,  as  shown  in  Figure  61b.  The  ratio  qb/q  did 


a.  l:15-scale 
wedge  hull. 


b.  l:15-scale 
tow  hull. 
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MEDIAN  DIMENSION,  d  (MM) 


CUMULATIVE  COOLING  TIME  (°C-DAYS) 


Figure  63.  Variation  of  median  plan  dimensions  of  brash  ice  with  cumulative  cooling  time. 
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not  vary  significantly  with  ft  for  the  tow  hull.  If 
anything,  qb/r)  decreased  slightly  with  increasing 
ft.  However,  for  the  tow  hull,  T|b/q  did  increase 
with  increasing  / .  For  both  hulls  the  error  margin 
in  measuring  t|b/q  and  T|r/r|  is  on  the  order  of  the 
mean  size  of  brash-ice  pieces  comprising  layers  or 
ridges.  For  both  hulls  the  ridges  thickened  more 
rapidly  with  increasing  /  than  did  the  layers  be¬ 
cause  more  ice  was  displaced  from  the  channel. 

Transit  frequency  influenced  the  mean  size  of 
the  brash  ice  formed  during  transits  by  either  the 
wedge  and  the  tow  hulls.  For/tgreater  than  about 
1/(4  hours).  Figure  62  and  Table  2  show  that  de 
decreased  almost  linearly  with  increasing/,  at  least 
for  the  range  of /  used  in  the  experiments.  The  few 
data  corresponding  to/t  less  than  about  1/(4  hours) 
suggest  that  deis  likely  to  be  independent  of/(.  The 
mean  plan  dimensions  of  brash  and  broken  ice 
were  measured  and  recorded  versus  ESd  during 
two  experiments  with  the  tow  hull,  those  for  which 
/t=  1/(2  hours)  and  1/(4  hours).  The  data,  pre¬ 
sented  in  Figure  63,  show  that  the  mean  plan  size 
of  brash  ice  decreased  initially  but  thereafter  at¬ 
tained  more-or-less  constant  values  with  increasing 
ZSd-  This  result  is  similar  to  that  reported  by 
Sandkvist  (1 980)  for  the  variation  of  brash-ice  thick¬ 
ness  with  number  of  transits  (Fig.  11). 

Influence  of  channel  depth 

Figures  25  and  57  show  that,  for  transits  of 
shallow  channels,  ridges  and  plugs  develop  more 


rapidly  for  smaller  values  of  y  /D.  Also,  as  is 
evident  in  Table  2  and  Figure  6l,  finer  brash  ice 
accumulated  in  the  shallower  channels.  The  rela¬ 
tionship  between  d  and  y  /D  is  almost  linear  for 
the  data  shown.  Onfy  smaller  ice  pieces  could  pass 
through  the  diminished  clearance  between  the 
hull  and  the  channel  bottom.  Larger  ice  pieces 
were  either  displaced  to  ridges  along  the  track  or 
were  shoved  to  plugs  in  the  track.  Additionally, 
transit  of  shallower  channels  resulted  in  consid¬ 
erably  more  abrasion  and  grinding  of  ice  pieces 
than  for  the  deeper  channels,  thereby  contributing 
to  the  greater  amount  of  fine  ice  in  the  track. 

Variations  of  P  and  e 

To  •  se  the  numerical  model  outlined  earlier, 
mathe  matical  expressions  must  be  developed  for 
the  \  ariations  of  P  and  £  throughout  a  navigation 
period.  This  task,  however,  is  made  complicated 
by  the  diverse  variables  that  may  influence  p  and 
£,  e.g.  hull  dimensions  and  shape,  hull  speed,  ice- 
sheet  thickness,  and  thickness  of  brash  ice  accu¬ 
mulated  in  a  vessel  track  and  ridges  flanking  it,  as 
well  as  the  element  of  randomness  that  pervades 
ice  formation  in  frequently  transited  channels. 
Nonetheless,  an  attempt  is  ventured  here  to  estab¬ 
lish  sets  of  curves  for  the  variations  of  the  ice-tank 
data  from  the  experiments  with  the  wedge  hull. 
Note  that  terms  P  and  e  are  not  necessarily  in  direct 
relationship  to  one  another,  because  ice  cleared 
from  the  surface  of  the  track  may  not  only  be 
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displaced  to  the  ridges,  but  it  may  collect  as  multi¬ 
layer  clusters  within  the  track  (Fig.  14). 

Figures  35,  40,  44  and  48  show  that  P  varied 
fairly  consistently  with  the  number  and  frequency 
of  transits  by  the  wedge  hull.  Generally,  lesser 
values  of  ft  resulted  in  initially  lower  values  of  p 
and  required  fewer  transits  before  the  track  be¬ 
came  fully  covered  by  a  layer  of  brash  ice.  This 
trend  becomes  more  apparent  when  values  of  P 
are  related  to  the  nondimensional  parameter 
a(«AT//t) 05  ID,  as  shown  in  Figure  65.  This  pa¬ 
rameter  is  significant  because  it  reflects  the  effect 
of  relative  magnitudes  of  ice-cover  thickness  as 
well  as  hull  draft  and  shape  on  Po,  the  initial  value 
of  P  after  the  first  transit  of  a  vessel.  Transits  of 
thicker  ice  covers  left  more  ice  in  the  track.  A 
family  of  curves  can  be  defined  in  accordance  with 
values  of  Po  and  the  decline  of  P  with  increasing  n. 
It  was  found,  as  shown  in  Figure  66,  that  the 
number  of  transits  required  for  a  virtually  full 
cover  to  develop  over  the  track  was  linearly  pro¬ 
portional  to  log  (/).  Actually,  throughout  the  ex¬ 
periments,  coverage  was  never  full.  There  always 
seemed  to  be  about  10%  of  the  track  that  was  open 
water. 

Variation  of  the  ice-displacement  parameter  e 
during  each  experiment  with  the  wedge  hull  was 
estimated  as  2AVr[(l— p/(r|.B)]  and  related  to  the 
parameter  (D'-t)  )/r\  ,  in  which  D'  -  [q.B/ 
({l-p}tan9w)j°  5,  with  6^  =  53°  (Fig.  67).  As  indi¬ 
cated  in  Figure  68,  the  latter  expression  is  intended 


to  take  into  account  the  relative  elevations  of  the 
brash-ice  layer  when  it  is  passed  around  the  hull 
and  the  ridge  thicknesses.  Similar  expressions  can 
be  devised,  provided  they  reflect  the  influences  on 
E  of  the  relative  thicknesses  of  the  ridge  and  the 
layer,  as  well  as  the  layer  thickness  relative  to  the 
hull  form  and  draft.  Generally,  e  increased  as  r|  / 
T]b  decreased  and  as  the  layers  thickened.  Ridges 
acted  to  partially  contain  such  ice  in  the  tracks.  The 
considerable  scatter  in  the  plot  of  e  versus  (O'— r\f)  / 
T]rcan  be  related  to  the  irregular,  and  noncontinu- 
ous,  manner  by  which  layers  and  ridges  actually 
thickened.  As  can  be  seen  from  Figure  34,  for 
example,  t|band  t|rdid  not  increase  smoothly  with 
either  time  or  number  of  transits.  No  influence  of 
/  or  n  was  evident  in  the  data. 

On  the  whole,  values  of  E  are  quite  small,  being 
typically  on  the  order  of  0.1  or  less.  The  curve  fitted 
as  an  envelope  to  the  data  is  used  for  selecting 
values  of  e  used  in  exercising  the  numerical  model. 

NUMERICAL  PREDICTIONS 

The  numerical  model  is  exercised  here  to  ex¬ 
tend  the  insights  obtained  from  the  ice-tank  ex¬ 
periments  and  to  obtain  a  broader  understanding 
of  ice  formation  in  frequently  transited  channels. 
The  model  is  used  first  to  determine  volumes  and 
distributions  of  brash  ice  produced  in  a  deep,  wide 
channel  transited  at  a  single,  sample  frequency  of 
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/  =  6  transits/ day  under  a  range  of  air  tempera¬ 
tures.  The  resulting  predictions  illustrate  the  in¬ 
formation  that  the  model  yields.  The  effects  of 
transit  frequency  and  scheduling  on  the  volume  of 
ice  grown  and  the  distribution  of  brash  ice  are 
considered  subsequently.  Investigated  finally  are 
the  influences  of  channel  geometry,  notably  chan¬ 
nel  width  and  depth,  on  ice  formation.  This  section 
finishes  with  some  reflections  on  the  validity  of  the 
numerical  predictions. 

Ice  formation  is  simulated  for  a  representative 
channel  that  is  transited  in  frigid  air  for  a  naviga¬ 
tion  period  of  30  days,  which  is  long  enough  so 
that  the  calculations  would  illustrate  the  salient 
features  of  ice  formation.  It  is  assumed,  somewhat 
conservatively,  that  no  snow  falls  during  this 
period.  The  channel  is  initially  taken  to  be  deep 
and  wide  so  as  to  exclude  the  influences  of  channel 
geometry  on  ice  formation.  All  predictions  are 
based  on  an  ice-growth  coefficient  of  a  =  24  mm/ 
(°C-day).  This  value  is  suggested  from  estimates 
made  for  ice  growth  on  the  Mississippi  River 
(Ashton  1974)  and  on  the  St.  Marys  River  (Vance 
1980).  Furthermore  it  is  comparable  to  the  value  of 
28  mm/(°C-day)  recommended  by  Michel  (1971) 
for  ice  growth  over  snow-free  lakes.  A  wide  river, 
such  as  the  Mississippi  River,  without  a  snow 
cover  can  be  likened  to  a  small  lake. 

As  stated  earlier  the  model  simulates  transits  by 
vessels  that  do  not  shove  or  move  significant 
amounts  of  ice  along  channels.  Therefore,  it  simu¬ 
lates  transits  by  tows,  or  other  bluff-bowed  flat- 
bottomed  vessels,  only  if  they  do  not  move  ice 
along  the  tracks.  This  assumption  may  be  reason¬ 
able  for  tows  that  either  have  a  ttained  equilibrium 
under-hull  accumulations  of  ice  or  are  moving 
fairly  rapidly  so  as  not  to  entrap  much  ice.  The  ice- 
tank  experiments  revealed  that  numerical  simula¬ 
tion  of  ice  formation  is  useless  if  tows  accumulate 
and  transport  ice  while  transiting  a  track. 

Values  for  the  parameters  f)  and  £  are  obtained 
from  curves  fitted  to  the  experimental  data  con¬ 
tained  in  Figures  65  and  67,  respectively.  When 
using  Figure  65  to  estimate  P  values,  experiment- 
scale  values  of  transit  frequency  f  are  related  to 
full-scale  values  through  the  geometric  scale  X,  = 
15,  in  line  with  the  reasons  explained  earlier.  In 
accordance  with  this  relationship,/^  1  / (0.5  hour), 
or  48/day,  at  the  ice  tank  is  taken  here  to  be 
equivalent  to  (48/  15)/day =3/day  at  full  scale.  All 
predictions  are  calculated  with  p  =  0.5,  which  is 
obtained  from  experimental  data.  As  channel  size 
(relative  to  vessel  size)  primarily  affects  the  space 
available  for  brash  ice  to  be  displaced  from  tracks 


opened  by  vessels  or  to  accumulate  as  ridges  flank¬ 
ing  tracks,  an  attempt  is  made  to  account  for  these 
influences  on  ice  formation  by  modifying  the  rela¬ 
tionships  for  ridge  geometry.  In  the  ice  tank,  ap¬ 
proximately  one  minute  was  required  for  water 
and  brash  ice  to  calm  sufficiently  after  a  transit 
such  that  ice  could  begin  to  regrow.  If  a  factor  of  1 5 
is  used  for  the  time  scale,  one  minute  translates  to 
0.25  hour  at  full  scale.  This  value,  which  is  used 
here  for  t^,  agrees  well  with  observed  times  for  ice 
regrowtn  to  commence  after  vessel  transit.*  A 
value  of  0f=  30°  obtained  from  experimental  data 
is  used  for  estimating  ridge  proportions. 

Influences  of  air  temperature  on  ice  formation 

The  influences  of  air  temperature  on  ice  forma¬ 
tion  in  frequently  transited  channels  are  examined 
here  by  estimating  ice  formation  under  three  tem¬ 
peratures  of  frigid  air:  Ta  =  -5,  -10  and  -20°C.  As 
noted  above,  all  predictions  are  based  on  an  as¬ 
sumed  value  of  a  =  24  mm/(°C-day).  The  tem¬ 
perature  at  which  water  freezes  is  taken  as  0°C. 
Transits  at  a  sample  frequency  of  /  =  6/day  are 
considered,  with  the  initial  ice  growth  period  P.  - 
P{  (Fig.  16),  and  navigation  period  Pn  =  30  days. 

The  volumes  of  ice  grown,  the  layer  thicknesses 
and  the  ridge  thicknesses  predicted  for  simulated 
transitsoftherepresentativechannelarepresented 
in  Figure  69  vs  time,  with  Taas  the  third  parameter. 
The  ice-cover  thickness  q  is  shown  only  for  Ta  = 
-10°C.  The  corresponding  ice-cover  thicknesses 
for  Ta  =  -5  and  -20°C  can  be  determined  by  using 
eq9  with  appropriate  factors.  For  these  figures  the 
ice  volumes  are  calculated  as  ice  volume  per  unit 
area  of  track  opened  by  transiting  vessels  q.;  the 
unit  area  of  track  is  taken  as  the  unit  length  of  track 
times  the  track  width,  which  here  is  taken  to  be 
equal  to  the  vessel  beam  B.  Figure  69  reveals  that, 
after  30  days  of  transiting  at  /  =  6/ day  with  Ta  = 
-10°C,  the  ice  cover  thickness  q  is  about  415  mm, 
q.  =  900  mm,  qr=  1000  mm  and  qr=  1 100  mm.  It  is 
also  evident  that,  if  normalized  with  q,  values  q;, 
qb  and  qr  are  not  fully  independent  of  Tfl.  The 
model  predicts  that  values  of  q;,  qb  and  qr  are  of 
nearly  comparable  magnitude  during  the  first  four 
days  of  transiting  when  Tg  =  -5  and  -10°C.  This 
somewhat  surprising  prediction  is  a  consequence 
of  the  relationship  adopted  for  (3,  the  proportion  of 
the  track  covered  by  open  waterduring  transiting. 
From  Figure  65,  it  follows  that  fewer  transits  are 
required  to  establish  a  full  layer  of  brash  ice  over  a 
track  transited  in  colder  air,  because  thicker  and 
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larger  ice  fragments  result  during  transiting.  As 
larger  ice  fragments  are  less  readily  displaced 
beneath  the  bordering  ice  sheet,  a  rapid  decrease 
in  P,  or  smaller  values  of  P,  result  and  less  ice  is 
produced.  The  example  presented  in  Figure  69 
reflects  the  sensitivity  of  the  predictions  to  values 
of  p. 


Transit  frequency  effects 


Consider  now  the  predicted  effects  of  transit 
frequency/,  on  ice  formation  under  a  single  value 
of  T  .  Figure  70,  which  generalizes  the  trends  pre¬ 
sented  in  Figure  69,  shows  the  effects  of  /  on 
volume  of  ice  grown,  normalized  as  q./q,  and  on 
the  thicknesses  of  brash  ice  accumulated  as  layers 
and  ridges  in  the  representative  channel,  normal¬ 
ized  as  qb/ q  and  qr/q,  respectively.  Values  of  r\  J 
q,  qb/q  and  qr/q  are  determined  for  a  navigation 
period  Pn  =  30  days  under  steady  conditions  of 


constant  T  =  -10°C,  with  a  =  24  mm/(°C-day). 
The  transit  frequency /  ranges  from  8/day  through 
1/(15  days)  to  zero. 


Though  the  volume  of  ice  produced  during 
frequent  transits  of  a  navigation  channel  increases 


with  increasing  duration  of  navigation  period  and 
with  increasing  /,,  values  of  q./q,  qb/q  and  qr/q 
follow  slightly  more  complex  trends.  For  each 
value  of/t,  q,/q  generally  increases  with  increas¬ 
ing  time  or  number  of  transits  (Fig.  70a).  Values  of 
q/q  increase  very  rapidly  during  initial  transits  at 
large  values  of  /,  because  of  the  initially  large 
values  of  P  for  larger  values  of  f  (as  indicated  in 
Fig.  65).  Once  P  values  have  decreased  ti  constant 


values  and  the  track  is  virtually  covered  by  a  layer 
of  brash  ice,  values  of  q,/q  increase  mildly  with 
increasing  time  or  n.  Larger  values  of  /  lead  to 
larger  values  of  q;/q  at  any  instant,  but  they  also 
cause  greater  rates  of  increase  of  qb/q.  The  vari¬ 
ation  of  qb/ q  with  duration  of  navigation  period 
(Fig.  70b)  is  similar  to  that  between  q./q  and  time 
or  number  of  transits.  Values  of  qb/q  initially 
increase  rapidly  with  time  as  a  full  layer  develops 


across  the  track,  and  the  rate  of  increase  is  larger 
for  larger  values  of/,.  Once  a  virtually  full  layer  has 
developed  across  a  track,  qb/q  increases  only 
mildly  with  increasing  time  or  n.  Values  of  q  /q 
(Fig.  70c)  begin  at  their  maximum  values  and  then 


gradually  diminish.  The  trends  for  q./q  and  qb/ q 
are  very  similar  to  those  from  the  experiments 
conducted  with  the  ice  tank  (the  pertinent  trends 
are  presented  as  Appendix  C). 


The  relationship  between  the  volumes  of  ice 
produced,  normalized  as  q./q,  and  the  transit 
frequency  /  after  30  days  of  transiting  under  air 


temperatures  of -5,  -10  and  -20°C  reveal  that  q./ 
q  increases  with  increasing  /  until  attaining  a 
maximum  value,  as  shown  in  Figure  71.  That  a 
maximum  value  for  q./q  should  occur  with  vari¬ 
ation  of  /,  is  evident  from  eq  19;  although  q./q 
increases  with  /.initially,  when  1  //,  decreases  to  t  , 
the  actual  time  of  ice  regrowth  is  decreased  such 
that,  eventually,  there  is  insufficient  time  for  ice  to 
regrow.  Here,  as  t  is  taken  to  be  0.25  hours,  no 
regrowth  occurs  when  /,  =  96  transits/day.  This 
frequency  is  an  order  of  magnitude  above  the 
frequencies  at  which  channels  are  likely  to  be 
transited.  Therefore,  for  practical  purposes,  it  can 
be  stated  that  increased  frequency  of  transiting 
generally  leads  to  increased  volume  of  ice  re¬ 
growth.  Figure  71  indicates  that  q./q  attains  a 
maximum  value  of  about  3,  which  occurs  when  / 
=  30-50  transits/day.  At  values  of  /  exceeding  1  / 
t ,  transiting  vessels  are  essentially  in  convoy  (if/, 
relates  to  transits  of  different  vessels,  rather  than  a 
single  vessel,  such  as  a  harbor  tugboat,  repeatedly 
transiting  a  track).  Transits  in  close  convoy  pre¬ 
vent  or  greatly  limit  regrowth  of  ice  over  a  track, 
with  the  possible  result  that  channel  water  may 
supercool  and  spawn  frazil  ice.  The  likely  effect  is 
channels  through  shallow  rivers.  The  present 
model  does  not  handle  frazil-ice  growth. 

Values  of  q./q  do  not  form  a  single  curve  for  the 
three  values  of  Ta  used  in  the  calculations.  By  and 
large,  lower  values  of  Ta  show  negligibly  larger 
values  of  q./q.  This  result  is  attributable  to  the 
larger  initial  values  of  P  that  result  when  using 
Figure  65  and  smaller  values  of  Ta. 

Effects  of  scheduled  transits 

Vessels  may  transit  navigation  channels  in 
complex  sequences  involving  several  transit  fre¬ 
quencies,  not  just  a  single  frequency  such  as  once 
or  a  dozen  times  daily.  A  general  schedule  is 
outlined  in  Figure  16.  It  may,  for  example,  consist 
of  daily  transits  that  are  limited  to  daylight  hours, 
during  which  transits  occur  at  some  constant  fre¬ 
quency.  Transits  also  may  occur  randomly.  A 
practical  motivation  for  examining  scheduled  tran¬ 
siting  is  to  determine  whether  there  are  any  sig¬ 
nificant  advantages  in  scheduling  transits  so  as  to 
minimize  ice  growth  due  to  transiting. 

The  numerical  model  was  exercised  to  ascertain 
the  effects  on  ice  formation  of  scheduled  transits 
involving  the  two  frequencies  /  and  /,  as  defined 
earlier.  For  a  total  period  PT  =  3(5  days,  the  transit- 
series  frequency  /  =  n/  day  was  varied,  with  n  =  2, 
4  and  6  transits.  Transits  are  taken  to  occur  in  a 
navigation  period  Pn  =  (n-1)  P,,  with  the  transit 
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period  P  =  0.5-12  hours;  the  period  between 
navigation  periods  Pq=  l(day)  -  P  Additionally, 
Ta  =  -10°C  and  a  =  24  mm/(°C-day). 

The  results  are  presented  in  Figure  72  for  fs=  2, 
4  and  6/ day.  Evidently  tv/t],  t|b/T)  and  rjr/r\  are 
mildly  affected  by  though  its  influence  increases 
with  increasing  n  and/s-  Smaller  amounts  of  ice  are 
produced  for  larger  values  of even  though  n  is 
constant.  At  larger  values  of  ft,  P  approaches  t 
such  that  the  net  period  of  ice  regrowth  decreases. 
A  limiting  condition  is  convoying  of  vessels  when 
P(<f .  Figure  72b,  for  example,  shows  that  a  daily 
sequence  of  four  transits,  0.5-hour  apart  results  in 
about  18%  less  ice  growth  than  does  a  daily  se¬ 
quence  of  four  transits  6-hours  apart.  The  com¬ 
mensurate  reduction  in  layer  thickness  predicted 
for  the  two  schedules  is  about  11%.  For  an  ice  cover 
thickness  of  about  416  mm,  the  additional  thick¬ 
ness  of  ice  growth  over  the  track  during  30  days  is 
less  than  about  75  mm.  The  practical  implication  of 
the  information  presented  in  Figure  72  is  that  there 
is  no  significant  advantage  in  scheduling  transits 
of  ice-covered  channels,  at  least  for  a  daily  sched¬ 
ule  involving  six  or  fewer  transits,  for  the  purpose 
of  minimizing  the  volume  of  additional  (or  brash) 
ice  growth. 


exceeds  1 0%  of  the  maximum  tow-barge  draft;  the 
minimum  navigation  depth  in  waterways  main¬ 
tained  by  the  Corps  of  Engineers  is  3.05  m  (10  ft), 
whereas  the  maximum  draft  of  tow-barges  is  2.75 
m  (9  ft).  Thus  far,  however,  the  numerical  model 
has  been  exercised  only  for  a  representative  chan¬ 
nel  that  was  wide  and  deep  enough  that  channel 
boundaries  did  not  have  to  be  taken  into  account. 
An  attempt  is  now  made  to  assess  the  effects  of 
channel  width  and  depth  on  the  volume  of  ice 
produced  by  frequent  transiting.  Also,  so  as  to  put 
the  effects  of  frequent  transiting  in  perspective,  the 
volumes  of  ice  produced  at  navigation  channels 
are  related  to  the  overall  volume  of  ice  produced  at 
channel  cross  sections. 

Channel  width 

The  approach  taken  here  to  evaluate  the  effects 
of  channel  width  on  ice  formation  is  to  relate  the 
ridge  geometry  to  the  available  space  adjacent  to  a 
track.  As  indicated  in  Figure  67,  the  ridge  width  VV. 
is  constrained  by  the  width  of  channel  such  that 
the  following  relationships  are  needed  to  determine 
W. and  the  ridge  thickness: 

Wr  =  2  r|r  tan  0r ,  if  VV  <  0.5  ( W-B)  (32) 


Effects  of  channel  geometry 

Many  channels,  especially  river  channels,  have 
reaches  that  are  narrow  and  shallow  compared  to 
vessel  size.  Approaches  to  locks,  for  instance,  can 
be  confined,  and  the  minimum  navigation  depth 
at  many  sections  of  navigable  U.S.  rivers  barely 


or 

0.5  (W-B) ,  if  2  t|r  tan  0r>  0.5  (W-B).  (33) 

It  is  assumed  that  the  curves  for  P  and  £  (Fig.  66, 68) 
remain  valid  for  estimating  (3  and  e. 
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a.  f.  =  2/day 


b.  fs  =  4/day. 


c.  fs  =  i;day. 


Figure  12.  Vari¬ 
ation  of  r\./x\,y\^/ 
t)  and  t|r/T|  unth 


'6 


W/B 


Figure  73.  Variation  ofx\Jr\  with  W/B  and  ff. 


5 


As  channel  boundaries  restrict  the  lateral  dis¬ 
placement  and  accumulation  of  brash  ice,  the 
channel  width  influences  volumes  of  ice  that  can 
be  grown  in  frequently  transited  channels  as  shown 
in  Figure  73,  where  predicted  values  of  rj./rj  are 
related  to  W/B.  For  confined  channels,  whose 
widths  marginally  exceed  vessel  beams,  brash  ice 
cannot  be  displaced  sideways  and  remains  in  the 


track.  Consequently,  as  indicated  in  Figure  73,  less 
ice  grows  in  them  than  in  wider  channels. 

Frequent  transits  of  an  ice-covered  channel 
increase  the  volume  of  ice  growth,  and  it  is  useful 
to  compare  the  increased  ice  volume  with  the 
overall  volume  of  ice  grown  in  a  channel.  Figure  74 
shows  that  the  volume  of  additional  ice  grown 
relative  to  the  overall  ice  volume  diminishes  with 
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increasing  channel  width,  or  (q ,-q)B/(q  W)  de¬ 
creases  as  W/B  increases,  as  is  to  be  expected  from 
simple  algebra.  Frequent  transits  of  a  channel  twice 
the  vessel  beam  results  in  about  a  50%  increase  in 
ice  grown,  for  the  given  calculation  scenario.  The 
increase  reduces  to  about  20%  for  a  channel  five 
times  the  vessel  beam. 

Channel  depth 

The  principal  effect  of  channel  depth  on  ice 
formation  is  essentially  to  restrict  the  displace¬ 
ment  of  brash  ice,  thereby  constraining  its  accumu¬ 
lation  along  a  channel.  A  depth  parameter  is  best 
included  in  the  predictive  model  by  modifying  the 
relationships  for  p  and  e,  the  parameters  describing 
icemovement.  Unfortunately  theice-tankexperiments 
do  not  facilitate  this  modification.  An  appropriate 
depth  parameter  should  be  related  to  the  draft  and 
shape  of  hull,  as  these  parameters  would  act  con¬ 
currently  in  determining  the  trajectories  of  brash- 
ice  displacement.  Observations  of  ice  displace¬ 
ment  around  the  wedge  hull  indicate  that  depth 
effects  would  likely  not  be  significant  because  its 
"keeled"  profile  provided  adequate  space  for  brash 
ice  to  pass  around  and  remain  in  the  track.  On  the 
other  hand,  the  tow  hull,  being  much  blunter  in 
profile,  blocked  the  passage  of  brash  ice  and  has¬ 
tened  ridge  and  plug  formation.  It  was  not  pos¬ 
sible  to  estimate  from  the  few  experiments  with 
the  tow  hull  whether  shallower  channels  transited 
by  tows  or  smaller  values  of  q  ,/q  result  in  less  ice. 
For  keeled  hulls,  such  as  the  wedge  hull  and  deep- 
draft  hulls,  channel  depth  likely  does  not  signifi¬ 
cantly  affect  q./q. 

The  wedge  fiull  has  been  described  as  a  physical 
analogy  of  a  short  tow  with  an  under-hull  accumu¬ 
lation  of  brash  ice,  the  accumulation  being,  in  form 
at  least,  a  sort  of  pseudo-keel.  This  analogy  breaks 
down  for  shallow  channels  (say,  y  JD  <  2),  because 
the  magnitude  of  the  accumulation  diminishes 
and  ice  formation  reverts  to  the  much  less  predict¬ 
able  pattern  depicted  in  Figure  25. 

Conclusion 

Although  circumspection  is  required  when 
appraising  the  predictions  of  numerical  models  of 
ice  formation  in  navigation  channels,  the  present 
predictions  do  usefully  extend  the  information  on 
ice  formation  gained  from  the  ice-tank  experi¬ 
ments.  Admittedly  they  are  an  idealized  simula¬ 
tion  of  ice  formation,  but  they  do  indicate  the 
sensitivities  of  ice  formation  to  air  temperature, 
transit  frequencies  and  channel  width.  Also,  the 
magnitudes  of  ice  growth  and  distributions  of 


brash  ice.  >  concur  reasonably  with  those  formed 
in  the  ice  tank  and  reported  from  field  studies  (e.g. 
Fig.  7, 8). 

As  pointed  out  earlier,  predictions  of  ice  forma¬ 
tion  are  influenced  by  values  assumed  for  P,  e  and 
p,  primarily  because  these  terms  affect  the  removal 
of  ice  from  the  track  or  the  accumulation  of  brash 
ice  within  it.  The  present  predictions  are  based  on 
a  comparatively  simple  predictive  model  sup¬ 
ported  with  data  on  p,  e  and  p  that  were  obtained 
from  ice-tank  experiments  using  a  simplified  hull 
form.  This  hull  is  taken  to  be  representative  of 
hulls  generally  and,  in  particular,  of  tow  hulls  with 
under-hull  accumulations  of  ice.  Different  rela¬ 
tionships  to  those  adopted  here  for  P  and  e  may 
prevail  for  different  hull  forms  and  for  different 
speed  ranges.  For  example,  since  vessels  would 
transit  thinner  ice  covers  at  greater  speeds,  it  is 
possible  that  substantially  different  relationships 
would  result  for  P  and  e  with  number  and  fre¬ 
quency  of  transits  at  such  high  speeds.  Further¬ 
more,  for  actual  channels  the  vagaries  of  transiting 
(e.g.  diverse  hull  sizes  and  shapes  as  well  as  speeds, 
ice  cover  and  weather  conditions)  are  likely  to 
cause  P,  e  and  p  to  vary  randomly  and  unpredicta- 
bly .  Therefore,  it  is  doubtful  that  a  numerical  model 
can  be  sufficiently  reliable  to  predict  precisely  the 
amount  of  ice  produced  in  a  channel  after  a  period 
of  navigation.  This  requirement  of  numerical 
simulation  may  be  moot,  as  the  problems  caused 
by  transiting  are  not  likely  to  be  mitigated  through 
a  predictive  scheduling  program. 

By  way  of  comparison  with  the  predictions  on 
ice  growth  obtained  with  the  present  model,  ice 
growth  predictions  obtained  using  Ashton's  model 
(eq  1 0)  are  briefly  discussed  in  Appendix  D.  Equa¬ 
tion  10  predicts  comparable  values  of  q./q  when /( 
is  on  the  order  of  1/day  or  less.  However,  for 
larger  values  of  f,  eq  10  predicts  much  larger 
values  of  q.  /q,  and  the  values  increase  monotoni- 
cally  with/{. 

ICE  ACCUMULATION  BENEATH 
TOW  BARGES 

The  objective  of  this  part  of  the  ancillary  study 
is  to  determine  the  equilibrium  thicknesses  and 
forms  of  ice  accumulations  beneath  and  ahead  of 
flat-bottomed  tows  moving  through  ice-covered 
channels.  The  scope  of  the  study  is  limited  to  tow 
transits  of  ice  covers  comprising  either  an  ice  sheet, 
a  thin  layer  of  brash  ice,  or  a  layer  of  ice  floes. 
Figure  75  illustrates  these  processes.  The  approach 
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b.  Moving  through  brash  ice. 


Figure  75.  Ice  accumulation  beneath  the  model  barge. 


taken  for  the  study  is  laboratory  experiments 
guided  by  dimensional  analysis. 

Factors  influencing  ice  accumulation 

Ice  accumulation  beneath  tows  is  analogous  to 
the  development  of  sedimentary  waveforms,  par¬ 
ticularly  to  large-scale  ripples  or  aeolian  dunes 
(notably  through  similarity  of  the  flow  field  caus¬ 
ing  wave-form  accumulations). 

The  maximum  equilibrium  thickness  and  the 
form  of  ice  accumulations  beneath  and  ahead  of  a 
moving  barge  are  influenced  by  several  factors, 
including  the  barge  speed  relative  to  the  water;  the 
form,  draft,  beam  and  roughness  of  the  barge  hull; 
the  water  depth;  and  the  dimensions  and  proper¬ 
ties  (e.g.  strength,  roughness,  porosity)  of  the  ice. 
Other  factors,  including  the  temperatures  of  the 
water,  air  and  ice,  may  affect  accumulation  thick¬ 
ness.  If  the  development  and  equilibrium  forms  of 


ice  accumulations  are  treated  as  wave-like  forma¬ 
tions  of  particles,  then  the  following  functional 
relationship  should  be  considered: 


where  h  =  height  of  the  wave 

X  =  wavelength  of  the  wave 
u*  =  shear  velocity  of  the  water  flow 
u*c  =  its  critical  value  associated  with  parti¬ 
cle  entrainment 

(  =  representative  length  dimension  of  ice 
pieces 

v  =  kinematic  viscosity  of  water 
yo  =  water  depth. 

Wave  forms  in  particulate  materials  are  appropri¬ 
ately  described  in  terms  of  their  wavelength  X. 
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Wave  height  h  (or  steepness  h/X)  is  more  complex 
than  X  because  it  is  determined  by  a  larger  number 
of  variables.  However,  for  ice  accumulation  be¬ 
neath  a  barge,  complications  arise  when  attempt¬ 
ing  to  predict  equilibrium  accumulation  forms 
(i.e.  X  and  h).  The  rate  of  ice  discharge  beneath  the 
barge  is  influenced  by  the  barge  draft  and  geome¬ 
try  (particularly  the  bow);  the  ice  characteristics; 
the  relatively  small  area  of  accumulation  bound¬ 
ary  (bottom  of  barges);  the  relatively  large  size  of 
ice  fragments  compared  to  the  barge  size;  and,  for 
rivers  and  many  harbors,  the  shallowness  of  the 
water.  These  factors  prevent  the  use  of  a  purely 
analytical  approach  to  determining  equilibrium 
thickness  of  ice  accumulations  beneath  a  barge. 
Therefore,  it  is  appropriate  to  go  directly  to  dimen¬ 
sional  analysis  to  arrive  at  a  set  of  parameters  for 
describing  the  accumulation  thickness. 

Buckingham's  Pi  theorem  yields 

f  [h_  X  V  *0  'B  ' 

V'y0-D'  V^i-p./p)Ti '  ^  '  y0-D'  e ' 


L  l  IL 
BVo-D't 


,  0,  p^  p;,  ice  properties)  =  0 


(35) 


need  not  be  explicitly  included  in  eq  35,  as  the  flow 
is  fully  turbulent.  For  barge  operators  who  wish  to 
know  the  thickness  of  an  accumulation  beneath  a 
barge,  the  following  functional  relationship  is  of 
interest: 


B  L  tj_ 
t '  B  '  l ' 


(36) 


The  wave-steepness  parameter  h/X  is  replaced  by 
h  /t|,  as  it  can  be  argued  for  ripple-like  waveforms, 
or  for  waveforms  in  infinitely  deep  water,  that  X  is 
proportional  to  q.  Also,  for  the  present  study,  X 
could  not  be  accurately  defined,  as  most  experi¬ 
ments  resulted  in  single  wave-like  accumulations 
that  commenced  beneath  the  rounded  bow  of  a 
model-scale  test  barge.  In  eq  35  the  parameter  y  / 
(yo-D)  is  a  blockage  term  and  is  included  in  eq  36 
as  F*l=  F^y  /(yQ-D)]  to  account  for  the  effects  on 
h  of  water  depth  relative  to  barge  draft;  l / ( yo~D) 
is  a  measure  of  waveform  roughness  relative  to 
water  depth;  and  B/i  relates  barge  beam  to  rubble 
size  and  thereby  to  potential  wave  size. 

A  formulation  exercise  similar  to  the  one  that 
led  to  eq  35  can  be  followed  to  identify  a  paramet¬ 
ric  relationship  for  accumulation  slope  0a: 


where  V  =  barge  speed  relative  to  the  water 
p.  =  density  of  the  ice  rubble 
p  =  density  of  water 
0  =  stem  angle  of  the  barge  bow 
p^  =  coefficient  of  static  friction  between 
the  ice  and  the  hull 

p.  =  coefficient  of  static  friction  between 
ice  and  ice 
B  =  barge  beam 
D  =  barge  draft 
L  ~  barge  length 

q  =  thickness  of  the  ice  fragments 
t  =  major  length  of  the  ice  fragments. 

Ice  properties  refer  to  the  strength,  temperature, 
shape,  etc.  of  ice  fragments.  Strictly  speaking,  for 
describing  wave  topography  on  extensive  planes, 
the  parameter 
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Pj  ,  0  and  ice  properties)  j  . 


Before  discussing  the  results  of  the  experiments,  it 
is  helpful  to  consider  briefly  the  influence  of  ice- 
rubble  proportions,  or  the  term  r\/f,  on  ice-rubble 
accumulations  beneath  a  moving  barge.  As  the 
critical  value  of  Fr  associated  with  the  entrainment 
of  ice-rubble  pieces  lodged  against  the  flat  bottom 
of  a  barge  Pfc  increases  with  decreasing  q //,  it  is 
likely  that  larger  values  of  q  /t  would  lead  to 
thicker  accumulations.  Tatinclaux  and  Gogus 
(1981),  for  example,  proposed  the  following  equa¬ 
tion  for  assessing  the  stability  of  ice  blocks  lodged 
beneath  a  flat  ice  cover: 


Vg(l-p/p)n  =Fr 

should  be  replaced  by  u./qV  However,  as  u.c  is 
not  readily  determined  for  ice  rubble  beneath 
barges,  it  is  more  convenient  to  use  F  .  Addition¬ 
ally,  the  second  term  in  eq  34,  a  Reynolds  number. 


V[(l-pj/p)CT]"a5  = 

[c,  M)2+C2  M)+c3]- 


where  V  is  the  mean  velocity  of  the  water  flow 
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(below  the  cover)  associated  with  the  entrainment 
of  blocks,  and  Cy  C2  and  C3  are  empirical  coeffi¬ 
cients. 

Experiments 

A  series  of  experiments  was  conducted  to  deter¬ 
mine  the  equilibrium  thicknesses  and  forms  of  ice 
accumulations  around  barges.  In  particular,  the 
relationship  between  h/r\  and  F*rwas  sought  for 
varying  forms  of  ice  cover.  The  experiments  en¬ 
tailed  the  use  of  one  of  the  l:30-scale  barges  de¬ 
scribed  earlier.  The  model  barge  was  1.8  m  long, 
0.36  m  wide  and  of  variable  draft,  with  0  =  21°.  It 
was  towed  at  speeds  of 0.03-0.45  m/s.  The  experi¬ 
mental  set-up  is  shown  in  Figure  76.  Two  sets  of 


experiments  were  conducted:  one  using  a  glass- 
sided  flume;  the  other,  ffi-IR's  ice  tank.  The  flume 
was  20  m  long,  0.3  m  deep  and  0.76  m  wide.  The  ice 
tank  is  depicted  in  Figure  17. 

For  the  flume  experiments,  ice  covers  were 
simulated  using  four  sizes  of  polyethylene  blocks 
or  beads.  In  the  ice  tank,  sheets  of  urea  ice  (4.5  and 
10  mm  thick)  were  used  to  simulate  ice  sheets  and 
ice  rubble.  The  dimensions  of  the  blocks  are  given 
in  Table  3.  Polyethylene  has  the  same  specific 
gravity  as  ice  (0.92).  The  coefficient  of  static  friction 
for  submerged  polyethylene  against  the  painted 
finish  of  the  model  barge  was  measured  to  be  0.15; 
for  urea  ice,  j Jt=  0.05.  These  values  were  obtained 
by  pulling  blocks  over  the  surface  of  the  model 
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Figure  76.  Experimental  set-up. 
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Table  3.  Properties  of  simulated  ice. 


Type 

Form 

Length 

(mm) 

Width 

(mm) 

Thickness 

(mm) 

B  H* 

A 

beads 

4.8 

4.8 

3 

74.1 

B 

large  block 

38 

32 

10 

9.4 

C 

medium  block 

38 

32 

3 

9.4 

D 

small  block 

13 

13 

3 

27.4 

E 

urea  ice 

200t 

50+ 

10 

1.8 

F 

urea  ice  sheet 

200t 

45+ 

4.5 

1.8 

G 

urea  ice  blocks 

50t 

50+ 

5 

7.1 

H 

urea  ice  blocks 

45t 

45+ 

4.5 

7.9 

*  B  =  356  mm 
t  Average 


barge.  We  realize  the  limitations  of  using  plastic  to 
simulate  ice:  somewhat  higher  values  of  friction 
coefficient,  the  lack  of  freezing  between  blocks,  a 
non- wetting  surface,  etc.  Tatinclaux  and  Lee  (1978) 
discussed  at  length  the  limitations  of  using  plastic 
to  simulate  ice.  However,  in  addition  to  using  urea 
ice,  polyethylene  blocks  constitute  an  appropriate 
model  ice  material  for  this  study  because  the 
strength  properties  of  ice  are  not  considered,  and 
the  accumulation  processes  mainly  involve  the 
movement  of  fully  submerged  blocks. 

The  38-  x  32-  x  1 0-mm  polyethylene  blocks  were 
arrayed  in  a  tight,  closely  packed  formation  of 
100%  areal  concentration  such  that  they  simulated 
an  ice  sheet  precut  into  the  sizes  of  broken  ice  likely 
to  result  when  a  barge  moves  through  it.  At  a  scale 
of  1:30  the  array  of  plastic  blocks  simulated  a  0.3- 
m-thick  ice  sheet  breaking  into  ice  blocks  with  a 
thickness-to-length  ratio  of  about  1 :3-4.  The  beads 
were  used  to  simulate  a  cover  of  relatively  small 
(0.09  m  diameter,  prototype)  ice  rubble.  The  areal 
concentration  of  the  beads  was  typically  about 
75-78%.  The  two  smaller  sizes  of  ice  block  were 
used  tosimulatea  cover  of  fairly  small  (0.4  to  1.0m 
square  x  0.09  m  thick,  prototype  dimensions)  ice 
rubble  arrayed  with  an  average  areal  concentra¬ 
tion  of  about  66%.  In  the  course  of  the  experiment 
it  was  found  that  the  areal  concentration  of  blocks 
or  beads  affected  the  length  of  barge  transit  re¬ 
quired  for  an  equilibrium  accumulation  to  de¬ 
velop  beneath  the  model  barge;  lower  concentra¬ 
tions  of  blocks  required  longer  transits.  This  was 
one  source  of  scatter  in  the  data.  Another,  less 
significant,  source  was  buoyancy,  which,  for  thick 
accumulations  of  ice  rubble,  slightly  altered  the 
mean  draft  and  trim  of  the  barge. 

For  each  experiment  the  barge  was  towed  at 
constant  speed  through  the  simulated  ice  covers. 
Once  ice  accumulation  below  the  barges  had  at¬ 


tained  an  equilibrium  thickness,  which  usually 
occurred  after  about  two  to  three  hull  lengths  of 
transit,  the  accumulation  profile  was  measured  by 
means  of  a  small-gauge  hook  probe.  As  can  be 
appreciated  from  the  photographs  accompanying 
this  section,  the  accumulation  thickness  could  be 
measured  to  an  accuracy  of  about  one  block  thick¬ 
ness  for  large  accumulations.  Measured  values  of 
h  and  estimates  of  X  are  tabulated  in  Table  4. 

A  brief  set  of  additional  experiments  was  per¬ 
formed  with  two  barges  moving  in  parallel  along 
the  flume.  It  was  found  that  the  narrowness  of  the 
flume — 7%  wider  than  the  beam  of  the  double 
barge  configuration — occasionally  caused  the 
barges  to  become  stuck  when  the  simulated  ice 
blocks  wedged  between  the  barges  and  the  side 
walls  of  the  flume.  This  is  also  a  problem  experi¬ 
enced  by  barges  moving  through  locks  or  confined 
channels. 

Observations 

The  manner  in  which  ice  rubble  accumulated 
beneath  and  ahead  of  the  model  barge  differed 
according  to  whether  or  not  it  transited  a  simu¬ 
lated  ice  sheet  or  a  simulated  layer  of  broken  (floes 
or  brash)  ice.  For  both  cases,  ice-rubble  accumula¬ 
tion  beneath  the  barge  was  analogous  to  alluvium 
wave  formation.  Figures  77  and  78  illustrate  accu¬ 
mulations  under  the  model  barge  moving  through 
simulated  ice  covers  of  sheet  ice  and  an  array  of 
floes,  respectively.  The  development  of  a  false 
bow  and  a  waveform  accumulation  of  simulated 
brash  ice  (rubble  with  full-scale  dimensions  of  less 
than  about  2  m)  are  depicted  in  Figures  79  and  80, 
respectively.  Figure  81  shows  the  especially  troub¬ 
ling  condition  of  an  accumulation  that  extends  to 
the  channel  bottom.  This  condition  occurred  when 
the  test  barge  was  towed  at  creeping  speed  close  to 
the  flume  wall. 
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Table  4.  List  of  experiments:  ice  accumulation  beneath  and  ahead  of 
the  l:30-scale  model  barge. 


Test 

no.* 

Ice 

type'* 

V 

(mis) 

y0 

(m) 

D 

(m) 

h 

(mm) 

A 

(m) 

i 

A 

0.109 

0345 

0.06 

33 

1.60 

?. 

0.163 

30.0 

0.39 

3 

0.190 

403 

0.43 

4 

0.210 

54.4 

055 

5 

0.244 

47.7 

0.80 

6 

0311 

34.6 

1.10 

7 

0.446 

7.4 

1.60 

8 

0.109 

0309 

0.06 

33 

1.60 

9 

0.163 

28.8 

037 

10 

0.176 

333 

0.49 

11 

0.190 

35.6 

054 

12 

0.210 

41.6 

0.60 

13 

0344 

35.5 

0.86 

14 

0311 

23.0 

1.60 

15 

0.446 

65 

1.60 

16 

0.109 

0.150 

0.06 

33 

1.60 

17 

0.176 

24.2 

058 

18 

0.190 

35.2 

0.74 

19 

0.210 

33.1 

0.90 

20 

0344 

26.7 

138 

21 

0311 

12.0 

1.60 

22 

0.446 

3.2 

1.60 

23 

0.163 

0.143 

0.06 

12.6 

0.72 

24 

0.176 

23.6 

0.78 

25 

0.190 

32.1 

0.84 

26 

0310 

29.1 

0.86 

27 

0344 

20.2 

1.37 

28 

0.163 

0.117 

0.06 

83 

0.71 

29 

0.176 

203 

0.64 

30 

0.190 

19.0 

1.08 

31 

0.210 

12.3 

134 

32 

0.190 

0.081 

0.06 

33 

1.60 

33 

0.176 

33 

1.60 

34 

0344 

33 

1.60 

35 

0.311 

33 

1.60 

36 

0.446 

33 

1.60 

37 

B 

0.041 

0.235 

0.088 

51.7 

055 

38 

0.057 

41.3 

0.44 

39 

0.086 

355 

0.62 

40 

0.143 

29.0 

0.62 

41 

0301 

9.6 

1.11 

42 

0.041 

032 

0.088 

59.9 

055 

43 

0.088 

325 

055 

44 

0.142 

24.2 

0.69 

45 

0.028 

0303 

0.088 

41.3 

0.43 

46 

0.086 

37.0 

050 

47 

0.143 

19.1 

0.75 

48 

0301 

9.6 

1.60 

83 


Table  4.  Continued. 


Test 

Ice 

V 

y0 

D 

h 

A 

no. 

type 

( m/s ) 

(m) 

(m) 

(mm) 

(m) 

49 

0.041 

0.178 

0.088 

49.0 

0.44 

50 

0.088 

32.5 

0.66 

51 

0.143 

19.1 

0.86 

52 

0201 

9.6 

155 

53 

0.041 

0.168 

0.088 

37.8 

0.56 

54 

0.088 

29.4 

0.62 

55 

0.142 

18.8 

0.86 

56 

0.041 

0.231 

0.06 

65.3 

0.69 

57 

0.088 

56.9 

0.75 

58 

0.142 

40.5 

0.83 

59 

0.203 

18.5 

0.95 

60 

0.028 

0.203 

0.06 

62.0 

052 

61 

0.143 

28.6 

0.74 

62 

0.201 

19.1 

0.95 

63 

0.041 

0.178 

0.06 

53.2 

0.45 

64 

0.086 

41.3 

0.65 

65 

0.143 

28.6 

0.70 

66 

0201 

9.6 

0.86 

67 

c 

0.041 

0.180 

0.06 

30.0 

027 

68 

0.075 

48.0 

050 

69 

0.109 

66.3 

057 

70 

0244 

20.8 

1.60 

71 

0.446 

13.4 

1.60 

72 

D 

0.142 

0240 

0.088 

21.2 

034 

73 

0.159 

37.7 

0.48 

74 

0.176 

34.0 

058 

75 

0244 

27.6 

0.86 

76 

0379 

19.4 

1.60 

77 

0.446 

12.0 

1.60 

78 

E 

0.042 

1.067 

0.088 

30.5 

79 

0.144 

51.9 

80 

0.163 

54.9 

81 

0.202 

47.3 

82 

0245 

30.5 

83 

0.447 

0.0 

84 

F 

0.042 

1.067 

0.088 

51.9 

85 

0.144 

12.2 

86 

0.221 

12.2 

87 

G 

0.144 

1.067 

0.088 

13.7 

88 

0.202 

41.2 

89 

0221 

27.5 

90 

0.245 

15.2 

91 

0.447 

10.0 

92 

H 

0.144 

1.067 

0.088 

33.6 

93 

0202 

40.3 

94 

0221 

70.2 

95 

0245 

48.8 

% 

0.447 

45 

*  Total  number  of  experiments:  258. 
**  Types  are  illustrated  in  Table  3. 
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Figure  77.  Simulated 
ice  accumulated  be¬ 
neath  the  model  barge 
using  medium-size 
blocks  arrayed  as  an 
ice  sheet. 


Figure  78.  Simulated 
ice  accumulated  be¬ 
neath  the  model  barge 
using  large  blocks. 


Fig  lire  79.  Falsc-bozc 
format  ion  ahead  of  the 
model  barge. 
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Figure  80.  Dune-like  accumulation  of  ice  beneath  the  model  barge  using  beads. 


Figure  81 .  Ice  accumulation  extending  to  the  channel  bed  when  the  barge  crept  along  the  flume  wall. 
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When  themodel  bargemoved  through  thelayer 
of  blocks  arrayed  to  simulate  an  ice  sheet  breaking 
into  a  more-or-less  regular  train  of  blocks,  it  over¬ 
rode  the  blocks  and  pushed  them  beneath  its  bow. 
On  rounding  the  bow  and  coming  to  the  flat  bot¬ 
tom  of  the  barge,  the  initial  train  of  blocks  was 
dislocated,  causing  blocks  to  pile  up  and  thereby 
initiating  an  ice  accumulation.  A  relatively  quies¬ 
cent  region  of  flow,  caused  by  flow  separation 
immediately  downstream  of  the  bow,  aided  the 
accumulation  process  such  that,  for  relatively  low 
values  of  F,  blocks  could  accumulate  as  a  more-or- 
less  randomly  oriented  jumble  beneath  the  bow. 
The  accumulation  thickened  until  fluid  forces 
strengthened  enough  to  balance  the  buoyancy 
forces  acting  through  the  blocks  resting  along  the 
accumulation  perimeter.  With  increasing  hull 
speed,  fluid  forces  acting  on  the  blocks  increased 
in  magnitude  such  that  the  maximum  thickness  of 
accumulation  diminished,  and  then  the  accumula¬ 
tion  spread  over  a  greater  length  of  the  barge.  Also, 
the  maximum  thickness  of  the  accumulation 
diminished  with  decreasing  depth  of  water  be¬ 
neath  the  model  barge.  Imbrication  (overlapping 
or  shingling)  of  plate  ice  fragments  may  increase 
the  resistance  to  erosion  of  an  accumulation,  as  can 
be  seen  from  Figure  78. 

Small  or  irregularly  shaped  pieces  of  ice  dislo¬ 
cated  from  one  another  and  did  not  move  together 
when  shoved  ahead  of  the  test  barge.  Instead, 
when  the  model  barge  moved  through  the  layer  of 
beads  simulating  a  cover  of  brash  ice  or  pans  of 
frazil  ice,  a  false  bow  developed  ahead  of  the  barge 
withoutsignificantaccumulationforming  beneath 
it  (Fig.  79).  The  false  bow  streamlined  the  bow, 
deflecting  ice  rubble  around  the  barge  and  allow¬ 
ing  only  a  trickle  of  ice  to  pass  beneath  it.  With 
increasing  speed  the  false  bow  shortened  and  thick¬ 
ened  such  that  the  beads  could  resist  the  fluid 
forces,  and  a  greater  discharge  of  beads  passed 
beneath  the  hull,  causing  a  dune-like  accumula¬ 
tion  to  develop,  as  illustrated  in  Figure  80.  For  high 
speeds,  possibly  not  of  practical  importance  for 
conventional  barge  traffic,  a  false  bow  did  not 
develop,  and  the  accumulation  thickness  decreased 
until  only  a  single  layer  of  beads  flowed  beneath 
the  barge.  For  the  range  of  speeds  associated  with 
declining  thickness  of  accumulation,  it  was  not 
uncommon  for  rubble  to  accumulate  as  two  or 
more  waves  beneath  the  test  barge,  though  in  most 
cases  only  a  single  wave  formed  under  the  single 
model  barge.  It  may,  therefore,  be  possible  that 
many  waveform  accumulations  may  develop 
beneath  the  extensive  bottom  area  of  full  barges  or 


a  tow  comprising  several  barges. 

The  smaller  blocks  that  simulated  a  cover  of 
randomly  arrayed  ice  pieces  behaved  similarly  to 
the  cover  of  beads.  Lacking  the  facility  to  readily 
interlock  when  shoved,  the  smaller  blocks  accu¬ 
mulated  as  a  conical  false  bow  that  deflected  blocks 
around  the  barge.  With  increasing  speed,  the  false 
bow  shortened  but  thickened  such  that  more  blocks 
could  pass  beneath  the  barge  and  collect  in  a 
waveform  accumulation. 

Equilibrium  thickness  of  the  ice  accumulation 

The  relationships  between  h/r\  and  Fr  differ 
according  to  whether  or  not  the  barge  transited 
simulated  ice  sheets  or  layers  of  ice  rubble.  In  each 
case,  for  polyethylene  blocks,  the  maximum  value 
of  h/x\  was  about  17-21  for  B/g  >  9.  Barge  transit 
through  ice  sheets  is  shown  in  Figure  82.  Figures 
83  and  84  show  the  variation  of  h/T)  with  Fr through 
broken  ice  covers.  In  all  three  figures,  g/(yo~D)  and 
B / 1  are  indicated. 

If  the  barge  moved  at  creeping  speed,  the  ex¬ 
perimental  set-up  was  truly  two  dimensional,  and 
the  length  of  simulated  ice  sheet  transited  was 
suitably  long,  ice  rubble  would  tend  to  accumulate 
to  the  full  water  depth  beneath  the  model  barge.  In 
other  words,  the  maximum  value  of  h/r\  should 
have  approached  (yo-D)/q  as  Fr  approached  zero, 
provided  the  barge  or  tow  was  large  enough  to 
accommodate  the  accumulation.  This  occurred 
when  the  test  barge  was  towed  close  to  the  flume 
wall  or  when  a  double  barge  configuration  was 
used,  because  the  accumulation  was  confined  by 
the  flume's  side  walls.  However,  for  a  single  barge 
moving  slowly  through  the  flume,  whose  width 
was  2.1  times  the  barge's  beam,  the  accumulation 
became  less  two  dimensional  as  blocks  could  slide 
or  be  swept  from  the  heaped  sides  of  the  accumu¬ 
lations.  Consequently,  for  a  barge  moving  through 
a  relatively  wide  channel,  the  values  of  h/r\  are 
likely  to  be  less  than  (yo~D)/r),  and,  for  the  same 
reason,  they  are  likely  to  diminish  somewhat  with 
decreasing  value  of  B/g. 

Figure  82  shows  that  h/r\  diminishes  with  in¬ 
creasing  valuew  of  F^  until  Ff  attains  some  critical 
value,  in  excess  of  about  8,  beyond  which  all  blocks 
swept  beneath  the  barge  move  as  a  single  layer  ( h / 
T|  =  1).  Although  thedata  collapse  reasonably  well, 
they  vary  consistently  with  t/(yo~D),  reflecting 
perhaps  the  influence  of  accumulation  roughness 
on  local  flow  velocity  and  fluid  forces. 

For  barge  transit  through  broken  or  brash  ice, 
the  development  of  a  false  bow  complicates  the 
relationship  between  h/r\  and  Frsuch  that  h/r\ 


87 


attains  a  maximum  value  for  a  nonzero  value  of  F^ 
as  shown  in  Figures  83  and  84.  The  relationship  is 
analogous  to  that  between  bed-form  height  and 
flow  velocity  for  flow  over  alluvial  beds  (e.g. 
Raudkivi  1976).  Bow  and  false-bow  forms  affect 
the  discharge  of  ice  beneath  a  hull.  As  the  flow  of 
ice  rubble  beneath  the  bow  increased,  so  did  the 
height  of  the  accumulation  until,  with  further 
increases  in  both  ice  and  water  flow,  the  accumu¬ 
lation  became  washed  out  and  rubble  moved  as  a 
single  layer  beneath  the  barge.  In  Figure  83  the 
larger  peak  value  of  h  /r\  that  was  measured  for  the 
larger  B/t  can  be  ascribed  to  imbrication,  or  layer¬ 
ing,  of  blocks;  imbrication,  which  did  not  occur 
extensively  for  the  smaller  blocks  (Fig.  77,  78) 
increases  the  resistance  of  the  block  to  entrain¬ 
ment.  Additionally,  in  accordance  with  eq  38, 
longer  blocks  lodged  beneath  a  planar  boundary 
are  more  resistant  to  entrainment  and  therefore 


may  enable  larger  accumulations  to  develop. 

Figure  84  should  be  considered  together  with 
Figure  85,  which  shows  the  variation  of  false-bow 
extent,  or  angle  0a,  with  Ff.  As  the  false  bow  dimin¬ 
ished  in  size  with  increasing  barge  speed,  more  ice 
was  transported  beneath  the  barge  so  as  to  form  an 
accumulation  wave.  The  geometry  of  the  false 
bow  formed  at  the  bow  of  the  model  barge  can  be 
described  in  terms  of  the  angle  8a,  as  depicted  in 
Figure  85.  With  increasing  hull  speed  F ^  0fl  in¬ 
creased  until  it  was  about  the  same  value  as  the 
bow's  stem  angle  8.  When  0a exceeded  8,  the  false 
bow  no  longer  developed  and  ice  was  swept  be¬ 
neath  the  barge.  As  described  in  the  foregoing 
section,  the  peak  accumulation  beneath  the  model 
barge  moving  through  brash  ice  occurred  when  8a 
approximately  equaled  0. 

An  interesting  result  not  readily  apparent  in 
Figure  84  is  that,  if  the  maximum  thickness  of 
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Figure  83.  Accumulation  thickness  beneath  the  model  barge  for  the  simulated  cover  of  ice  floes. 


accumulation  h  is  normalized  as  h/(yo~D),  the 
peak  values  of  h/ (yo~D)  range  from  0.26  to  0.36, 
with  an  average  of  0.31 .  This  result  is  remarkably 
similar  to  the  maximum  values  of  dune  amplitude 
normalized  with  water  depth  for  alluvial  river 
channels  (Raudkivi  1976). 

A  typical  barge  (9  m  wide  with  a  2.6-m  draft) 
moving  at  a  speed  of  0.5  m/s  through  a  0.2-m-thick 
ice  sheet  on  a  6-m-deep  pool  of  the  Mississippi 
River  (F  =  2.2,  £/(yo~D)  =  0.25)  could  possibly 
accumulate  about  3  m  of  ice  rubble  beneath  its 
bottom  (Fig.  81).  If  the  ice  cover  comprised  brash 
ice  with  piece  0.1  m  thick,  and  the  barge  were 
moving  at  speed  of  about  1.0  m/s  (F  =  6.1),  an 
accumulation  about  1.5  m  thick  may  develop. 
These  examples  indicate  the  potential  magnitudes 
of  ice  accumulations.  Obviously  many  factors  af¬ 
fect  accumulation  and  its  precise  estimation.  For 
example,  if  the  barge  were  moving  slowly  through 
a  confined  channel,  such  as  a  lock  approach,  the 


accumulations  may  become  much  thicker.  These 
estimates  do  not  take  into  account  the  effect  of  ice 
fragments  freezing  to  one  another  or  the  lower 
values  of  friction  coefficient  that  occur  between  ice 
than  occur  between  plastic.  Freezing  would  in¬ 
crease  the  resistance  of  an  accumulation  to  erosion 
by  fluid  forces  and  would  likely  result  in  larger 
accumulations.  Lower  values  of  friction  coeffi¬ 
cient  may  have  the  opposite  influence. 

Discussion 

Several  field  observations  cited  earlier  attest 
that  significant  quantities  of  ice  do  accumulate 
beneath  flat-bottomed  tows  moving  through  ice- 
covered  channels.  The  present  study  shows  that 
the  manner  of  ice  accumulation  varies  according 
to  whether  an  ice  cover :  an  ice  sheet  or  a  layer  of 
ice.  If  a  cover  is  an  ice  sheet,  broken  ice  can  be 
directly  submerged  beneath  a  transiting  tow,  and 
the  maximum  thickness  of  the  accumulation  de- 
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creases  with  increasing  tow  speed,  or  F  .  If  an  ice 
cover  comprises  comparatively  small  ice  fragments 
(e.g.  brash  ice),  a  false  bow  may  form  ahead  of  a 
slow-moving  tow,  deflecting  ice  around  the  tow. 
With  increasing  F,  the  false-bow  size  diminishes 
such  that  the  flow  of  ice  beneath  the  tow  increases 
and  waveform  accumulations  develop.  The  accu¬ 
mulation  thickness  peaks  with  increasing  F.  and 
decreases  thereafter. 

The  thickness  of  an  ice  accumulation  beneath  a 
tow,  expressed  as  h/r\,  may  attain  maximum  val¬ 
ues  of  about  17-21.  However,  accumulations  ex¬ 
tending  to  the  channel  bed  may  form  beneath 
large,  slow-moving  tows  in  narrow  or  confined 
channels. 

The  use  of  q  for  normalizing  h  is  not  altogether 
satisfactory,  but  neither  is  any  other  dimension, 
such  as  yo-D  or  £.  A  drawback  in  using  q  is  that 
measurement  accuracy  is  on  the  order  of  tv  The 
use  of  yo-D  seems  attractive,  especially  for  tow 
transits  of  loose  covers  of  brash  ice  over  shallow 
channels,  for  which  h/(yo~D)  attains  a  maximum 
value  of  about  0.31.  However,  for  deep  channels 
(yo»D ),  the  water  depth  does  not  directly  affect 
the  local  flow  field  around  a  hull,  and  therefore  the 
use  of  yo-D  is  not  appropriate.  The  quandary  here 
is  the  same  as  that  faced  in  the  appropriate  de¬ 
scription  of  sediment  waveforms.  Whereas  they 
are  often  conveniently  described  and  calculated  in 
terms  of  waveform  steepness  (amplitude  divided 
with  wavelength),  tow  skippers  are  earnestly  in¬ 
terested  in  knowing  the  maximum  accumulation 
thickness  relative  to  some  familiar  dimension  such 
as  ice  thickness,  channel  depth  or  hull  draft.  They 
are  not  interested  in  knowing  it  relative  to  any¬ 
thing  as  nebulous  as  accumulation  wavelength. 


IMPROVED  TRANSITING  OF  ICE-COVERED 
CHANNELS 

Frequent  transiting  of  navigation  channels 
under  frigid  conditions  incurs  two  fundamental 
problems.  One  problem  is  simply  that  frequent 
transiting  causes  more  ice  to  grow  than  would 
grow  without  transiting.  The  other  problem  arises 
because  transiting  causes  most  of  the  ice  grown 
along  tracks  opened  by  vessels  to  be  broken  and 
metamorphosed  to  brash  ice,  which  may  collect  as 
thick  accumulations  that  eventually  impede  tran¬ 
siting.  The  ice-tank  experiments  indicate  that,  for 
tow  vessels,  the  latter  problem  is  generally  more 
severe,  as  transiting  leads  to  the  formation  of  thick 
localized  accumulations  that  pose  major  impedi¬ 


ments  to  sustained  transiting. 

A  practical  objective  of  the  present  study  was  to 
identify  approaches  that  could  improve  transiting 
of  ice-covered  channels  and  to  avoid  or  mitigate 
the  problems  caused  by  frequent  transiting.  The 
dual  nature  of  these  problems  suggests  two  gen¬ 
eral  approaches  for  ice  control  and  problem  miti¬ 
gation.  One  approach  entails  minimizing  ice 
growth  and  requires  that  the  number  and  sched¬ 
ule  of  transits  be  optimized  in  accordance  with 
prescribed  economic  constraints.  This  approach 
gives  rise  to  the  concept  of  optimum  scheduling  of 
transits  and  convoying  of  vessels.  The  other  ap¬ 
proach  entails  the  use  of  mechanical  methods  for 
controlling  brash-ice  accumulations  at  specific 
channel  locations  where  especially  troublesome 
accumulations  perennially  form.  Before  outlining 
approaches  to  problem  control  or  mitigation,  it  is 
useful  to  summarize  the  troublesome  consequences 
that  transiting  may  incur. 

Troublesome  consequences  of  transiting 

The  results  obtained  from  the  ice-tank  experi¬ 
ments  and  the  numerical  model,  together  with 
observations  cited  from  the  field  studies  (reviewed 
earlier),  show  that  tracks  opened  by  transiting 
vessels  become  covered  with  rather  porous  layers 
of  brash  ice  that  are  approximately  1 .5-3  times  the 
ice-cover  thickness.  The  higher  thickness  value  is 
associated  with  higher  frequencies  of  transit.  The 
equivalent  solid  volume  of  brash  ice  grown  at  a 
typical  channel  section  is  about  1. 5-2.5  times  the 
volume  of  ice  cover  that  would  have  formed  over 
vessel  tracks.  Both  the  experiments  and  the  nu¬ 
merical  model  indicate  that  frequent  transiting 
does  not  lead  to  as  dramatic  an  increase  in  the 
volume  of  ice  grown  as  some  of  the  literature 
would  suggest.  The  experiments  do  show  that 
higher  frequencies  result  in  finer  sizes  of  brash  ice 
and  less  porous  accumulations.  In  addition  to 
increasing  the  resistance  to  transit,  increased  vol¬ 
umes  of  ice  growth  in  river  channels  may  exacer¬ 
bate  problems  associated  with  ice-jam  formation 
or  ice  passage  at  dams.  The  results  of  the  present 
study  show  that,  although  transiting  certainly  leads 
to  additional  ice  growth,  the  increase  in  the  vol¬ 
ume  of  ice  grown  on  most  river  channels  is  not 
great.  For  example,  transiting  a  channel  that  is  five 
times  the  width  of  vessel  beam  produces  about 
20%  more  ice  than  would  form  without  transiting 
(Fig.  73).  On  the  other  hand,  transiting  a  narrow 
channel  that  is  only  twice  the  width  of  transiting 
vessels  may  double  the  amount  of  ice  grown.  For 
tow  transiting  of  channels  through  shallow  rivers. 
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possibly  of  greater  concern  than  the  volume  of  ice 
grown  is  the  hazard  that  plugs  of  accumulated 
brash  ice  may  precipitate  ice-jam  formation  by 
being  more  resistant  to  break-up  than  are  ice  cov¬ 
ers. 

Although  brash-ice  layers  in  vessel  tracks  in¬ 
crease  the  resistance  to  transiting,  ridges  of  brash 
ice  formed  beneath  ice  covers  flanking  the  tracks 
pose  greater  impediments  to  transiting.  For  river 
channels  transited  by  tows,  plugs  of  brash  ice  that 
accumulate  locally  in  tracks  are  an  additional  and 
formidable  impediment.  Ridges  and  plugs  maybe 
considerably  thicker  than  brash-ice  layers  in  ves¬ 
sel  tracks.  The  ice-tank  experiments  show  that 
ridges  can  be  about  three  to  five  times  the  ice-cover 
thickness  and  that  the  plugs  can  be  even  thicker  at 
about  ten  or  more  times  the  ice-cover  thickness. 
The  ice-tank  experiments  also  show  that  ridges 
and  plugs  form  more  rapidly  in  shallower  chan¬ 
nels  and  may  extend  the  entire  channel  depth.  It  is 
to  be  expected,  therefore,  that  ridge  and  plug 
formation  can  create  especially  troublesome  con¬ 
ditions  for  transiting  ice-covered  shallow  chan¬ 
nels.  Additionally,  brash-ice  accumulations  are 
likely  to  create  severe  difficulties  for  the  winter  use 
of  small  confined  harbors,  even  if  they  are  deep. 

A  further  accumulation-related  problem  occur¬ 
ring  for  tows  is  their  propensity  to  entrap  and 
transport  brash  ice  beneath  their  flat-bottomed 
hulls.  Such  accumulations  not  only  increase  the 
tow  resistance,  but  they  may  cause  tows  to  ground 
when  transiting  shallow  channels  or  attempting  to 
pass  through  locks. 

Approaches  to  problem  control  and  mitigation 

Two  general  approaches  can  be  implemented: 
minimization  of  ice  formation  through  optimal 
scheduling  of  transits,  or  mechanical  redistribu¬ 
tion  (and  possibly  removal)  of  troublesome  accu¬ 
mulations  of  brash  ice.  The  two  approaches  could, 
of  course,  be  implemented  jointly.  However,  the 
present  study  indicates  that  the  latter  approach  is 
probably  more  practicable  and  promising. 

Optimal  schedules  for  transiting  of  ice-covered 
channels  can  ideally  be  identified  by  mathemati¬ 
cal  optimization  analyses.  Such  analyses  are  be¬ 
yond  the  scope  of  the  present  study  as,  in  addition 
to  the  constraint  of  minimizing  volume  of  ice 
growth,  constraints  related  to  navigation  econom¬ 
ics  need  to  be  taken  into  account;  these  constraints 
include  minimizing  the  power  requirements  of 
transiting  vessels,  minimizing  delays,  and  for  some 
river  channels,  minimizing  costs  associated  with 
locking  operations.  The  results  from  both  the  ice- 


tank  experiments  and  the  exercises  with  the 
numerical  model,  however,  indicate  that  sophisti¬ 
cated  transiting  schedules  are  not  called  for,  espe¬ 
cially  if  the  vessels  are  river  tows.  This  is  fortunate 
because  scheduling  of  transits  is  fraught  with 
numerous  practical  difficulties,  the  least  of  which 
is  the  unpredictable  nature  of  severe  weather. 

The  experiments  and  the  numerical  model 
indicate  that  a  basic  rule  for  minimizing  volume  of 
ice  grown  in  a  frequently  transited  channel  is  to 
minimize  the  total  number  of  icebreaking  transits. 
The  numerical  model  showed  further  that,  for  a 
fixed  number  of  transits  per  day  (or  a  fixed  se¬ 
quence  frequency /s),  varying  the  interval  between 
individual  transits  (or  transit  frequency  /t)  had  no 
significant  effect  on  the  volume  of  ice  grown  (Fig. 
72).  For  example,  if  it  is  required  that  a  channel  be 
transited  daily  by  four  vessels  (j f  =  4/day),  then 
there  is  negligible  difference  in  the  volume  of  ice 
grown  in  the  channel  whether  it  is  transited  in  a 
brief  sequence  during  which  vessels  pass  at  a 
frequency  of  one  per  two  hours  [/t  =  1  /  (2  hours)], 
with  no  transits  for  the  remainder  of  the  day,  or  if 
the  transits  are  equally  spread  over  an  entire  day  [/( 
=  1/(6  hours)]. 

Convoying  of  vessels  is  the  special  case  of  sched¬ 
uling  that  prevails  when/t= 0.  The  chief  advantage 
of  convoying  is  that  each  convoy  involves  only  one 
icebreaking  event.  Icebreaking  and  the  brunt  of 
transiting  resistance  are  borne  by  the  lead  vessel, 
which  could  be  used  especially  for  those  tasks. 
Furthermore,  because  it  reduces  the  number  of 
icebreaking  transits,  convoying  is  likely  to  reduce 
the  volume  of  ice  grown  in  a  navigation  channel. 
An  advantage  of  convoying  tows  is  that,  except  for 
the  lead  tow,  tows  within  a  convoy  would  likely 
experience  less  under-hull  accumulation  of  brash 
ice  or  broken  ice,  because  they  wc  old  not  be  tran¬ 
siting  intact  ice  covers  or  layers  of  brash  ice.  Con¬ 
voying,  however,  does  have  its  practical  draw¬ 
backs.  A  fundamental  one  .s  that  convoys  may  be 
difficult  to  form  and  maintain  if  vessels  arrive  at 
convoy-staging  locations  at  unpredictable  and 
irregular  intervals  and  if  vessels  then  have  diverse 
and  multiple  destinations.  A  further  drawback 
exists  for  river  channels  that  are  comparatively 
short  and  are  punctuated  with  locks,  such  as  many 
channels  comprising  the  inland  waterways  of 
North  America.  Erstwhile  convoys  transiting  these 
channels  may  be  disrupted  by  lock  passage,  as,  in 
some  extreme  situations,  the  time  required  for  a 
convoy  to  pass  through  a  lock  may  exceed  the  time 
required  to  transit  between  locks.  These  draw¬ 
backs  could  be  overcome  by  limiting  the  size  of 
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convoys  according  to  the  circumstances  of  the 
channel  to  be  transited.  For  example,  convoys 
formed  of  two  to  four  tows  may  be  effective  in 
reducing  ice  formation  in  relatively  short  chan¬ 
nels.  Generally,  however,  convoying  seems  more 
attractive  for  longer  channels  such  as  navigable 
rivers  of  northern  Soviet  Union,  where  convoying 
and  scheduling  of  transits  are  the  rule  rather  than 
the  exception  (Tronin  et  al.  1984). 

Because  the  more  severe  problem  created  by 
transiting  is  the  potential  for  thick  accumulations 
of  brash  ice  to  develop  at  specific  troublesome 
locations  (such  as  shown  in  Fig.  4)  in  navigation 
channels,  an  alternate  approach  is  to  use  mechani¬ 
cal  methods  for  controlling  brash-ice  accumula¬ 
tions.  This  approach  is  probably  more  likely  co 
minimize  the  volume  of  ice  grown  than  is  schedul¬ 
ing  of  transits.  Although  the  present  study  did  not 
directly  examine  mechanical  methods  of  improv¬ 
ing  transiting,  the  insights  gained  from  it  suggest 
three  general  methods:  vessel-based  techniques 
for  preventing  excessive  ridge  formation;  excava¬ 
tion  and  removal  of  brash  ice;  and  modification  of 
channel  alignment  and  cross  section,  especially  for 
relatively  shallow  or  narrow  channels.  All  three 
methods  require  further  investigation  to  ascertain 
their  effectiveness. 

Vessel-based  methods  could  be  developed  for 
controlling  ice-cover  formation,  breaking  ridges 
or  redistributing  accumulated  brash  ice.  Solve 
(1986),  for  example,  described  the  use  of  small 
boats  fitted  with  rakes  for  removing  brash  ice  from 
entrances  to  lock  approaches  on  the  Trollhatte 
Canal.  Mellor  et  al.  (1978)  and  Vance  (1980)  as¬ 
sessed  the  physical  and  economic  feasibilities  of 
various  methods  for  removing  brash  ice  from 
navigation  channels.  The  impetus  for  their  studies 
was  the  need  to  control  brash-ice  accumulations  in 
portions  of  theSt.  Marys  River.  The  most  attractive 
technique  that  both  studies  identified  involved  the 
use  of  a  vessel  for  mechanically  excavating  brash 
ice  from  channels.  Alterations  to  channel  align¬ 
ment  and  geometry,  which  may  be  required  for 
recurring  problem  sites,  are  a  last  resort.  Never¬ 
theless,  channel  modifications,  dredging  activities 
and  new  lock  approaches  and  river-training  works 
should  be  planned  with  consideration  of  their 
effects  on  brash-ice  accumulation. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  major  contributions  of  this  study  are  de¬ 
tailed  descriptions  of  ice  formation  in  frequently 


transited  navigation  channels  and  the  develop¬ 
ment  of  a  predictive  model.  A  practical  objective  of 
the  study  was  to  assess  the  merits  of  scheduling 
vessel  transits  so  as  to  control  or  mitigate  potential 
problems  caused  by  heavy  traffic.  The  study  in¬ 
volved  extensive  laboratory  experiments,  which 
were  conducted  using  an  ice  tank  and  model  ves¬ 
sels,  and  numerical  simulation  using  the  predic¬ 
tive  model.  Although  the  descriptions  of  ice  for¬ 
mation  apply  to  frequently  transited  channels  in 
general,  the  study  was  primarily  aimed  at  river 
channels  transited  by  tows  such  as  those  that  ply 
the  navigable  inland  waterways  of  North  Amer¬ 
ica.  The  principal  conclusions  derived  from  the 
study  are  summarized  here. 

A  conclusion  tentatively  drawn  from  the  dispa¬ 
rate  and  somewhat  sketchy  field  descriptions  of 
ice  formation  in  navigation  channels,  and  con¬ 
firmed  subsequently  during  the  ice-tank  experi¬ 
ments,  is  that  ice  formation  is  significantly  affected 
by  the  hull  form  and  length  of  transiting  vessels, 
the  channel  geometry  and  the  presence  of  signifi¬ 
cant  water  currents.  For  example,  quite  different 
patterns  of  ice  formation  may  evolve  in  river  chan¬ 
nels  that  are  plied  by  bluff-bowed,  flat-bottomed 
tows  or  have  significant  water  currents  than  occur 
in  harbor  and  coastal  channels  that  are  plied  by 
ship-form  hulls  and  are  subject  to  comparatively 
weak  currents.  The  essential  difference  in  the  pat¬ 
terns  is  that  transiting  tows  and  water  currents 
may  cause  broken  and  brash  ice  to  be  conveyed 
along  channels  and  may  consequently  lead  to 
markedly  nonuniform  ice  formation.  Two  general 
patterns  of  ice  formation  are  illustrated  in  Figures 
14  and  25.  They  illustrate  the  main  differences  in 
ice  formation  in  channels  transited  by  ship-form 
hulls  and  tows  that  move  ice,  respectively. 

Frequent  transiting  of  ice-covered  navigation 
channels  results  in  increased  volume  of  ice  growth, 
but  the  increase  is  perhaps  not  quite  as  dramatic  as 
might  be  expected.  The  results  from  both  the  ice- 
tank  experiments  (Fig.  60)  and  the  numerical  model 
(Fig.  71 )  indicate  that  intense  transiting,  involving 
frequencies  in  excess  of  about  two  transits  per  day, 
produces  about  two  to  three  times  the  volume  of 
ice  than  would  grow  over  the  same  area  of  channel 
opened  by  the  track  if  it  were  not  transited.  Lesser 
volumes  of  ice  grow  with  fewer  than  about  two 
transits  per  day. 

Transiting  breaks  ice  and  converts  it  to  brash 
ice,  which  accumulates  as  porous  layers  in  vessel 
tracks  and  ridges  formed  beneath  the  ice  cover 
bordering  tracks.  The  ice-tank  experiments  and 
the  numerical  modei  »how  that  intense  transiting 


93 


of  a  track  leads  to  layers  that  average  about  two  to 
three  times  the  thickness  of  the  ice  cover  riven  by 
the  track.  This  range  of  layer  thickness  corresponds 
to  ridge  thicknesses  that  are  three  to  four  times  the 
ice  cover  thickness. 

Two  essential  problems  are  created  by  frequent 
transiting.  One  problem  is  increased  volume  of  ice 
growth,  which  increases  the  resistance  to  transit¬ 
ing  vessels,  produces  more  ice  at  ice  jams  and 
exacerbates  difficulties  in  passing  ice  at  dams.  The 
other  problem,  which  is  potentially  more  severe, 
arises  because  brash  ice  may  accumulate  as  thick 
layers  along  vessel  tracks,  as  thicker  ridges  flank¬ 
ing  them,  and  for  channels  transited  by  tows,  as 
especially  thick  local  accumulations  termed  plugs. 
Both  problems  (but  particularly  the  latter  one) 
become  more  troublesome  for  shallower  channels 
with  depths  marginally  greater  than  the  vessel 
draft,  because  the  lack  of  clearance  prevents  brash 
ice  from  passing  beneath  the  vessels.  Ridges  and 
plugs  may  extend  over  the  entire  depth  of  shallow 
channels. 

The  ice-tank  experiments  and  exercises  with 
the  numerical  model  indicate  that  sophisticated 
scheduling  of  transits  to  minimize  ice  problems  is 
not  warranted.  A  general  conclusion  is  that  fewer 
ice-breaking  transits  result  in  less  ice  growth. 
Convoying,  therefore,  appears  attractive  because 
it  reduces  the  number  of  ice-breaking  transits. 
However,  an  approach  that  controls  ice  growth 
and  brash-ice  accumulation  at  troublesome  chan¬ 
nel  locations  is  more  likely  to  mitigate  problems 
incurred  by  frequent  transiting.  This  approach 
would  require  development  of  appropriate  me¬ 
chanical  means  of  redistributing  or  removing  brash 
ice. 

A  problem  peculiar  to  tow  transit  of  ice-cover 
channels  is  the  accumulation  of  broken  ice  and 
brash  ice  beneath  along  the  flat  bottoms  of  tows. 
The  experiments  show  that  accumulations  on  the 
order  of  15-20  times  the  ice-cover  or  brash-ice 
thickness  may  accumulate  beneath  tows.  For  the 
usual  range  of  tow  speeds,  the  accumulation  thick¬ 
ness  decreases  with  increasing  speed  for  tows 
transiting  ice  sheets.  However,  for  transits  through 
thin  covers  of  brash  ice,  the  accumulation  thick¬ 
ness  increases  with  increasing  speed .  The  accumu¬ 
lation  thickness  is  affected  by  tow  dimensions  and 
bow  geometry. 

As  a  logical  extension  of  the  present  study,  we 
recommend  that  sections  of  an  actual  tow  or  ship 
track  be  monitored  so  as  to  confirm  the  major 
conclusions  drawn  from  the  present  study.  Ideally 
such  field  monitoring  might  comprise  field  tests 


conducted  with  sections  of  navigable  channels 
being  transited  at  regular  frequencies.  A  further 
extension  would  be  to  survey  important  naviga¬ 
tion  channels  to  identify  locations  where  trouble¬ 
some  accumulations  of  brash-ice  form.  A  practical 
extension  would  then  be  to  investigate  effective 
methods  for  moving  and  "storing"  brash  ice  accu¬ 
mulated  at  these  channel  locations.  Among  the 
more  fundamental  aspects  of  ice  formation  to  be 
studied  further  is  ice  grown  through  porous  layers 
of  brash  ice. 
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APPENDIX  A.  BRASH-ICE  RESISTANCE  TO  TOW  TRANSIT 


Introduction 

As  part  of  their  study  on  the  icebreaking  performan  o  of  tows  on  the  Mississippi 
River,  Ashton  et  al.  (1973)  estimated  mean  values  of  resist  ,  nee  experienced  by  tows 
transiting  brash-ice  covers.  Their  resistance  data,  reproduced  in  Figure  Al,  show  a 
curious  if  not  controversial  trend:  the  total  resistance  generally  decreased  with 
increasing  tow  speed  and  brash-ice  thickness,  and  the  ice  resistance  (evaluated  as 
total  resistance  minus  open  water  resistance)  decreased  with  increasing  tow  speed. 
Ashton  et  al.  suggested  that  with  increasing  tow  speed,  less  ice  was  accumulated 
either  as  a  false  bow  or  beneath  its  hull.  Less  ice  accumulation  would  result  in  less 
resistance.  However,  the  trend  shown  by  their  data  is  counter  to  those  reported  for 
other  hulls  (e.g.  Ettema  et  al.  1985). 

A  brief  series  of  experiments  was  conducted  to  assess  the  levels  of  resistance  and 
power  requirements  for  tow  transit  of  navigation  channels  that  are  covered  with 
brash  ice. 

Experiments 

The  experiments  were  conducted  using  IIHR's  ice  tank  and  the  l:30-scale  tow 
model  described  earlier.  The  model  was  towed  at  speeds  ranging  from  0.05  to  0.75 
m/s  along  a  channel  containing  a  preformed  layer  of  brash  ice,  as  illustrated  in  Fig¬ 
ure  A2.  Two  configurations  of  tow  model  were  used.  For  one  the  towboat  pushed 
two  barges  aligned  in  parallel  (model  length  and  beam  =  3.05  and  0.71  m,  respec¬ 
tively).  For  the  other  the  towboat  pushed  two  barges  in  series  (model  length  and 
beam  =  4.88  and  0.36  m,  respectively).  For  both  configurations  the  model  draft  was 
92  mm.  The  tow  models  were  connected  by  a  towing  shaft  to  a  dynamometer 
mounted  on  the  towing  carriage.  A  guide  fitted  to  the  lead  barge(s)  ensured  that  the 
model  held  course  along  the  channel.  In  this  towing  harness  the  tow  model  was  free 
to  pitch  and  heave  as  well  as  undergo  slight  roll  and  yaw.  Towing  forces  were 
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a.  Hull  towed  through  a  preformed  layer  of  brash  ice 


b.  View  of  the  track  left  by  the  tow.  The  markers  indicate 
the  shear  deformation  of  the  layer. 


Figure  A2.  l:30-scale  tow  hull. 


measured  using  a  45-N  load  cell  fitted  to  the  towing  dynamometer. 

For  most  of  the  experiments  the  channel  width  was  2.13  m.  A  subset  of  experi¬ 
ments  was  performed  with  the  channel  width  at  0.84  m.  In  other  words,  the  channel 
was  1.2  and  3.0  times  the  beam  of  the  tow  model  configured  with  two  barges  in 
parallel,  and  2.4  and  6.0  for  the  series  configuration.  Since  the  water  depth  for  all  the 
experiments  was  maintained  at  1.07  m,  such  that  y  JD  =  12,  the  effects  of  shallow 
water  on  resistance  were  not  investigated.  The  tow  resistance  was  examined  for 
three  layer  thicknesses  h:  4  (single  layer),  50  and  1 00  mm.  Values  of  h/D  ranged  from 
0.05  to  1.14.  The  brash  ice  was  formed  by  fragmenting  4-mm-thick  sheets  of 
unseeded  urea.  The  average  dimension  of  the  brash  ice  was  1 9  mm,  with  a  geometric 
standard  deviation  of  o  =1.31.  The  brash  ice  was  confined  as  a  layer  over  a  channel 
formed  by  means  of  two  floating  side  panels,  which  are  visible  in  Figure  A2. 

For  each  experiment  the  brash-ice  layer  was  leveled  such  that  it  was  of  the 
requisite  average  thickness;  at  some  points  the  thickness  varied  by  as  much  as  4  mm, 
the  thickness  of  a  single  piece  of  brash  ice.  The  model  was  then  accelerated  t-  -  peed 
and  towed  at  a  constant  speed  for  about  12  m  along  the  ice  tank.  Each  transit  of  the 
model  through  the  channel  was  videotaped  to  record  the  manner  of  brash-ice 
movement  around  the  model.  Of  particular  interest  was  the  influence  of  model 
speed  on  the  formation  and  areal  concentration  of  open  water  in  the  track  left  by  the 
model. 

Results 

The  relationships  between  the  mean  total  resistance  R  and  the  hull  speed  V  for  the 
tow  in  parallel  and  series  configurations  are  shown  in  Figures  A3  and  A4,  respec¬ 
tively.  The  third  variable  in  these  figures  is  h/D.  The  general  trend  apparent  in  these 
figures  is  that  R  increases  with  increasing  V,  though,  for  the  model  configured  with 
tows  in  parallel,  R  decreased  mildly  as  V  increased  to  about  0.2  m2  (1.1  m2  at  full 
scale).  No  such  decrease  occurred  for  the  model  configured  with  the  two  barges  in 
series.  The  difference  in  the  trends  for  R  can  be  attributed  to  the  influence  of  false- 
bow  formation.  The  broader  parallel-barge  arrangement  created  a  larger  false  bow 
than  did  the  series-barge  arrangement.  Figure  A2b  illustrates  the  lateral  extent  to 
which  false-bow  formation  generates  shear  stresses  through  a  layer  of  brash  ke  and 
thereby  increases  layer  resistance.  The  extent  of  false-bow  formation  decreased  with 
increasing  V,  and  so  did  the  component  of  resistance  attributable  to  false-bow  for¬ 
mation. 

Figures  A3  and  A4  show  that  the  relationship  between  mean  ice  resistance  R;  (= 
R-R  )  and  V  is  similar  to  that  between  R  and  V.  On  the  whole,  R  increased,  albeit 
mildly,  with  increasing  V.  However,  it  did  not  decrease  with  increasing  V  as  it  did 
in  the  data  presented  by  Ashton  et  al. 

When  the  channel  was  narrowed  from  W /B  =  3.0  to  W/B  =  1.2,  the  tow  with  two 
barges  in  parallel  encountered  increased  resistance,  except  when  transiting  the 
thinnest  layer  ( h/D = 0.05)  at  creeping  speeds.  For  those  transits  the  width  of  the  false 
bow  was  limited  to  the  channel  width,  causing  less  ice  to  be  shoved  and  therefore 
resulting  in  diminished  R. 

Other  than  indicating  that  R  generally  increased  with  increasing  h/D  and  V, 
further  discussion  of  the  data  presented  in  Figures  A3  and  A4  is  complicated  by 
variations  in  the  processes  contributing  to  resistance. 


99 


Conclusions 

The  following  conclusions  concerning  the  relative  resistance  for  each  tow  con¬ 
figuration  can,  however,  be  made  from  the  data.  The  tow  comprising  two  parallel 
barges  generally  encountered  greater  resistance  than  did  the  tow  comprising  two 
barges  in  series.  For  transits  at  creeping  or  slow  speed,  the  difference  of  R  values 
diminished  with  increasing  h/D.  The  difference  in  R  values  also  decreased  with 
increasing  Vsuch  that,  at  the  highest  speed  tested,  the  resistances  are  of  comparable 
magnitude.  Essentially,  when  brash  ice  significantly  accumulated  beneath  the  hull 
or  formed  as  a  false  bow,  the  wider  hull  encountered  about  three  times  more 
resistance. 


APPENDIX  B:  PHOTOGRAPHS  OF  BRASH-ICE  SAMPLES 


Figure  Bl.  l:15-scale 
tow ;  f(  =  1/(2  hours). 


Figure  B2.  l:15-scale 
tow;  f(  =  1/(4  hours). 


Figure  B3. 1:1 5-scale 
tow;  ff  =  1 1(0.5  hour) 

andi  =  1/(7. 5  hours) 

n  1 
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Figure  B4.  l:15-scale 
tow;  f(  =  1/(8  hours). 


Figure  B5.  l:15-scale 
tow;  (f  =  1/(4  hours) 
and  yJD  =  2.78. 


Figure  B6.  l:15-scale 
tow;  f(  =  ?  1(2  hours) 
and  yJV-1.16. 
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Figure  B7.  l:15-scale 
wedge;  f(  =  1/(4  hours). 


Figure  B8.  l:15-scale 
wedge;  ff  =  1/(2  hours). 


Figure  B9  1:1 5-scale 
wedge;  =  l/(  0.5 
hour)  and  f  =1/(7. 5 

n  1 

hours). 
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Figure  BIO.  l:15-scale  wedge;  ft  =  1/(8  hours). 


Figure  Bll.  l:15-scale  wedge ;  f(  =  l/(0.5  hours). 
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APPENDIX  C:  VOLUMES  OF  ICE  GROWN  ALONG  EACH  TRACK  TRANSITED 
BY  THE  WEDGE  HULL  AND  THE  THICKNESS  OF  BRASH-ICE  ACCUMULA¬ 
TION  IN  EACH  TRACK,  BOTH  NORMALIZED  WITH  r\. 


APPENDIX  D.  PREDICTIONS  FROM  ASHTON'S  FORMULATION 


Ashton's  formulation  of  ice  formation  in  frequently  transited  navigation  channels 
is  exercised  here  to  compare  its  predictions  with  those  obtained  with  the  formulation 
presented  earlier.  For  Ashton's  formulation,  a  =  24  mm/(°C-day)  and  P  =  0.25. 
Figures  D1-D3  are  to  be  compared  with  Figures  69-71. 

For  transit  frequencies  less  than  about  1/day,  Ashton's  formulation  predicts 
values  of  q./q  that  are  comparable  with  those  predicted  using  the  authors'  formu¬ 
lation.  For  /t  exceeding  1  /day,  eq  10  predicts  even  increasing  values  oih.Jh.  The  pre¬ 
dictions  obtained  using  the  authors'  formulation  are  also  indicated  for  comparison 
on  Figure  D2. 


Figure  Dl.  Variation  of  q  .  with  time. 
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